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Abstract 


The  electric  glow  discharge  in  oxygen  was  experimen- 
tally studied  and  mathematically  modeled  for  the  parameter 
space  of  one  to  ten  torr  of  pressure  and  for  current  densi- 
ties up  to  35  milliamperes  per  square  centimeter.  Experi- 
mentally it  was  observed  that  the  oxygen  discharge  existed 
in  one  of  two  forms:  a high-electric-field  form  or  a low- 
electric-field  form.  The  high  field  form  was  present  at 
high  currents  and  low  pressures,  the  low  field  form  at  low 
currents  and  high  pressures.  The  high  field  form  exhibited 
stable  characteristics  while  the  low  field  form  was  observed 
to  be  nonstable,  exhibiting  periodic  or  aperiodic  oscilla- 
tory behavior.  The  periodic  oscillations  existed  only  at 
certain  frequencies  or  modes,  and  the  discharge  impedance 
was  a function  of  the  frequency.  The  optical  output  of  the 
oscillatory  form  of  the  discharge  was  highly  modulated  with 
short  duty  cycle  pulses.  The  oscillations  exhibited  a small 
amount  of  normal  dispersion.  The  point  of  transition  between 


forms  and  the  electrical  impedance  of  the  forms  was  unaf- 


fected by  the  addition  of  an  inert  gas,  argon,  to  the  dis- 
charge. Adding  an  electron-detaching  gas,  carbon  monoxide, 
increased  the  low  field  impedance  until  the  difference  be- 
tween the  forms  essentially  ceased  to  exist.  Mass  spectra 
measurements  indicated  that  the  atomic  oxygen  density  in- 
creased significantly  when  the  discharge  transitioned  from 


xviii 


the  low  to  high  field  forms.  Solutions  of  the  Boltzmann 
transport  equation  indicated  that  the  electron  energy  dis- 
tribution function  exhibited  a characteristic  intermediate 
between  Maxwellian  and  Druyvesteyn.  Chemistry  model  results 
showed  the  dominant  positive  ion  to  be  02  and  the  dominant 
negative  ion  to  be  0 . The  transition  from  the  stable  to 
unstable  forms  was  determined  to  be  related  to  two  criteria: 
the  dominance  of  the  derivative  of  the  attachment  coeffi- 
cient with  respect  to  E/N  over  that  of  the  ionization  coeffi 
cient  and  the  equivalence  of  the  negative  ion  and  electron 
number  densities.  The  former  was  satisfied  for  all  values 
of  E/N  below  67  TD.  The  later  condition  was  determined  to 

be  a function  of  pressure  and  accurately  predicted  the  point 
. * 1 

of  transition.  The  metastable  02(a  A)  was  important  in  the 
discharge  for  two  reasons.  Two  step  ionization  was  required 
in  addition  to  direct  ionization  to  match  the  experimental 
electron  density  measurements.  Also,  the  inclusion  of  the 
metastable  molecule  as  the  primary  electron  detacher  was  re- 
quired to  match  the  transition  points  determined  by  experi- 
ment with  those  determined  from  the  model. 
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DISCHARGE  PROCESSES  IN  THE  OXYGEN  PLASMA 


I.  Introduction 


This  chapter  begins  with  a statement  of  the  objectives 
of  the  study  and  then  describes  the  general  approach  that  was 
used  to  achieve  the  objectives.  Next,  in  order  to  provide  a 
background  for  the  reader,  the  general  characteristics  and 
criteria  of  a plasma  are  presented;  and  it  is  shown  that  the 
parameter  space  covered  by  the  study  met  the  listed  criteria. 
The  chapter  then  focuses  on  the  particular  type  of  plasma 
that  was  studied  - the  electrical  glow  discharge  positive 
column.  Finally  the  chapter  discusses  the  oxygen  positive 
column:  its  characteristics,  its  unique  properties,  and  the 

theories  that  have  been  developed  to  explain  them. 


Objectives 


>*■*  The  objectives  of  this  study  were  to  investigate  the 
plasma  existing  iri  the  positive  column  of  an  oxygen  glow  dis- 
charge and  to  determine  the  physical  and  chemical  processes 
involved. 

Recently  there  has  been  a renewed  interest  in  the  oxy- 
gen plasma,  an  interest  generated  by  the  development  of  the 
laser.  Oxygen  has  been  postulated  and  used  as  both  the  las- 
ing medium  and  the  pumping  medium  in  electrical  and  electri- 
cal-chemical lasers  -1+4)— 16 , Ref  ? • at ^-414 ; vi-ig) 
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V 

4:2243-2245,  Ref  5:834-8384-  Also,  there  has  been  a con- 
tinuing interest  in  efficient  production  of  ozone  (Ref  6:1-4). 
Effective  utilizations  of  oxygen  electric  discharges  require 
a knowledge  of  the  oxygen  plasma,  a plasma  which  differs  from 
other  plasmas  in  many  of  its  characteristics.  For  example, 
electrical  glow  discharges  in  oxygen  have  been  observed  to 
exist  in  two  different  forms:  a low  electrical  impedance 
form  which  exhibits  an  oscillatory  behavior  and  a high  elec- 
trical impedance  form  which  is  relatively  stable.  The  transi- 
tion between  the  two  forms  occurs  as  a function  of  current 
density  and  gas  pressure  .y  'The  characteristics  of  the  two 
forms  and  the  transition  phenomenon  have  been  studied  for 
many  years,  and  several  different  theories  have  been  devel- 
oped  - theories  which  have  not  been  entirely  consistent  with 
each  other.  Therefore  it  was  a purpose  of  this  study  to 
determine  a consistent  explanation  of  the  characteristics  of 
the  oxygen  discharge. 

Specifically  the  study  was  designed  to  determine: 

1.  A consistent  set  of  measurements  of  the  physical  proper- 
ties of  the  two  forms  of  the  oxygen  discharge. 

2.  Measurements  of  the  effects  on  the  discharge  of  adding 
inert  and  electron-detaching  gases. 

3.  The  energy  distribution  function  of  the  electrons. 

4.  The  relative  densities  of  the  various  charged  and  neutral 
species  present  in  the  discharge  and  the  dominant  chemi- 
cal and  physical  processes  associated  with  their  produc- 
tions and  losses. 
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5.  A resolution  of  the  role  that  the  metastable  molecule 

O2  (a-^A)  plays  in  the  ionization  and  electron-detachment 
processes . 

6.  An  explanation  of  the  two  forms  of  the  oxygen  discharge 
and  the  criteria  for  transitioning  between  them. 

The  study  of  the  oxygen  discharge  was  approached  from 
two  different  aspects:  an  experimental  investigation  and  a 
mathematical  modeling.  During  the  experimental  phase  the 
observable  characteristics  of  the  oxygen  positive  column  were 
investigated  for  a gas  pressure  range  of  one  to  ten  torr  and 
a current  density  range  of  .35  to  35  milliamperes  per  square 
centimeter.  Experimentally  determined  were  the  impedance 
characteristics,  gas  temperatures,  electric  field  shapes  and 
dispersion  characteristics  of  the  low  field  oscillations, 
electron  densities,  neutral  mass  spectra,  and  the  effects  of 
adding  other  selected  gases.  The  modeling  phase  consisted 
of  a determination  of  the  electron  energy  distributions  and 
electron-impact  pumping  rates  for  reduced  field  strengths 
(E/N)  from  one  to  one  hundred  Townsends*  and  a determination 
of  the  steady  state  charged  and  neutral  species  densities 
over  the  same  reduced  field  range  for  initial  molecular  den- 
sities from  3.2  x 10 ^ to  5.6  x 10^  molecules  per  cubic  cen- 
timeter. The  results  of  the  experimental  phase  were  used  to 
validate  the  mathematical  models;  the  results  of  computations 
using  the  mathematical  models  were  then  used  to  develop  the 


*E/N  is  the  Electric  Field  divided  by  the  density  of 
the  neutral  species.  One  Townsend  equals  10~17  volt-cm^. 
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theory  of  the  oxygen  plasma. 


Background 

General  Characteristics  of  a Plasma.  (Ref  7:2-8,  Ref 
8:1-8,  Ref  9:1-2,  Ref  10:205-246,  Ref  11:217-251).  The  term 
plasma  was  first  utilized  by  Irving  Langmuir  in  1928  to  des- 
cribe the  inner  region  of  an  electrical  discharge.  Since  • 
then  the  definition  has  been  broadened  to  describe  a state  of 
matter  in  which  a significant  number  of  the  atoms  and/or  mole- 
cules are  electrically  charged  or  ionized.  Characteristically 
in  a plasma  the  charged  particles  arrange  themselves  in  such 
a manner  as  to  effectively  shield  externally  and  internally 
generated  electrostatic  fields  within  a distance  of  a Debye 
length,  Ad  , and  to  achieve  local  charge  neutrality  over  re- 
gions greater  than  a Debye  length.  For  an  electron  density 
N at  temperature  Te , the  Debye  length  is  given  by: 


For  a confined  plasma  the  satisfaction  of  the  above  criteria 
requires  that  the  physical  dimensions  of  the  confining  system 
be  large  compared  to  the  Debye  length  and  that  there  be  a 
sufficient  number  of  electrons  within  a Debye  sphere*  to  pro- 
duce shielding.  For  a cylindrical  geometry  of  radious  RQ  the 
plasma  criteria  are: 


*Sphere  with  radius  \6 . 
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(2) 


A.cJ  R0 

Ne  « 1.9  xio^  7e3 

For  the  parameter  space  of  this  study  the  electron  tempera- 
ture  varied  between  2 x 10J  and  4 x 1(P  K;  the  tube  radius 

q 11 

was  .95  cm,  and  the  electron  density  varied  from  10  to  10 
cm~3.  This  resulted  in  a maximum  Debye  length  of  .05  cm. 

Thus,  the  criteria  of  equations  2 and  3 were  satisfied. 

The  observable  characteristics  of  a plasma  are  strong 
functions  of  the  species  that  constitute  the  plasma.  Colli- 
sions between  species  play  a particularly  important  role  in 
ionized  gases.  For  a weakly  ionized  gas  (such  as  used  in 
this  study)  the  dominant  collisions  are  those  between  charged 
and  neutral  species  and  those  between  different  neutral  spe- 
cies; the  collisions  between  charged  species  can  often  be 
neglected.  The  properties  of  the  species  and  their  collision 
processes  play  important  roles  in  the  plasma  characteristics. 
Since  molecular  gases  have  many  more  inelastic  electron-impact 
cross  sections  than  atomic  gases,  plasmas  in  the  molecular 
gases  have  very  different  characteristics  from  those  in  the 
atomic  gases.  Electronegative  gases  form  stable  negative 
ions  upon  electron  impact;  electropositive  gases  do  not.  For 
plasmas  in  electropositive  gases  the  only  charged  species 
that  must  be  considered  are  the  positive  ions  and  the  elec- 
trons; for  plasmas  in  electronegative  gases  the  negative  ions 
must  also  be  considered.  The  characteristics  of  plasmas  with 
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three  charged  species  may  be  expected  to  differ  from  those 
with  only  two  charged  species.  Because  of  these  factors  oxy- 
gen is  a very  interesting  gas  to  study:  it  is  a molecular 
gas;  it  forms  stable  negative  ions;  its  electron  collision 
cross  sections  have  been  generally  well  defined;  and  it  is 
easily  obtained  in  a highly  purified  state. 

The  majority  of  the  universe  exists  in  the  plasma  state. 
This  includes  the  stars,  which  are  almost  completely  ionized 
because  of  their  high  temperatures;  this  type  of  plasma  is 
an  example  of  an  equilibrium  plasma  in  which  the  ionization 
is  thermally  produced  and  the  temperatures  of  the  neutral  and 
charged  species  are  in  equilibrium.  The  equilibrium  plasma 
is  relatively  uncommon  in  laboratory  experiments.  The  usual 
laboratory  techniques  involve  some  nonequilibrium  process  by 
which  the  ionization  is  maintained  by  raising  some  of  the 
charged  species  to  a higher  temperature  than  the  neutrals. 

The  most  common  example  of  this  is  the  gas  discharge  in  which 
an  electrical  potential,  applied  across  a gap,  provides  the 
selected  excitation  of  the  charged  species.  If  the  majority 
of  the  ionization  is  by  direct  electron  impact  rather  than  by 
thermal  channels,  the  discharge  is  a glow  discharge.  The 
positive  column  of  a glow  discharge  is  a region  in  which  the 
physical  and  chemical  processes  are  generally  uniform  or  of 
a bounded  periodic  nature.  The  positive  column  interfaces 
on  the  cathode  end  of  the  discharge  with  the  Faraday  dark 
space  and  on  the  anode  end  with  the  anode  glow.  Because  the 
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processes  are  well  ordered*  in  the  positive  column,  it  pro- 
vides a convenient  medium  for  the  investigation  of  the  plasma 
characteristics . 

The  electrical  current  through  the  positive  column  is 
maintained  primarily  by  the  highly  mobile  electrons.  Suffi- 
cient ions  are  created  to  maintain  space  charge  neutrality.. 

In  the  uniform  column  the  electric  field  is  constant  and  of 
relatively  low  value  compared  to  other  regions  of  the  dis- 
charge. However,  the  electric  field  is  sufficient  to  raise 
the  electron  temperature  to  as  high  as  40,000°K,  while  the 
heavy  particles  remain  within  several  hundred  degrees  of  room 
temperature. 

The  characteristics  of  the  positive  column  vary  with 
the  constituents  of  the  gas,  its  pressure  and  temperature, 
the  current  density,  and  the  boundary  conditions.  For  exam- 
ple, the  light  emitted  and  electric  field  are  dependent  on 
the  diameter  of  the  confining  tube.  Smaller  diameters  in- 
crease diffusion  losses,  and  therefore  a higher  field  is  re- 
quired to  sustain  the  discharge.  This  results  in  a more 
luminous  column.  The  distance  between  the  electrodes  essen- 
tially determines  the  length  of  the  column  - as  the  tube 
length  is  increased,  there  is  an  equivalent  increase  in  the 
length  of  the  positive  column;  the  other  regions  of  the 


*Well  ordered  in  the  sense  that  the  electric  field  and 
charged  particle  densities  are  relatively  constant  in  the 
axial  dimension.  Refer  to  Figure  1 for  a graphical  presenta- 
tion of  these  parameters  in  the  glow  discharge. 
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discharge  remain  constant.  As  the  pressure  of  the  discharge 
is  increased,  the  effects  of  collision  processes  become  more 
important  and  eventually  dominate  over  the  diffusion  losses. 
Gas  heating  leads  to  radial  density  gradients  and  higher  re- 
duced electric  fields,  E/N,  in  the  center  of  the  column  than 
at  the  boundary.  This  is  evidenced  by  the  apparent  radial 
constriction  of  the  column  as  the  pressure  or  current  is  in- 
creased. The  electric  field  in  molecular  gas  discharges  is 
higher  than  in  atomic  gas  discharges  by  an  order  of  magnitude. 
This  is  caused  by  the  greater  number  of  inelastic  loss  pro- 
cesses of  the  molecular  gas  which  dissipate  the  electron 
energy  into  nonionizing  channels. 

The  positive  column  can  be  either  uniform  or  oscilla- 
tory in  its  properties.  Oscillations  appear  as  alternate 
dark  and  bright  regions  which  may  be  either  stationary  or  may 
move  along  the  column.  Often  their  motion  is  so  rapid  as  to 
result  in  the  appearance  of  a uniform  column. 

General  Characteristics  of  the  Positive  Column  in  Oxy- 
gen. The  electrical  and  optical  characteristics  of  the  oxy- 
gen positive  column  are  unique;  they  occupy  an  intermediate 
place  between  the  normal  electropositive  discharges,  such  as 
nitrogen  or  the  rare  gases , and  the  strongly  electronegative 
discharges  exemplified  by  iodine  (Ref  12:519).  This  fact  has 
caused  the  oxygen  discharge  to  be  studied  by  many  investiga- 
tors through  the  years.  However,  only  recently  have  the  ex- 
perimental and  theoretical  techniques  been  available  to  fully 
define  and  explain  the  discharge  characteristics. 
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In  1927,  Guntherschulze  observed  that  the  oxygen  posi- 
tive column  existed  in  two  forms:  one  with  a high  voltage 
gradient  and  a greyish  green  optical  output  and  another  which 
had  a low  voltage  gradient  and  was  violet  in  color.  The  exis- 
tence of  one  or  the  other  form  was  noted  to  be  a function  of 
the  pressure  and  current  (Ref  13:763).  In  1951,  Seeliger  and 
Wichman  defined  the  High  field  form  as  the  H form  and  the  low 
field  form  as  the  T form*.  Their  impression  of  the  color 
difference  was  that  the  optical  emissions  from  the  H form 
were  blue  with  a light  violet  core;  whereas  in  the  T form,  a 
grey  blue  light  was  emitted  (Ref  14:235-244). 

Other  investigators  further  defined  the  electrical  and 
optical  characteristics  of  the  column.  It  was  noted  that  the 
H form  existed  at  higher  currents  and  lower  pressures  than 
the  T form.  If  the  pressure  was  held  constant,  an  increase 
in  current  would  cause  the  column  to  transition  from  the  T 
to  H forms.  Holding  current  constant  and  increasing  pressure 
would  cause  it  to  switch  from  the  H to  T forms . Pekarek  pub- 
lished a paper  in  which  the  existence  regions  for  the  two 
forms  were  defined  (Ref  15:750).  Figure  2 depicts  these  re- 
gions. The  occurances  of  the  T and  H forms  were  not  neces- 
sarily exclusive;  they  have  both  been  observed  in  the  column 
at  the  same  time.  When  both  forms  were  present,  they  were 
axially  separated  with  the  H portion  always  situated  on  the 


*H  and  T had  their  origins  in  the  German  words  for  High 
and  Low:  Hoch  and  Tief. 
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anode  side  of  the  column  and  the  T portion  on  the  cathode 
side.  As  these  parameters  were  swept  through  the  transition 
region  from  the  H form  to  the  T form,  it  was  observed  that 
the  T form  appeared  at  the  cathode  face  of  the  positive  col- 
umn and  progressed  to  the  anode.  Conversely  when  the  para- 
meters were  swept  from  their  T form  to  H form  values,  the  H 
form  started  at  the  anode  and  progressed  to  the  cathode. 

Also,  the  transition  between  the  forms  exhibited  a hysteresis 
effect.  At  any  position  in  the  positive  column,  if  either 
the  discharge  current  or  gas  pressure  was  swept,  the  transi- 
tion between  the  forms  was  dependent  on  the  direction  of  the 
sweep  (Ref  16:305-317,  Ref  17:3.1.5.1).  This  hysteresis  ef- 
fect is  shown  in  Figure  3. 

Of  the  two  forms  it  is  the  T form  which  is  considered 
to  differ  from  other  molecular  gases;  the  voltage  gradient 
in  this  form  is  anomalously  low  when  compared  with  other 
molecular  gases  (Ref  17:3.1.5.1).  Practically  all  of  the  in- 
vestigators have  noted  that  the  T form  is  oscillatory,  with 
the  oscillations  traveling  from  the  cathode  to  the  anode. 
These  oscillations  have  been  noted  as  either  being  periodic 
or  aperiodic  depending  on  the  discharge  parameters  (Ref  18: 
11-13) . The  H form  has  generally  been  noted  as  homogeneous 
with  a gradient  typical  of  other  molecular  gases;  this  has 
led  to  it  being  described  as  the  normal  form  of  the  oxygen 
discharge.  Sabadil,  however,  has  noted  that  under  certain 
conditions  the  H form  can  also  be  oscillatory  with  low  am- 
plitude moving  or  stationary  strata  (Ref  16:13-18). 
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Species  in  the  Oxygen  Positive  Column.  The  concentra- 
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tion  of  atomic  oxygen  and  ozone  have  been  observed  to  be 
functions  of  the  form  of  the  discharge.  Ozone  has  been 
claimed  only  to  exist  in  the  T form  and  not  in  the  H form 
(Ref  14:235-244,  Ref  19:1).  On  the  other  hand  the  concentra- 
tion of  atomic  oxygen  has  been  determined  to  be  much  greater 
in  the  H form  than  in  the  T form.  This  was  initially  deter- 
mined by  analysis  of  the  optical  output  of  the  column;  the 
molecular  band  structure  dominated  in  the  T form  and  the  ato- 
mic line  structure  in  the  H form  (Ref  14:235-244,  Ref  15:749). 
Theoretical  estimates  of  up  to  8%  dissociation  have  been  made 
(Ref  12:520).  The  major  source  of  atomic  oxygen  has  been 
postulated  to  be  from  electron  impact  excitation  of  the 
ground  state  oxygen  molecule  to  the  repulsive  vibrational 
levels  of  the  B3Xuu  electronic  state: 

e+ o2[X3I“>v<a)  -»  e+ o2  (B3r^,v> it) ->  e+ 20 

The  excitation  to  the  dissociative  vibrational  levels  is 
caused  by  the  offset  between  the  potential  energy  curves  of 
these  two  electronic  states  and  the  fact  that  it  is  a sudden 
collision,  and  consequently  the  Franck-Condon  principle  ap- 
plies (Ref  20:510-511,  Ref  21:11).  Figure  4 is  the  term  dia- 
gram for  oxygen  and  is  included  to  illustrate  the  potential 
energy  offsets.  Depending  on  the  pressure,  atomic  oxygen  was 
postulated  to  be  lost  by  diffusion  to  the  walls,  where  recom- 
bination takes  place,  or  by  collisional  processes.  The  col- 
lisional  processes  were  favored  due  to  the  low  recombination 
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Figure  4.  Term  Diagram  for  Oxygen  (From  Ref.  79) 
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rate  at  the  wall*  and  the  resulting  flat  radial  profile  (Ref 
22:3.1.5.9,  Ref  23:235-251). 

Because  of  the  electronegative  nature  of  oxygen,  the 
oxygen  plasma  is  a three  component  plasma  consisting  of  three 
types  of  charge  carriers:  electrons,  positive  ions,  and  nega- 
tive ions.  In  the  glow  discharge  the  predominant  positive 
ion  was  reported  to  be  02+  (Ref  24:820,  Ref  12:520,  Ref  19: 
8-9).  At  PRq  = .18,  P = .04  torr  and  j = 60  uA/cm^ , Thompson 
reported  a molecular  ion  to  atomic  ion  ratio  of  70  to  1 (Ref 
12:519).  At  a pressure  of  2.5  torr,  the  ratio  of  0^  to  0+  de- 
pended on  the  discharge  form?  in  the  H form,  0^  overwhelm- 
ingly dominated;  and  in  the  T form,  the  ratio  was  approxi- 
mately 2 to  1 (Ref  19:8-9).  The  predominant  negative  ion 
was  determined  to  be  0“  (Ref  25:1068,  Ref  26:299-308,  Ref  24: 
820).  For  the  .04  torr  observation  referenced  previously, 
Thompson  reported  a ratio  of  10  to  1 for  the  0~  to  O2-  ion 
(Ref  12:519).  The  main  production  mechanism  for  the  0-  ion 
was  postulated  to  be  the  dissociative  attachment  process: 
e + O2  “♦  0~  + 0 (Ref  24:820,  Ref  23:235-251).  In  a three 
component  plasma  with  a large  concentration  of  negative  ions, 
ambipolar  diffusion  was  considered  to  be  a weak  process;  and 
the  radial  electric  field  was  reported  to  be  very  low.  A 
weak  ambipolar  diffusion  process  reduces  the  radial  loss  to 
the  walls  of  the  negative  ions  and  requires  that  the  dominant 
loss  be  collisional  (Ref  24:818,  Ref  22:3.1.5.9,  Ref  23: 

*Estimated  at  10-^  for  pyrex  (Ref  23:234-251). 
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235-251) . The  associative  detachment  process  0"  + 0 02  + e 

was  suggested  as  the  major  loss  mechanism  for  the  atomic  ne- 
gative ion  (Ref  12:522,  Ref  27:53-66,  Ref  24:821,  Ref  22:3. 
1.5.9).  The  molecular  negative  ion  O2-  was  claimed  to  be 
primarily  formed  by  the  charge  exchange  process  with  0~ : 0"  + 

0 2 02_  + 0 (Ref  12:522).  Under  certain  conditions  in  the 

discharge,  the  negative  ion  concentration  has  been  determined 
to  dominate  the  electron  concentration  by  up  to  20  to  1 (Ref 
24:818) . 

There  are  two  low-lying  metastable  states  of  molecular 
oxygen:  the  02  (a*A)  and  02  (b1Z)  states.  The  role  that 

these  states  play  in  the  characteristics  of  the  oxygen  plasma 
has  been  unresolved.  There  are  investigators  who  postulated 
that  ionization  from  these  metastable  levels  must  play  an  im- 
portant and  necessary  role  in  the  oxygen  discharge  in  order 
to  match  experimentally  observed  ionization  rates  (Ref  17: 

3. 1.5.1,  Ref  19:2).  Other  investigators  discounted  stepwize 
ionization  as  insignificant  when  compared  to  direct  processes, 
since  the  energy  difference  between  the  ground  and  metastable 
states  are  small  compared  to  the  ionization  energy  (Ref  12: 
520,  Ref  27:53-66) . 

Oscillatory  Behavior  in  the  Oxygen  Discharge.  The  os- 
cillations in  the  T form  of  the  oxygen  discharge  have  been 
claimed  to  differ  from  those  found  in  N2,  H2,  or  noble  gas 
discharges.  The  oxygen  oscillations  are  very  highly  modulated 
and  have  been  observed  to  exist  in  the  T portion  of  the  col- 
umn when  both  forms  are  present  (Ref  15:750-751,  Ref  16: 
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305-317,  Ref  26:299-308).  The  oscillations  in  other  gases 
are  typically  less  modulated  and  appear  in  the  total  column. 
The  oxygen  oscillations  were  observed  to  be  forward  waves 
moving  from  cathode  to  anode  (Ref  16:305-317)  with  a velocity 
approximately  that  of  the  negative  ion  drift  (Ref  23:235-251). 
This  velocity  was  on  the  order  10^  to  10^  cm/sec,  and  the 
frequency  was  between  10  KHz  and  1 MHz  (Ref  18:11,  Ref  29: 

85,  Ref  15:751).  The  oscillations  have  been  reported  as  dis- 
persionless, i.e.,  the  velocity  was  independent  of  frequency 
(Ref  16:305-317) . 

The  oscillation  frequency  in  the  T form  was  found  to 
vary  with  the  parameters  of  the  discharge  (Ref  28:409-420). 

As  the  current,  the  length  of  the  column,  the  pressure,  the 
anode  diameter,  or  the  tube  diameter  was  increased,  the  fre- 
quency decreased  (Ref  16:305-317,  Ref  28:409-420).  The  am- 
plitude of  the  oscillations  was  observed  to  be  inversely  re- 
lated to  the  frequency  (Ref  28:409-420). 

The  oscillations  in  the  T form  are  not  always  periodic; 
for  certain  conditions  of  current  and  pressure,  they  have 
been  found  to  be  aperiodic  or  noiselike.  The  aperiodic  T 
form  has  been  observed  to  have  a slightly  lower  gradient  than 
the  periodic  one  (Ref  18:12).  It  was  possible  to  use  exter- 
nal modulation  of  the  column  to  synchronize  the  column  to 
one  of  its  periodic  modes  of  oscillation  (Ref  26:299-308). 

The  oscillations  in  the  H form  were  noted  to  be  similar 
to  those  found  in  N2 , H2,  and  noble  gas  discharges  (Ref  23: 
235-251) . The  H form  oscillations  have  been  observed  to  move 
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in  the  forward  direction,  backward  direction,  or  to  remain 
stationary  (Ref  16:305-317).  Sabadil’s  observations  were  on 
discharges  with  pressures  in  the  .1  to  .35  torr  range.  Other 
investigators,  whose  investigations  were  normally  in  higher 
pressure  discharges,  failed  to  determine  any  oscillations  in 
the  H form. 

Theories  Regarding  the  Two  Forms  of  the  Oxygen  Dis- 
charge . The  characteristics  of  the  H form  of  the  oxygen  dis- 
charge were  observed  to  be  very  similar  to  other  molecular 
gases;  therefore,  it  was  generally  acknowledged  that  this  is 
the  normal  form  of  the  discharge,  and  it  generally  exhibits 
two-component-plasma  characteristics.  It  is  the  T form  that 
was  listed  as  abnormal,  and  several  theories  have  been  devel- 
oped to  explain  the  characteristics  of  this  form  and  the  rea- 
son it  exists  in  only  a portion  of  the  pressure-current  para- 
meter space. 

In  a 1966  article  Sabadil  noted  that  the  T form  oscil- 
lations showed  no  dispersion  and  that  the  velocity  of  the  os- 
cillations was  independent  of  the  discharge  current  and  tube 
diameter  and  that  there  was  no  indication  of  damping  by  ion- 
neutral  collisions.  This  led  him  to  postulate  that  the  os- 
cillations were  ionic  sound  waves  resulting  from  an  ioniza- 
tion-recombination instability.  Based  on  previously  pub- 
lished data  on  the  recombination  rates  and  ion  densities,  he 
was  able  to  calculate  a positive  ion-negative  ion  recombina- 
tion time  constant  which  agreed  with  his  T oscillation  period 

» 

within  an  order  of  magnitude  (Ref  16:305-317).  In  a follow-on 
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paper  in  1968,  Sabadil  noted  that  the  coefficient  for  the  re- 
combination process  between  the  negative  and  positive  ions 
in  oxygen  was  inversely  related  to  the  gas  temperature 
(Kr  ~ Tg“5/2) . Experiments  with  various  oxygen  discharges 
at  temperatures  between  -170°C  and  +300°C  indicated  very 
slight  changes  in  the  oscillation  period  over  the  temperature 
range.  This  led  Sabadil  to  abandon  the  ionization-recombin- 
ation theory  for  the  oscillations  (Ref  26:299-308). 

In  a 1968  paper  Sabadil  noted  the  similarity  between 
the  Funn-Effect  instabilities  that  occur  in  gallium  arsenide 
semiconductors  and  the  oscillations  in  the  T form  of  the  oxy- 
gen discharge.  In  the  semiconductor,  differential  negative 
mobility*  arises  due  to  the  transfer  of  electrons  from  a high 
mobility  central  valley  in  the  conduction  band  to  satellite 
valleys  which  are  at  a higher  energy  but  where  the  mobility 
is  lower.  For  the  oxygen  discharge  the  atomic  negative  ion 
0~  was  suggested  as  the  low  mobility  charge  carrier  and  the 
electron  as  the  high  mobility  carrier.  To  exhibit  this  ef- 
fect the  production  of  low  mobility  carriers  must  depend  on 
some  process  which  has  an  energy  threshold.  For  GaAs  this 
threshold  is  the  energy  difference  between  the  central  and 
satellite  valleys  of  the  conduction  band;  for  oxygen  this 
threshold  is  the  minimum  energy  associated  with  the  dissocia- 
tive attachment  process  - approximately  4.4  ev.  For  both  the 
semiconductor  and  the  oxygen  discharge,  under  certain 


*As  a function  of  the  electric  field 
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conditions  the  charge  carrier  distribution  becomes  unstable; 
and  there  are  established  high  and  low  field  domains  which 
propagate  from  cathode  to  anode.  The  Sabadil  paper  proposed 
a qualitative  model  of  transitioning  from  the  low  field  to 
the  high  field  domain  but  failed  to  describe  a mechanism  for 
the  converse  (Ref  26:299-308). 

In  a 1974  report  Nighan  and  Wiegand  developed  a theore- 
tical model  of  the  causes  and  conditions  for  instabilities 
in  molecular  discharges,  particularly  discharges  in  which 
negative  ion  processes  play  a significant  role.  Their  analy- 
sis indicated  that  the  oscillations  in  the  T form  of  the  02 
discharge  were  a manifestation  of  an  ionization  instability 
caused  by  the  dissociative  attachment  process.  The  occur- 
rence of  this  instability  was  stated  to  be  a function  of  the 
electron  temperature  dependence  of  the  attachment  coefficient, 
the  electron  temperature  dependence  of  the  ionization  coeffi- 
cent,  and  the  relative  magnitude  of  the  negative  ions  and 
electrons.  Specifically,  a necessary  condition  for  the  in- 
stability to  occur  is  that  the  following  inequality  be  true: 


STe  ^ STe 


(4) 


where  Ka  is  the  attachment  coefficient  for  the  process: 
e + 02  ->  0~  + 0,  Ki  is  the  direct  ionization  coefficient; 
and  Te  is  the  electron  temperature.  The  dominance  of  the  at- 
tachment temperature  dependence  over  that  of  ionization  al- 
lows for  an  inverse  relationship  between  disturbances  in 
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electron  temperature  and  electron  density,  thereby  providing 
a positive  feedback  path  in  the  electron  kinetics.  A suffi- 
cient condition  for  the  oscillations  to  occur  is  that  the  ne- 
gative ion  concentration  must  be  significant  compared  to  the 
electron  concentration.  This  requires  that  detachment  pro- 
cesses play  a minor  role.  The  transition  from  the  T form  to 
the  H form  was  postulated  to  occur  when  the  detaching  species 
concentration  increases  to  the  point  where  the  detachment 
process  dominates  the  negative  ion  balance  equations.  For 
the  pure  oxygen  discharge  the  dominant  detaching  species  was 
suggested  to  be  either  the  neutral  atom  0 or  the  metastable 
neutral  molecule  02  (a1*)  - The  buildup  of  another  stable 
negative  ion  O3-  by  the  process:  0~  + 02  + M — * O3”  + M 

was  considered  as  the  method  of  delaying  the  transition  from 
the  T to  H forms.  These  two  trends  tend  to  explain  the  cur- 
rent and  pressure  dependence  of  the  transition  between  the 
forms  of  the  oxygen  discharge  (Ref  30:1-40). 
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II.  GENERAL  EXPERIMENTAL  SETUP 


The  general  setup  described  in  this  chapter  was  uti- 
lized throughout  the  experimental  phase  of  this  study.  With 
one  exception,  only  minor  changes  were  made  in  progressing 
from  one  experiment  to  another.  The  one  exception  was  the 
discharge  tube;  a different  design  was  utilized  for  the 
mass  spectra  experiments.  The  general  setup  was  designed 
to  provide  the  necessities  of  the  experiments;  controllable 
pressure,  gas  flow,  gas  mixture,  D.C.  electrical  power  with 
variable  potential  and  current,  and  cooling  for  the  discharge 
tube.  The  following  sections  describe  the  discharge  tube 
design  and  the  systems  employed  in  the  general  experimental 
setup. 

Discharge  Tube  (Figure  5) 

The  discharge  tube  was  constructed  entirely  of  quartz. 
This  material  was  chosen  over  glass  due  to  its  high  soften- 
ing temperature  and  low  electrical  dissipation  factor  at 
microwave  frequencies;  the  latter  characteristic  was  required 
for  the  electron  density  measurements.  With  the  exception 
of  the  anode  region,  the  tube  was  jacketed  for  cooling  pur- 
poses. The  inside  diameter  of  the  inside  tube  was  19mm  with 
the  exception  of  the  cathode  region.  In  this  region  the 
diameter  was  increased  to  25mm  to  accommodate  the  large  cath- 
ode. The  length  between  the  column  edges  of  the  anode  and 
cathode  was  505mm.  Ground  taper  joints  were  installed  every 
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Figure  5. 


DIMENSIONS  IN  mm 


100mm  along  the  column  to  allow  probes  to  be  inserted  into 
the  tube.  During  the  majority  of  the  experimentation,  elec- 
tric field  probes  were  inserted  into  the  taper  joints.  These 
probes  consisted  of  .25mm  diameter  platinum  wire  encapsulated 
in  glass,  except  for  2mm  at  the  end  near  the  center  of  the 
tube.  For  the  plasma  temperature  experiment,  a thermocouple 
was  inserted  through  one  of  the  taper  joints. 

The  electrodes,  anode  and  cathode,  were  constructed 
of  oxygen-free  copper  to  minimize  water  vapors.  The  lead 
pins  were  constructed  of  Kovar.  An  epoxy  substrate  was  used 
to  mechanically  support  and  to  seal  the  pins  at  the  points 
where  they  passed  through  the  glass.  During  the  break-in 
period  for  the  tube,  these  substrate  joints  had  to  be  ther- 
mally cycled  several  times  before  they  formed  good  mechani- 
cal and  vacuum  seals.  The  electrodes  were  broken-in  by  a 
gradual  increase  in  discharge  current  to  slowly  create  an 
even  oxidation  of  the  surfaces.  Too  fast  a break-in  caused 
localized  heavy  oxidation  which  resulted  in  surface  flaking 
and  surface  unevenness.  In  order  to  avoid  flaking  at  high 
currents, the  cathode  design  had  to  be  as  massive  as  the  tube 
design  would  allow.  An  80mm  long  by  20mm  diameter  hollow 
cylinder  was  experimentally  determined  to  be  sufficient  for 
currents  to  lOOma.  The  anode  size  was  relatively  insensi- 
tive to  discharge  current.  Therefore  a smaller  design,  40mm 
long  by  16mm  diameter  hollow  cylinder,  was  used  for  this 
electrode. 
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Gas  System  (Figure  6) 


Throughout  the  experimentation  the  primary  gas  used 
was  oxygen.  Other  gases  - helium,  argon,  or  carbon  monox- 
ide - were  used  to  initiate  the  discharge  and  as  additives 
to  the  oxygen.  Power  supply  limitations  precluded  discharge 
initiation  in  pure  oxygen.  All  the  gases  used  were  ultra 
high  purity  and  were  dispensed  from  size  A bottles. 

The  portion  of  the  gas  system  upstream  of  the  dis- 
charge tube  was  constructed  of  stainless  steel  components, 
with  the  exception  of  a short  section  of  tygon  tubing  con- 
necting the  mixing  chamber  to  the  tube  itself.  The  portion 
of  the  system  between  the  pressure  regulators  and  the  flow 
controlling  needle  valves  was  maintained  at  a small  positive 
pressure  relative  to  ambient  to  minimize  the  effect  of  small 
leaks  should  they  occur.  High  and  low  range  flowmeters  in 
the  oxygen  and  additive  gas  lines  allowed  mixing  to  one  part 
in  five  hundred  and  measurement  of  the  gas  flow.  The  flow 
and  mixtures  were  controlled  by  needle  valves  in  each  line. 
Mixing  occurred  in  a chamber  filled  with  small  diameter 
glass  beads. 

Downstream  of  the  discharge  tube,  coarse  and  fine 
needle-valves  were  used  to  precisely  control  system  pressure. 
A large  capacity,  17.7  CFM  Duo-Seal  #1397,  mechanical  vacuum 
pump  was  used  to  maintain  system  vacuum.  A 6.5  Cu.Ft.  glass- 
lined  water  tank  was  used  as  a ballast  for  the  system  to 
smooth  pump  pulses. 

System  pressure  was  measured  near  the  anode  end  of 
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the  discharge  tube  by  a Baratron  capacitance  manometer  sy- 
stem. Both  absolute  and  relative  100  torr  heads  were  uti- 
lized during  different  phases  of  the  experimentation. 

The  overall  gas  system  was  capable  of  operating  at  a 
minimum  pressure  of  one  micron  with  no  flow.  During  opera- 
tions with  a discharge  the  pressure  was  adjusted  to  between 
one  and  ten  torr  in  the  tube.  A positive  flow  was  main- 
tained throughout  all  experiments  to  minimize  the  buildup 
of  contaminants  in  the  discharge  volume  that  could  have  re- 
sulted from  small  leaks  or  from  out-gassing.  Flows  ranged 
between  10  and  500  SCCM. 

Electrical  System  (Figure  7) 

The  electrical  power  for  the  discharge  was  provided 
by  a Spellman  Model  HP5PN500SRX  power  supply.  This  supply 
was  capable  of  providing  a regulated  voltage  between  zero 
and  five  kilovolts  D.C.  or  a regulated  current  between  zero 
and  five  hundred  m.a.  The  current-regulating  circuit  of  the 
Spellman  was  not  capable  of  responding  to  the  rapid  change 
of  impedance  of  the  oxygen  discharge  during  transition  be- 
tween the  high  and  low  field  forms.  Therefore,  a custom- 
built,  fast-reaction,  current-regulator  was  placed  in  ser- 
ies between  the  Spellman  power  supply  and  the  discharge; 
and  the  Spellman  power  supply  was  operated  in  the  constant 
voltage  mode.  The  current  regulator  was  capable  of  being 
controlled  manually,  by  a motor  driven  sweep,  or  by  an  elec- 
tronic signal  of  arbitrary  waveform.  Pulses  were  used 
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during  one  experiment  to  synchronize  the  oscillations.  The 
utilization  of  a good  current  regulator  was  key  to  removing 
the  external  circuit  effects  from  the  discharge  characteris- 
tics. The  custom  regulator  used  in  this  study  outperformed 
any  commercially  available  at  the  time. 

The  current  through  the  discharge  tube  and  the  poten- 
tial across  the  tube  were  monitored  on  digital  voltmeters 
and  recorded  on  a Hewlett  Packard  7004A  x-y  recorder.  The 
electric  field  was  determined  from  the  potential  difference 
between  probes  inserted  in  two  ports  of  the  discharge  tube. 
Two  systems  for  measuring  the  electric  field  were  utilized 
during  various  phases  of  the  experimentation.  A laboratory- 
constructed  differential  amplifier  system  with  an  input  im- 
pedance in  excess  of  10^®  ohms  was  used  for  D.C.  and  low  fre 
quency  measurements.  For  high  frequency  measurements  of  the 
oscillatory  electric  fields,  frequency  compensated  Tektronix 
probes  with  input  impedances  of  10  ohms  were  used  in  con- 
junction with  a Tektronix  1A7A  differential  amplifier.  The 
outputs  of  these  systems  were  variously  connected  to  a digi- 
tal voltmeter,  a Hewlett-Packard  7004B  x-y  recorder,  a Tek- 
tronix 555  oscilloscope,  and/or  a Hewlett-Packard  8556A 
spectrum  analyzer. 

Cooling  System  (Figure  8) 

It  was  observed  early  in  the  experimentation  that  the 
oxygen  discharge  ran  very  hot;  a significant  portion  of  the 
input  energy  was  converted  to  heat.  In  order  to  insure 
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experimental  repeatability  and  to  protect  the  discharge  tube, 
a cooling  system  was  utilized  to  maintain  the  discharge  tube 
wall  temperature  near  that  of  tap  water.  Carbon  tetrachlor- 
ide was  utilized  as  the  coolant  due  to  its  low  dissipation 
factor  at  microwave  frequencies;  this  was  important  for  the 
electron  density  experiment. 

The  coolant  was  pumped  through  the  tube  at  a pressure 
of  approximately  five  to  ten  psi.  After  passing  through  the 
tube,  the  heat  was  removed  from  the  carbon-tetrachloride  in 
a water-cooled  heat  exchanger.  The  system  was  so  efficient 
that  the  outside  wall  of  the  outer  jacket  never  exceeded  a 
temperature  of  18°C,  as  measured  by  a platinum  resistance 
probe  taped  to  the  tube. 
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III.  Impedance  Characteristics 


This  chapter  is  the  first  in  which  experimental  results 
are  presented.  It  covers  the  measurements  of  the  electrical 
impedance  of  the  oxygen  discharges  both  with  and  without  ad- 
ditive gases.  Electrical  impedance  experiments  are  rela- 

» 

tively  easy  to  perform,  and  they  tell  a lot  about  the  over- 
all characteristics  of  a particular  plasma.  However,  it  is 
very  difficult  to  determine  from  these  results  the  explana- 
tions of  the  characteristics.  In  this  chapter  both  quanti- 
tative and  qualitative  results  are  presented.  Some  of  these 
are  explained  in  future  chapters  after  further  experimenta- 
tion or  modeling. 

The  chapter  starts  with  a brief  explanation  of  the 
factors  that  determine  the  impedance  of  a discharge.  This 
is  followed  by  a description  of  the  experimental  techniques 
used.  Then  the  experimental  results  are  presented.  Speci- 
fically, the  discharge  potential  and  electric  field  values 
are  depicted  as  functions  of  current  and  pressure  for  pure 
oxygen,  oxygen  plus  argon,  and  oxygen  plus  carbon  monoxide. 
Also,  mappings  of  the  high  and  low  field  forms  of  the  dis- 
charge and  the  effects  of  inert  and  detaching  additive  gases 
are  presented. 

Theory 

When  an  electric  field  is  applied  to  a plasma,  the 
charged  particles  are  accelerated  in  a direction  parallel 
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or  anti-parallel  to  the  field.  The  charged  particle  current 
density, j , can  be  determined  from: 


where  q is  the  electronic  charge,  N;  and  Ne  are  the  ion 
and  electron  number  densities;  and  V;  and  Ve  are  the  ion 
and  electron  drift  velocities.  Due  to  its  high  mobility  the 
electron  is  the  major  contributor  to  the  current  flow  in  a 
discharge.  The  electron  drift  velocity  can  be  calculated 
from: 


JE 


(6) 


This  relationship  combines  the  effect  of  the  accelerating 
force,  the  electric  field  (E) , and  the  effects  of  the  decel- 
erating force,  the  frequency  of  electron-neutral  collisions 
( Vew  ) • Combining  the  above  equations  and  disregarding  the 
ion  current  leads  to  the  following  relationship  for  the 
plasma  impedance : 
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The  discharge  current  versus  voltage  relationship  can  be 
determined  from  the  above  equation  by  integrating  the  cur- 
rent density  throughout  the  cross  section  of  the  discharge 
and  the  electric  field  along  the  length  of  the  discharge. 
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In  equation  7 the  electron  number  density  and  effec- 
tive collision  frequency  are  complex  functions  of  the  elec- 
tric field,  gas  temperature,  gas  pressure,  collisional  cross 
sections  of  the  constituent  species,  and  the  discharge  tube 
geometry-  Although  the  theoretical  determination  of  the 
impedance  is  complicated,  the  measurement  of  the  discharge 
impedance  is  normally  an  elementary  experiment.  However, 
in  the  case  of  oxygen,  the  measurement  of  the  impedance  in 
the  transition  regions  is  complicated  by  the  requirement  for 
a very  good  current  regulator.  Previous  low-field-form  in- 
vestigations without  good  current  regulation  are  subject  to 
the  suspicion  of  the  discharge  acting  as  a circuit  element 
in  an  oscillatory  circuit  rather  than  just  exhibiting  posi- 
tive column  oscillations. 


The  general  experimental  setup  described  in  Chapter  II 
was  utilized  during  the  determination  of  the  discharge  impe- 
dance characteristics.  X-Y  recorders  were  used  for  record- 
ing the  tube  potential  and  electric  field  as  functions  of 
the  discharge  current.  The  discharge  current  was  slowly 
swept  from  both  low  to  high  current  and  from  high  to  low, 
thus  allowing  the  hysteresis  effects  to  be  displayed.  The 
electric  field  was  measured  by  probes  inserted  through  ports 
3 and  4 . 

The  gases  were  flowed  through  the  tube  at  rates  calcu- 
lated to  maintain  a constant  residence  time  regardless  of 


pressure.  Temperature  effects  were  not  taken  into  account. 
The  flow  was  adjusted  according  to  the  following  relation: 

F=P*6.7S  (8) 

where  F is  the  gas  flow  in  standard  cc  per  min.  The  velo- 
city of  the  gas  in  the  tube  Vg  was: 


( 


(91 


where  Fa  is  the  actual  gas  flow  in  cc  per  sec  and  is  deter- 
mined by: 


Fa*4.32x  id2-  ,10) 

The  residue  time  of  the  gas  in  the  tube,'!",.^  , was: 

^res  ~ (11) 

where  L is  the  length  of  the  discharge.  For  T = 300°K, 

L = 50.5  cm,  R„=  .95  cm:  Vg  = 31.02  cm/sec,  “Y  res  = 1.63 
sec. 

Experimental  Results 

Pressure  Effects  and  Hysteresis.  Figures  9,  10,  and 
11  depict  the  tube  potential  and  electric  field  variations 
as  functions  of  the  tube  current.  These  figures  depict 
typical  curves  for  the  parameters  specified.  Experimental 
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repeatibility  was  good  between  runs,  especially  in  regard 
to  the  potential  and  electric  field  values.  The  greatest 
variation  occurred  in  the  current  values  at  which  the  dis- 
charge transitioned  between  forms.  Generally  this  variation 
was  less  than  10%  for  the  high  current  transition.  The  low 
current  transition  variation  was  greater.  However,  it  is 
postulated  that  this  variation  was  due  to  the  reduced  stabi- 
lity of  the  current  regulator  at  low  currents  and  the  high 
negative  impedance  of  the  discharge.  Often,  at  or  near  the 
point  of  the  low  current  transition,  the  regulator  would  be- 
gin to  oscillate;  and  the  discharge  would  extinguish.  It 
was  impossible  to  map  this  transition  at  pressures  in  excess 
of  6 torr. 

Figure  9 depicts  the  potential  variation  as  a function 
of  current  at  various  pressures.  The  high  current  transi- 
tion between  the  high  and  low  field  forms  was  a very  strong 
function  of  both  current  and  pressure;  the  low  current  tran- 
sition was  a strong  function  of  current  only.  The  transi- 
tion between  the  forms  was  very  rapid  at  any  one  axial  posi- 
tion in  the  tube.  The  two  forms  often  existed  in  the  tube 
at  the  same  time.  When  this  occurred  the  high  field  form 
existed  on  the  anode  side  of  the  positive  column  and  the 
low  field  form  on  the  cathode  side.  Discontinuities  in  the 
tube  often  provided  the  interface  between  the  two  forms  when 
they  coincided  in  the  column.  This  simultaneous  existence 
of  the  two  forms  is  illustrated  by  the  stair-stepped  nature 
( of  the  transitions  in  figure  9.  As  different  axial  sections 
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of  the  tube  transitioned,  the  potential  trace  stepped  to  a 
new  value.  Often  the  electric  field  probes  acted  as  the 
spatial  discontinuities  for  the  interfaces  between  the  forms. 

In  the  high  field  form  the  discharge  was  stable;  there 
were  no  perceivable  oscillations  or  random  behavior.  In  the 
low  field  form,  however,  there  were  large  amplitude  oscilla- 
tions and/or  noiselike  behavior.  This  was  observed  from 
oscilloscope  and  spectrum  analyzer  traces  of  the  tube  cur- 
rent. In  the  low  field  form  there  was  considerable  jumping 
among  different  potential  values  as  the  current  was  swept. 
This  is  shown  only  on  the  7 torr  trace  of  Figure  9 but  ap- 
peared on  all  the  raw  data.  The  difference  between  the  low 
field  potentials  was  small  and  was  related  to  the  frequency 
of  oscillation.  As  the  oscillation  frequency  jumped  between 
different  stable  and/or  astable  modes,  the  potential  would 
jump. 

Figure  10  depicts  the  hysteresis  nature  of  the  transi- 
tion. The  discharge  tended  to  stay  in  the  form  in  which  it 
existed.  The  high  current  transition  always  exhibited  a 
predictable  hysteresis  phenomenon.  The  low  current  transi- 
tion did  not;  the  limitations  of  the  experiment  equipment 
precluded  a predictable  mapping  of  this  transition. 

Figure  11  is  a typical  electric  field  variation  for 
the  discharge.  The  transitions  of  the  electric  field  tracked 
those  of  the  tube  potential.  Due  to  the  fact  the  electric 
field  measurements  encompassed  only  a small  section  of  the 
axial  length  of  the  tube,  the  stair-stepping  during 
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transition  was  much  less  evident.  Usually  only  two  steps 
were  included  in  the  spatial  span  between  the  two  electric 
field  probes. 

Existence  Regions.  Figure  12  is  a plot  of  the  exis- 
tence regions  for  the  high  and  low  field  forms  as  functions 
of  the  discharge  current  and  pressure.  The  transition  points 
wer.e  taken  from  the  electric  f ield-vs-current  traces.  Dur- 
ing a high  current-to-low  current  sweep,  transition  was  con- 
sidered to  have  taken  place  when  the  first  portion  of  the 
region  between  the  two  probes  transitioned.  For  the  low-to- 
high  sweep  the  final  transition  was  used  for  the  transition 
point.  Points  from  several  different  experimental  runs  were 
plotted,  and  a least  squares  power  curve  fit  was  made.  The 
resulting  boundary  curves  for  the  high  current  transitions 
showed  good  correlation  (r^  = .96  to  .98)  with  the  data;  the 
low  current  transitions  were  poorly  correlated.  The  experi- 
mental results  of  Pekarek  are  shown  for  comparison*  (Ref  15: 
750). 

Flow  Effects.  An  experiment  was  conducted  to  deter- 
mine the  effects  of  gas  flow  rate  on  the  characteristics  of 
the  discharge.  Flow  rate  was  varied  over  approximately  an 
order  of  magnitude.  It  was  observed  that  the  impedances  in 
the  low  and  high  field  forms  of  the  discharge  did  not  vary 
significantly  as  functions  of  the  flow  rate.  However,  the 


*A  direct  comparison  is  difficult  to  make.  Pekarek' s 
experiments  were  conducted  with  a nonflowing  gas  in  an  un- 
cooled discharge  tube. 
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HIGH  FIELD  FORM 


point  of  transition  between  the  forms  did  change.  Figure  13 
depicts  the  low-to-high  and  high-to-low  transition  points 
as  functions  of  the  flow  rate.  As  shown  on  the  figure  the 
transitions  occurred  at  lower  currents,  and  the  hysteresis 
became  more  pronounced  as  the  flow  was  increased. 

Effects  of  Additives.  Other  gases  were  mixed  with  -the 
oxygen  to  determine  the  effect  on  the  discharge  characteris- 
tics. Carbon  monoxide  and  argon  were  chosen  as  the  additive 
gases  - carbon  monoxide  because  it  is  a strong  detacher  of 
the  electron  from  the  0 ion  (Ref.  31:597-599)  and  argon  be- 
cause it  is  inert  and  should  have  only  minor  effects  on  the 
characteristics.  During  the  experimentation  the  oxygen  flow 
was  set  at  a high  rate*,  and  the  mix  percentage  was  con- 
trolled by  adjustment  of  the  additive  gas  flow  rate.  The 
partial  pressure  of  the  oxygen  was  maintained  constant  while 
the  total  pressure  was  allowed  to  vary  with  the  additive. 

The  mixture  percentages  specified  in  figures  14  through  17 
refer  to  volume  and  were  calculated  by  dividing  the  additive 
flow  rate  by  the  oxygen  flow  rate. 

Figure  14  depicts  the  effect  of  the  addition  of  argon 
on  the  impedance  characteristics  of  the  discharge.  As  shown 
in  the  figure  the  effect  was  minor.  The  impedances  of  the 
high  and  low  field  forms  remained  basically  unchanged.  The 
transition  point  between  the  two  forms  was  slightly  effected. 


*400-500  seem,  depending  on  the  vacuum  system's  abi- 
lity to  maintain  the  desired  pressure. 
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Figure  15.  Argon  Additive  Effects  on  Transition 
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as  depicted  in  figure  15. 
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The  effect  of  carbon  monoxide  on  the  discharge  was 
much  more  pronounced,  as  depicted  in  figure  16.  As  carbon 
monoxide  was  added,  the  DC  impedance  of  the  low  field  form 

i 

increased  to  the  point  where  the  difference  between  the  two 
forms  essentially  ceased  to  exist.  The  transition  between 
the  two  forms  generally  decreased  to  lower  currents,  and 
the  hysteresis  became  less  pronounced  as  the  additive  was 
increased.  Figure  17  depicts  the  carbon  monoxide  effect  on 
the  transition  point.  For  both  the  4 and  8 torr  experimen- 
tal runj  there  was  no  mixture  that  completely  eliminated  the 
transition;  however,  as  illustrated  by  the  18.6%  run  in 
figure  16,  the  difference  in  impedance  across  the  transition 
was  minor. 


IV.  Discharge  Heating  and  Gas  Temperature 


The  number  density  of  the  neutral  particles  N is  one 
of  the  key  parameters  of  a partially  ionized  gas.  This  para- 
meter, used  in  conjunction  with  the  electric  field,  was  re- 
quired for  comparison  of  the  experimental  results  of  this 
study  with  the  modeling  results  and  with  other  experiments. 
The  number  density  was  calculated  by  using  the  perfect  gas 
law* : 


N = 


The  determination  of  the  pressure,  P,  was  easy;  the  determin- 
ation of  the  temperature,  T,  was  very  difficult. 

This  chapter  addresses  the  determination  of  the  gas 
temperature.  In  the  first  portion  of  the  chapter  a relation- 
ship for  the  radial  and  axial  gas  temperature  is  determined 
using  the  assumption  that  all  the  input  energy  is  deposited 
in  gas  heating.  Next  the  setup  used  to  experimentally  deter- 
mine the  gas  temperature  is  presented.  This  is  followed  by 
the  experimental  results.  Finally,  these  results  are  used 
to  develop  a relationship  between  gas  temperature,  input 
electrical  power,  and  gas  pressure.  Once  determined,  this 
relationship  was  used  in  reducing  all  of  the  experimental 
data.  As  an  example,  the  relationship  between  E/N  and 


*Valid  for  the  pressures  used  in  this  study. 
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current  is  presented  for  various  pressures. 

Theoretical  Determination  of  Gas  Temperature 

The  following  analysis  is  based  on  the  assumptions 
that  all  of  the  electrical  energy  input  ends  up  in  gas  heat- 
ing and  that  the  heating  is  uniform  in  the  axial  and  radial 
directions.  The  first  assumption  is  obviously  not  true;, 
some  of  the  energy  goes  into  excitations  which  decay  into 
optical  channels*  or  into  long-lived  metastable  levels  which 
flow  out  of  the  tube  before  decaying.  The  second  assumption 
for  the  radial  direction  is  good  for  the  low  field  form 
where  all  of  the  experimental  measurements  were  made.  Thomp- 
son has  shown  that  the  energy  pumping  species,  the  electrons, 
have  a generally  constant  density  profile  when  the  negative 
ion  to  electron  concentration  ratio  is  high,  as  is  true  of 
the  low  field  form  of  the  discharge  (Ref  20:514).  Due  to 
oscillations  the  axial  heating  is  definitely  not  uniform  at 
any  one  point  of  time.  Time  averaging,  however,  should  al- 
low this  assumption  to  be  reasonably  good. 

The  theoretical  determination  of  the  gas  temperature 
is  based  on  balancing  the  rates  of  energy  gain  and  losses. 
Consider  the  incremental  volume  between  the  two  concentric 
cylinders  of  figure  18. 

*This  is  a very  small  amount  as  evidenced  by  the 
limited  optical  output  from  the  discharge. 


Figure  18.  Geometry  for  Temperature  Determination 
The  energy  rate  balance  equation  is*: 


LLj  = ^Q.co"  ^ci]  ■*" 


Rate  of  energy  Net  radial  Net  flow 

gain  in  the  volume  conduction  loss  loss 

The  energy  gain  from  the  electric  field  is: 


(13) 


I*.  _ $78-  E- 1 - r-c/r  , 

'^3  "*  (cal  sec 


(14) 


The  radial  conduction  of  energy  into  or  out  of  a cylindrical 
surface  with  height  dz  is: 


Uc=  K'\VT*dA  = 2'IT*  K 


.21 

r 


(15) 


Where  K is  the  thermal  conduction  coefficient  for  the  gas. 
Using  the  first  two  terms  of  a Taylor  series  expansion  to 
relate  the  temperature  at  r + dr  to  the  temperature  at  r, 
the  net  rate  of  radial  conduction  loss  is: 


flow. 


*Axial  conduction  was  assumed  to  be  small  compared  to 
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The  axial  flow  of  energy  into  or  out  of  the  concentric  cir- 
cular areas  is: 


UF  = N-k-T-'H'- r-c^r.  Vg  (17) 

1 

where  k is  Boltzmann's  constant.  The  net  rate  of  axial  flow  1 

loss  is:  | 

i 

Uf0-Urt-  f-N’K-H'-r-dr-Vg (18) 

C»  ; 

where  the  first  two  terms  of  a Taylor  series  expansion  were 
used  to  relate  the  temperature  at  z + dz  to  the  temperature 
at  z.  Equation  12  was  used  to  determine  the  value  of  N;  and 
equations  8,  9 and  10  were  used  to  determine  the  value  of  Vg. 

After  substituting  these  equations  into  equation  18  and  com- 
bining it  with  equations  14  and  16  into  equation  13: 


.530 *E  -I  = - 


3.88v  \ 6 


where  the  value  of  6.18  x 10~~*  cal/(sec)  (cm^)  (°c/cm)  was 
used  for  K (Ref  32:E-2).  In  solving  equation  19  for  T,  it 


was  assumed  that  the  radial  and  axial  functional  relation- 


T Cfr'jZ')  — T(iO  -Tcz) 


(20) 


T(r)  represents  the  radial  profile  at  any  position  z;  T(z) 
is  the  axial  temperature  at  r = o.  Substituting  equation 
20  into  equation  19: 

,530-E-I  = (21) 

The  solution  for  the  radial  profile  is  a parabolic  function: 


TCr)  = 


I + 


n\  [Ik 

Rt 


(22) 


where  Tw  is  the  wall  temperature  and  it  was  assumed  to  be 
constant-  The  solution  for  the  axial  dependence  is  an  ex- 
ponential : 


T0(z)=  Tv  + 3.08  X IO-E-I 


(23) 


where  it  was  assumed  that  the  gas  entered  the  tube  at  z = o 
with  temperature  TQ(z)  = Tw.  For  any  significant  axial  po- 
sition along  the  positive  column,  the  exponential  term  in 
equation  23  is  negligible  and  can  be  neglected.*  Therefore, 
the  centerline  temperature  can  be  calculated  by: 


♦Experiments  were  run  with  an  order  of  magnitude 
change  in  flow  rate.  No  appreciable  effects  on  the  center- 
line  temperature  were  noted. 
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(24) 


To(2)  = Tw  + 3.08X  I02.  E-I 

The  average  temperature  can  be  calculated  by  combining  equa- 
tions 22  and  24  and  integrating  over  the  radius  of  the  tube. 

<T>  = + 2.05xl0Z-  E-I  (25)' 


Experimental  Determination  of  Gas  Temperature 

Experimental  Setup.  The  general  experimental  setup 
described  in  Chapter  II  was  utilized  with  one  addition:  a 
small  thermocouple  was  inserted  into  port  number  two  in  such 
a position  that  it  could  be  used  to  measure  centerline  tem- 
perature. An  additional  identical  thermocouple  was  immersed 
in  ice  water  and  was  utilized  as  a reference.  The  potential 
difference  between  the  two  junctions  was  read  out  from  a 
digital  voltmeter  and  was  used  to  determine  the  gas  temper- 
ature . 

Experimental  Results.  Runs  were  made  at  pressures  be- 
tween two  and  nine  torr.  Only  readings  in  the  low  field 
form  were  possible.  When  the  discharge  transitioned  to  the 
high  field  form,  the  probe  would  begin  to  glow,  and  the  tem- 
perature reading  would  rapidly  increase  to  a value  above  the 
theoretical  maximum.  This  incandescence  was  possibly  caused 
by  the  recombination  of  atomic  oxygen  at  the  probe  surface. 

The  raw  data  at  any  one  pressure  was  insufficient  to 


allow  temperature  prediction  for  the  total  range  of  power 
inputs  experienced  at  that  pressure.  It  was  important  that 
all  of  the  data,  for  all  pressures,  be  utilized  to  allow 
analytic  interpolation  and  reasonable  extrapolation  to  cover 
the  full  parameter  space.  The  data  was  reduced  to  a rela- 
tionship of  centerline  temperature  rise*  as  functions  of 
power  input  per  molecule  and  system  pressure.  The  wall  tem- 
perature was  measured  throughout  the  experiment  and  averaged 
290°K  with  very  little  variation.  The  centerline  tempera- 
ture rise  was  calculated  from: 

A Tc  = T thermocouple  - 290°K  (26) 

The  power  input  per  molecule  was  calculated  from: 

p - JT'X  . _1 (27 ) 

<n>  tr*R£ 

where  equation  12  was  used  to  calculate  the  average  gas  den- 
sity, <nX 


<N> 


9.&£xiOle-P 

<T> 


(28) 


The  experimental  relationship  between  temperature  rise,  AT^  / 
to  reduced  power  input  is  depicted  in  figure  19. 

Analytic  Relationship  for  Gas  Temperature.  The  shapes 


*Centerline  temperature  rise  was  the  difference  be- 
tween the  measured  centerline  temperature  and  the  measured 
wall  temperature. 
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of  the  curves  in  figure  19  indicated  that  an  analytic  ex- 
pression for  centerline  temperature  rise  could  be  developed 
from  independent  functions  of  reduced  power  input  and  pres- 
sure. 

ATc.  = f (Pm)-  <3  (P)  (29) 

The  reduced  power  input  function.  Pm,  was  fit  to  a power 
curve  and  the  pressure  function,  P,  to  a quadratic.  The 
following  relation  was  determined  to  fit  all  of  the  data 
with  good  correlation  (r^  = .99): 

ATc  = [3.23  +-  2 1 . 1 * P - j.05-P2)'(PM]O:7^  (30) 

Equations  25,  27,  and  28  were  used  for  the  determination  of 
Pm.  Equation  30  is  transcendental  and  was  solved  by  itera- 
tive computations. 

Figure  20  depicts  the  results  calculated  for  9,  6 and 
3 torr  using  the  above  expression.  Also,  the  theoretical 
maximum  temperature  and  the  experimental  data  are  plotted. 

It  is  to  be  noted  that  the  calculated  curves  show  initially 
a negative  second  derivative,  due  possibly  to  an  increase  in 
the  amount  of  energy  going  into  nonheating  channels  as  the 
power  input  is  increased.  As  the  power  input  is  increased 
further,  however,  the  second  derivative  becomes  positive. 
This  is  possibly  caused  by  the  rapidly  decreasing  gas  den- 
sity as  the  temperature  rises. 
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Figure  20.  Analytic  Temperature  Rise 
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Reduced  Electric  Field  Variation 


The  analytical  expression  for  the  temperature  rise  in 
the  discharge,  equation  30,  was  used  in  equation  28  to  cal- 
culate the  gas  density,  (N) , as  a function  of  input  power 
and  pressure.  This  calculation  was  used  to  determine  the 
reduced  electric  field,  (E/N) , as  functions  of  pressure  and 
current.  Figure  21  depicts  these  relationships  for  three 
pressures:  8,  6 and  4 torr.  The  current  sweep  was  from 

high-to-low. 
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V.  Low  Field  Oscillations 
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The  unique  characteristics  of  its  low  field  form  sets 
the  oxygen  plasma  apart  from  the  plasmas  of  other  gases. 

This  chapter  reports  the  results  of  experiments  that  inves- 
tigated the  oscillations  that  existed  in  the  low  field  form. 
The  results  of  these  measurements  were  not  used  in  valida- 
ting the  models,  since  the  models  were  limited  to  only  the 
high  field  form  of  the  discharge.  These  experiments  were 
run  to  provide  a better  insight  into  the  oxygen  plasma. 

Two  aspects  of  the  low  field  oscillations  were  inves- 
tigated. First  the  dispersion*  was  determined  for  various 
pressures  and  discharge  electrical  currents.  Secondly,  the 
nonsteady  electric  field  was  measured  as  a function  of  the 
same  parameters. 

Dispersion  of  the  Oscillations 

Experimental  Setup  and  Procedure.  In  order  to  deter- 
mine the  characteristics  of  the  low  field  oscillations , the 
experimental  setup  depicted  in  figure  22  was  used  in  con- 
junction with  the  general  experimental  setup  described  in 
Chapter  II.  A function  generator  operating  in  a pulse  mode 
was  utilized  to  pulse  the  current  regulator  at  the  frequency 
of  interest.  The  low  field  oscillations  were  locked  into  a 


*Functional  relationship  between  the  wave  number,  k, 
and  the  oscillation  f requency , u>  . 
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single  frequency  stable  mode  by  this  method.  The  output  of 
the  function  generator  was  also  connected  to  a frequency 
counter  for  a backup  determination  of  the  oscillation  fre- 
quencies. 

A photomultiplier  tube  was  mounted  inside  a casing 
with  a narrow  input  slit.  The  optical  output  from  the  dis- 
charge was  focused  on  the  slit.  A Gaertner  micrometer  tra- 
verse was  used  to  move  and  accurately  measure  the  position 
of  the  photomultiplier  along  the  length  of  the  discharge. 

The  output  of  the  photomultiplier  was  used  in  conjunction 
with  the  variable  time  delay  feature  of  the  oscilloscope  to 
accurately  determine  the  time  between  oscillation  pulses. 

This  was  the  primary  method  of  determining  the  oscillation 
frequencies.  The  wavelengths  of  the  oscillations  were  deter- 
mined by  synchronizing  the  oscilloscope  to  the  function 
generator  and  using  the  traverse  to  determine  the  distance 
between  pulse  peaks. 

Experiment  Results.  The  output  from  the  photomulti- 
plier showed  that  the  oscillations  consisted  of  very  short 
pulses  of  light  with  a relatively  long  dark  period  between 
pulses.  This  large  intensity  ratio  indicated  a high  level 
of  modulation.  Figure  23  depicts  a typical  photomultiplier 
output  signal. 

Figures  24  and  25  are  the  dispersion  characteristics 
for  pressures  of  4 and  8 torr  respectively.  The  observed 
oscillation  frequencies  ranged  from  20KHz  to  90KHz.  This 
falls  within  the  2 to  lOOKHz  range  reported  by  other 


( k-raq/  sec.) 


(K- RAD. /SEC) 


investigators  (Ref  15:751,  Ref  18:11,  Ref  29:85,  Ref  16:305- 
317).  For  comparison  the  results  of  Pekarek’s  experiments 
at  2 torr  in  a 0.9  cm  radius  tube  (Ref  15:749-753)  and  Saba- 
dil’s  experiments  at  .65  to  3.2  torr  in  a 2.5  cm  radius  tube 
(Ref  16:305-317)  are  shown  in  figure  24.  As  indicated  by 
the  positive-value  intercepts  of  the  dispersion  curves  with 
the  to  axis,  the  oscillations  observed  during  this  study 
usually  exhibited  a small  amount  of  normal  dispersion.*  The 

C 

phase  velocities  varied  between  2.6  x 10  cm/sec  and  3.5  x 

C C 

103  cm/sec  and  the  group  velocities  between  2.5  x 10  cm/sec 
and  3.1  x 10^  cm/sec.  Sabadil's  results  indicated  the  os- 
cillations to  be  dispersionless  with  a velocity  of  2.4  x 10^ 
cm/sec  for  observations  between  .65  and  3.2  torr  (Ref  16: 
305-317).  Pekarek's  observed  velocities  varied  between  2.7 
and  3.6  x 10^  cm/sec.  Since  the  smount  of  dispersion  deter- 
mined by  this  study  was  so  small  and  exhibited  some  scatter, 
the  experimental  evidence  is  considered  to  be  indicative  but 
not  conclusive  in  contradicting  Sabadil's  nondispersion  ob- 
servations. 

As  stated  in  Chapter  III,  it  was  observed  that  the 
voltage  across  the  tube  in  the  low  field  form  would  jump 
back  and  forth  between  several  different  values  and  that 
these  values  were  related  to  the  oscillation  frequency. 
Figures  26  and  27  depict  the  variation  in  tube  voltage  with 
changes  in  the  oscillation  frequency.  It  is  to  be  noted 


*The  phase  velocity  decreased  with  increasing  fre- 
quency. 
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Figure  26.  Potential  Variations  at  4 Torr 


POTENTIAL  (K-V.) 


that  it  was  possible  to  lock  the  oscillations  only  into  cer- 
tain stable  modes,  as  indicated  by  the  data  points.  Frequen- 
cies of  oscillation  between  these  stable  modes  were  impossi- 
ble to  achieve.  Therefore,  the  lines  are  drawn  on  the  fig- 
ures only  to  indicate  the  trend  of  the  data. 

Electric  Field  Shape 

Experiment  Theory.  The  objective  of  the  experiment 
was  to  measure  the  electric  field  variations  in  the  low 
field  form  of  the  discharge.  To  do  so  a method  was  required 
to  transform  the  observable  measurements  of  potential  differ- 
ence into  a determination  of  the  electric  field.  At  any 
axial  point,  z,  in  the  discharge  the  potential  referenced 
to  one  of  the  electrodes  can  be  represented  by  a Fourier 
series  expansion: 

oo  . 

V (z,t)  = -§x5+'X  (a-teVcos^-wt-kN-z}  + (31) 

buCzV  SIN  -KwZ 

The  electric  field,  E,  evaluated  at  an  arbitrary  point  z = o 

is : 


e Cz,t)  = 


3V(z»t) 

oZ 


■Z-O 


(32) 
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oo 
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N=  I 

(^-  + kN-aN)-SlN,{KI'CO-i}]  ( 


Assume  that: 


dz  ^ KM'dM 


^-1  < k-bN 


This  assumption  is  valid  if  the  oscillation  electric  field 
pulse  does  not  change  with  position.  In  the  region  of  mea- 
surement in  the  discharge  the  optical  pulse  shape  did  not 
change  significantly  as  a function  of  position,  which  is  con- 
sistent with  the  assumption.  The  electric  field  evaluated 
at  an  arbitrary  point  z = o is  therefore: 

OO 

E(zvO= -j- gf5  4-y(-KN'bN<0S{w-£J-t]  + (36) 

\ 

Kn'&w  ] 

One  method  of  determining  the  electric  field  is  to  use 
two  axially  separated  probes  located  in  the  plasma  and  mea- 
suring the  potential  difference.  This  method  has  serious 


disadvantages  if  the  probes  are  separated  by  distances 
greater  than  the  wavelength  of  any  of  the  significant  har- 
monics. For  the  range  of  wavelengths  observed  in  this  ex- 
periment, the  distance  between  the  probes  would  have  had  to 
have  been  less  than  .3mm  to  faithfully  reproduce  just  the 
fundamental  frequency. 

The  method  that  was  utilized  involved  a time  deriva- 
tive of  the  signal  from  a single  probe  referenced  to  an  elec 
trode  of  the  discharge.  Taking  a time  derivative  of  equa- 
tion 31  at  an  arbitrary  point,  z = o,  yields: 


CJ> 

= hl'OJ  *3.n*  SlN  ^KJ  *60 -t  j + N'Cd-bK*  } 


If  the  medium  is  dispersionless,  or  very  nearly  so: 


Kkj  - (38) 

1XK1  V<1> 

where  V<j)  is  the  phase  velocity  of  the  oscillations.  There- 


fore : 


o© 

= - V<j>^T  ('  Kn‘  bw' C05  (xJ * t j 4 (39) 

U=l  ' 


The  optical  dispersion  experiments  reported  previously  in 
this  chapter  showed  that  the  oscillations  were  just  slightly 


dispersive.*  The  maximum  variation  in  phase  velocity  at  any 
one  pressure  and  current  was  less  than  23%  over  a frequency 
variation  of  4.1  to  1. 

In  measuring  the  potential  in  the  discharge,  the  char- 
acteristics of  the  interaction  of  the  probe  with  the  plasma 
had  to  be  considered.  Any  object  immersed  in  a plasma  will 
develop  a sheath  around  it.  This  sheath  can  be  considered 
to  present  a parallel  resistance  and  capacitance  coupling 
between  the  probe  and  the  plasma.  A rigorous  analysis  of 
the  coupling  for  a cylindrical  probe  is  very  complex  and  is 
beyond  the  scope  of  this  study.  However,  a simplified  ana- 
lysis can  be  used  to  determine  if  this  coupling  is  important 
in  determining  the  overall  measurements.  Consider  the  elec- 
trical circuit  of  figure  28. 
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Figure  28.  Equivalent  Circuit  of  Plasma  Sheath 


Rs  and  Cs  are  the  sheath  resistance  and  capacitance 


*Refer  to  figures  24  and  25. 
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respectively.  For  a probe  operating  at  the  floating  poten- 
tial, where  no  current  is  being  drawn  from  the  plasma,  the 
plasma  dynamic  resistance  can  be  approximated  by  (Ref  9:142- 
143)  : 


dVp  _ k-Te  _ 4- 

d I i=o  A*<]'  Ne*<V>f 


(40)  • 


where  A = area  of  the  probe  = 1.6  x 10  2 cm^ 

<V>=  mean  velocity  of  the  ions 

= 4.8  x 104  cm/s ec  for  molecular 
oxygen  ions  at  300°K 

...  mlO  electrons  for  low  field  form 

Ne  = 10  -.3 

cm  3 

Te  = 2 x 104  °K  for  low  field  form 
Rs  = 6 x 106  ohms 

The  capacitance  can  be  approximated  by  (Ref  9:128-129,  Ref 
33:21)  : 


CS 


A*  £o 

Xs 


(41) 


where  X s = sheath  thickness  ***  Ai>  546  10~^  cm 
Cs  = 1.4  x 10"13  farad 

The  combination  of  Rs  and  Cs  had  a time  constant  equivalent 
to  approximately  2 x 105  hertz.  Even  at  this  frequency  and 
higher,  their  values  had  little  effect  on  the  overall  mea- 
surements due  to  the  high  impedance  of  the  probe . 

Experimental  Setup.  Figure  29  illustrates  the  setup 
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utilized  to  measure  the  time  dependent  electric  field  in  the 
low  field  form  of  the  discharge.  To  determine  the  steady 
(d.c.)  portion  of  the  field,  the  differential  voltage  be- 
tween two  probes  was  passed  through  a low  pass  filter  and 
displayed  on  a digital  voltmeter.  The  varying  (a.c.)  portion 
of  the  field  was  determined  by  differentiating  the  potential 
of  one  of  the  probes  referenced  to  ground.  A simple  resis- 
tor and  capacitor  differentiator  was  constructed  and  cali- 
brated for  the  purpose.  The  output  of  the  dif ferentiator 
was  connected  to  the  vertical  axis  of  a Tektronix  555  oscil- 
loscope; the  horizontal  trace  was  synchronized  to  the  oscil- 
lation frequency.  The  deflection  plates  of  the  oscilloscope 
were  connected  through  a Tektronix  coupler  to  a HP70043  x-y 
recorder  to  provide  a permanent  record. 

Experimental  Results.  Figures  30  and  31  depict  elec- 
tric field  shapes  at  the  four  different  port  locations  for 
two  different  discharge  currents.  The  d.c.  values  of  the 
field,  determined  from  the  potential  differences  between  two 
probes,  were  added  to  the  a.c.  shapes  determined  at  the  spe- 
cified probe.  Therefore,  the  a.c.  values  were  measured  at 
a single  port,  and  the  d.c.  values  were  averaged  over  the 
region  between  two  ports  (10cm) . This  resulted  in  an  inac- 
curacy which  is  evidenced  by  the  electric  field  trace  pass- 
ing below  the  zero  line  on  some  of  the  traces.  E/N  values 
corresponding  to  the  electric  field  values  are  listed  in 
parenthesis  on  °ach  of  the  traces.  The  reduced  electric 
field  was  determined  utilizing  equation  30. 
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VI . Electron  Density 


The  number  of  free  electrons  in  a plasma  determines 
many  of  the  characteristics  of  that  plasma.  For  example  - 
since  the  electrons  are  the  principal  electric  current  car- 
riers, their  density  is  a key  parameter  in  determining  the 
impedance  of  a discharge.*  Also,  the  relative  density  of 
electrons  to  negative  ions  is  one  factor  in  determining  the 
stability  of  a discharge.**  Therefore,  to  describe  a dis- 
charge it  is  necessary  to  have  a knowledge  of  the  electron 
density.  The  objective  of  the  experiments  described  in  this 
section  was  to  determine  the  number  density  of  electrons  in 
an  oxygen  discharge  by  utilizing  a microwave  cavity.  The 
results  of  these  experiments  were  compared  to  values  calcu- 
lated from  current  continuity  equations  and  were  used  then 
to  validate  the  computer  modeling. 

The  chapter  starts  with  a description  of  the  theory 
associated  with  the  utilization  of  a microwave  cavity  for 
electron  density  measurements.  Next,  the  experimental  setup 
and  procedures  are  described.  This  section  is  followed  by 
the  calculations  of  the  factors  used  to  determine  the  elec- 
tron density  values  from  the  changes  in  the  cavity  resonant 
frequency  and  quality  factor. 


*Refer  to  equation  7. 

**This  aspect  is  discussed  further  in  Chapter  X. 
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Finally,  the  experimental  results  are  presented  along  with 
electron  density  values  calculated  from  simple  current  con- 
tinuity relationships. 


Experiment  Theory 

A conventional  method  of  determining  plasma  electron 
densities  and  collision  frequencies  consists  of  measuring 
the  changes  of  resonant  frequency  and  quality  factor*  of  a 
microwave  cavity  caused  by  the  introduction  of  the  plasma 
(Ref  35:196).  The  basic  measurable  characteristics  of  a 
cavity  (resonant  frequency,  wo,  and  quality  factor,  Q)  are 
not  only  functions  of  the  physical  dimensions  but  also  of 
the  properties  of  the  medium  contained  within  it.  If  the 
changes  in  these  measurables  are  sufficiently  small  when  a 
plasma  is  introduced  into  the  cavity,  perturbation  theory 
can  be  used  and  the  following  relations  utilized**  (Ref  36: 
441)  : 


Ab) 

OJ 


1 

2 • €o*  bj© 


CT;  - E *JV 

Sv  €,'Ea-JV 


(42) 


*The  quality  factor  is  a measure  of  the  losses  of  a 
cavity.  It  can  be  defined  as  a ratio  of  the  energy  stored 
in  the  cavity  to  the  energy  loss  per  cycle  (Ref  34:146). 

**Slater  did  not  include  the  relative  permittivity  in 
the  normalization  integral.  Due  to  the  particular  geometry 
of  this  experiment,  it  must  be  added  to  account  for  the 
highly  dielectric  quartz  tubes  and  carbon  tetrachloride 
coolant. 
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(Th-  E *dV 


where : 


5v^'E’- 


C 0©-  resonant  frequency  of  the  cavity  without  the  plasma 
C*Jp-  resonant  frequency  of  the  cavity  with  the  plasma 
CJ  = U)  p — Cc )o 

Qo  = quality  of  the  cavity  without  the  plasma 
Qp  = quality  of  the  cavity  with  the  plasma 

= T5r  " 'So 

V = cavity  volume  including  the  plasma  interaction  region 
V /=  plasma  interaction  portion  of  the  cavity  volume 
= complex  plasma  admittance 
£r=  relative  permittivity  of  dielectrics  in  the  bore 
E = electric  field 

The  validity  of  the  perturbation  theory  is  based  on 
the  electron  density  being  low  enough  to  permit  interaction 
of  the  electromagnetic  field  with  the  plasma  electrons  (Kef 
35:196)  : 


Ne 


where  Me  and  q are  the  mass  and  charge  on  the  electron  re- 
spectively. As  the  electron  density,  Ne,  approaches  the 
value  on  the  right  hand  side  of  the  inequality  44  the  plasma 
begins  to  shield  its  interior  from  the  externally  applied 


electric  field;  and  the  theory  becomes  invalid.  Addition- 
ally, the  pressure  should  be  sufficiently  low  that  (Ref  35: 
196)  : 


VeN  < 


(45) 


A high  value  of  collision  frequency,  V en,  could  result  in 
the  lowering  of  the  cavity,  Q,  to  the  point  where  higher  or- 
der cavity  modes  may  be  excited. 

Based  on  the  assumption  of  a Lorentz  gas  model*,  the 
complex  admittance  can  be  expressed  as  (Ref  8:194): 


(rr+jg;  = 


N&*  ^ • VeN  _ t Ne  * ^ • 6l)  p 

me*CVe2N4  6Jp)  me'CVeVGJ^) 


(46) 


Substituting  equation  46  into  42  and  43: 


CQp 

U)o  ~ 2'eo'CJo'rvie'(VeH+0Jp)  * 


Me*  E • dV 


s 

. 


v'Me-E'dV/ 


6p*  E * dV 


(47) 


(48) 


*A  Lorentz  gas  is  a plasma  in  which  the  electron  and 
ion  number  densities  are  less  than  the  neutrals;  only  the 
electrons  participate  in  the  phenomena  of  interest;  electron- 
neutral  collisions  are  the  only  collisions  of  concern;  and 
the  neutrals  are  assumed  to  have  infinite  mass  (Ref  8:191). 
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Solving  these  equations  for  electron  density  and  collision 


O 


frequency : 


k , r2  „ / _ 2-  me 

Ne  E-*  dv  “ 


/ A(q)  ' CJp  ^ _ 

1 + v 2^r) 


2. 

\j  = A f.  i ) r Qe— 

Veu  &\q)  2-acJ 


The  difficulty  in  solving  equation  49  lies  in  evaluating  the 
integrals.  Some  knowledge  of  the  spatial  variations  of  the 
electric  field  and  of  the  electron  number  densities  is  re- 
quired. These  integrals  are  evaluated  in  a following  sec- 
tion of  this  chapter. 


Experimental  Setup 

Figure  32  depicts  the  setup  utilized  to  determine  the 
change  in  cavity  resonance  and  the  quality  factor  as  func- 
tions of  discharge  pressure  and  current  (Ref  37:2-3,  Ref  38: 
1-6) . The  output  from  a microwave  sweep  frequency  oscilla- 
tor was  fed  through  a splitter  and  directional  coupler  to  a 
toroidal  cavity.  The  latter  was  mounted  concentric  to  the 
discharge  tube  between  probe  ports  numbers  three  and  four. 

A portion  of  the  microwave  oscillator  output  was  taken  from 
the  splitter  and  fed  into  an  absorption  wavemeter.  A detec- 
tor was  then  utilized  to  measure  the  strength  of  the  signal; 
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OSCILLOSCOPE  1 OSCILLOSCOPE  2 


its  output  was  connected  to  one  input  of  a differential  am- 
plifier. The  reflected  energy  from  the  cavity  was  passed 
through  the  directional  coupler  to  another  detector;  its  out 
put  was  connected  to  the  other  channel  of  the  differential 
amplifier.  The  output  of  the  differential  amplifier  con- 
trolled the  vertical  axis  of  one  beam  (A)  of  a dual  beam  os- 
cilloscope (1) ; the  horizontal  axis  was  driven  by  the  sweep 
signal  output  from  the  microwave  sweep  oscillator.  There- 
fore, the  beam  "A”  trace  depicted  the  cavity  response  versus 
frequency  of  the  cavity  resonance.  The  narrow  resonance  of 
the  tunable  wavemeter  provided  a mark  on  the  trace  which 
permitted  the  accurate  determination  of  the  resor.int  fre- 
quency and  Q of  the  cavity. 

The  setup  described  in  the  previous  paragraph  was  suf- 
ficient for  determining  the  characteristics  of  the  non- 
oscillatory  high  field  form  of  the  discharge.  However,  the 
varying  electron  densities  of  the  low  field  form  caused  the 
"A"  trace  to  jitter  excessively,  thus  negating  the  capabi- 
lity of  measuring  anything  but  the  average  values.  To  mea- 
sure the  characteristics  of  the  discharge  at  any  desired 
phase  point  of  an  oscillation,  the  experimental  setup  was 
modified  to  highlight  the  portion  of  the  "A"  trace  that  cor- 
responds to  that  phase  point.  A small  diameter  quartz  rod, 
acting  as  a light  pipe,  was  passed  through  the  cavity  to  the 
discharge  tube,  thus  connecting  the  optical  signal  of  the 
plasma  at  the  center  of  the  cavity  to  a photomultiplier  tube 
The  output  of  the  tube  was  used  to  control  the  vertical  axis 


of  oscilloscope  2 and  the  "B"  trace  of  oscilloscope  1.  The 
horizontal  axis  of  oscilloscope  2 was  driven  by  its  internal 
sweep  which  was  synchronized  through  an  adjustable  delay  to 
a capacitive  pick-up  looped  around  the  discharge  tube.  The 
horizontal  axis  of  the  "B"  trace  of  oscilloscope  1 was  dri- 
ven by  one  of  the  sweep  oscillators  in  the  scope;  this  os- 
cillator was  synchronized  to  the  gate  output  of  oscilloscope 
2,  thus  assuring  that  the  output  of  the  photomultiplier  was 
centered  on  the  trace.  The  other  sweep  of  oscilloscope  1 
was  operated  in  a delayed  mode  and  was  used  to  intensify 
both  the  "A"  and  "B"  traces.  Therefore,  using  the  delay 
control  of  oscilloscope  1,  a particular  phase  point  could  be 
selected  on  the  "B"  trace;  and  the  corresponding  resonance 
curve  on  the  "A"  trace  would  be  intensified.  For  all  runs 
in  the  low  field  form,  the  electron  density  measurements 
were  made  at  the  peak  and  minimum  points  of  the  optical  out- 
put. Figure  33  depicts  a typical  oscilloscope  1 trace  ac- 
quired during  a .low  field  measurement. 

Determination  of  Cavity/Plasma  Coupling 

The  early  microwave  cavities  utilized  for  measuring 
plasma  parameters  were  of  cylindrical  design.  This  type  of 
resonator,  although  easy  to  analyze,  was  restricted  to  rela- 
tively long  thin  discharge  tubes  for  which  the  tube  diameter 
was  smaller  than  either  the  length  or  diameter  o^  the  cavity. 
The  results  of  such  measurements  represented  the  values 
averaged  over  an  axial  length  of  the  discharge  which  was 
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several  tines  the  diameter  of  the  discharge  tube;  thus,  the 
cylindrical  cavity  was  incapable  of  providing  the  spatial 
resolution  necessary  to  investigate  oscillatory  plasmas. 

Toroidal  resonators,  however,  do  not  exhibit  such  a 
restriction*.  Due  to  the  small  region  of  interaction  be- 
tween the  cavity  and  the  plasma,  the  length  of  the  measured 
plasma  region  can  be  less  than  the  diameter  of  the  discharge 
tube.  This  factor  made  the  toroidal  cavity  the  device  of 
choice  for  measuring  the  parameters  of  an  oscillatory  column. 
The  disadvantage  associated  with  the  toroidal  cavity  was  the 
difficulty  of  analysis  (Ref  39:3979,  Ref  40:48-49). 

The  analysis  of  the  toroidal  cavity  leads  to  the  eval- 
uation of  the  integrals: 


• b * JV 


(51) 


E*  JV 


(52) 


The  analysis  of  integral  51  was  not  trivial  but  was  straight- 
forward. Its  value  was  determined  by  the  geometry  of  the 
cavity  and  discharge  tube;  the  accuracy  of  t a calculation 
depended  on  the  method  of  evaluation  and  the  assumptions 


*The  toroidal  resonator  used  in  this  study  was  not  the 
conventional  toroidal  resonator  but  a special  design  origi- 
nated by  M.V.  Sicha  to  give  good  spatial  resolution  (Ref  41: 
256).  Figure  35  illustrates  the  cross  section  of  the  reso- 
nator. 
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that  went  into  that  method.  The  analysis  of  integral  52 
would  have  been  impossible  if  assumptions  had  not  been  made 
concerning  the  spatial  distribution  of  the  electron  density. 
In  order  to  evaluate  this  integral  and  treat  it  as  a linear 
function  of  the  electron  density,  the  electron  spatial  den- 
sity was  assumed  to  vary  between  two  different  profiles. 

For  the  two-component-plasma  the  electron  density  was  assumed 
to  have  a Bessel  profile  in  the  radial  direction  and  a uni- 
form profile  in  the  axial  direction  (Ref  11:241): 

Ne  = <53) 

For  the  three-component-plasma  the  radial  and  axial  profiles 
were  both  assumed  to  be  uniform.  The  radial  uniformity  was 
shown  to  be  good  for  oxygen  plasmas  in  which  the  negative 
ion  to  electron  densities  ratio  was  greater  than  ten  (Ref 
24:819-820).  The  axial  uniformity  assumption  was  reasonable 
for  the  minimum  of  the  oscillation;  however,  it  greatly  re- 
duced the  accuracy  of  the  measurement  of  the  electron  den- 
sity at  the  peak  of  the  oscillations.  For  example  - using 
the  electrical  field  trace  of  figures  30  and  31  and  ports 
numbers  three  and  four  and  the  measured  phase  velocity  of 
figure  25,  the  axial  extent  of  the  oscillation  minimum  was 
29mm;  and  the  peak  was  1.3mm*.  Therefore,  the  measurement 
at  the  oscillation  peak  must  be  considered  as  an  averaged 


*As  compared  to  a cavity  gap  of  1.6mm. 
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value;  all  other  values  should  be  actual. 
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The  calculation  of  the  integrals  51  and  52  was  based 
on  the  assumption  that  the  electric  field  pattern  set  up  by 
the  pole  tips  of  the  resonator  was  quasi-stationary . With 
this  assumption  the  resonator  was  treated  as  a capacitor 
formed  by  the  pole  tips  in  parallel  with  a coaxial  line  (Ref 
40:49).  This  method  of  analyzing  the  toroidal  cavity  and 
its  interaction  with  the  plasma  was  realistic  as  long  as  the 
gap  width  was  small  compared  to  the  discharge  tube  diameter. 
For  a gap  width  that  was  7%  of  the  tube  diameter,  Sicha 
showed  excellent  agreement  between  theory  and  experiment 
(Ref  39:3984-3989).  For  the  cavity  utilized  in  this  study 
the  gap  width  was  also  7%  of  the  tube  diameter. 

In  computing  the  normalizing  integral  51  the  method 
of  Sicha  was  employed  (Ref  40:50-56).  The  resonator  space 
was  divided  into  four  sections:  the  slot  between  the  pole 
tips  forming  the  capacitor,  the  stray  field  space  (caused 
by  the  capacitor)  inside  the  cavity,  the  coaxial  line  approx- 
imation of  the  cavity  itself,  and  the  bore  where  the  dis- 
charge tube  and  plasma  were  introduced. 


Integral  51  was  determined  by  summing  the  integrals 
for  each  of  these  spaces*: 


€ r E * civ  = \€r*ErdV  -4-  \ erc*E2.*dV+ 


'V 


/Vp 

between  pole  tips  stray  field 
\ €rc*  4-  \ €re'E4'dV 

A/c 


(54) 


cavity 


bore 


where  €ri  and  Ei  were  the  relative  permittivities  and  elec- 
tric fields  of  the  defined  regions.  The  evaluation  of  the 
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*The  general  solution  for  the  toroidal  cavity  also  in- 
cluded a term  associated  with  interaction  of  the  cavity  and 
stray  fields: 


2 • Es  • COS  • dV 


Because  of  the  axial  symmetry  of  the  cavity  utilized  in  this 
study,  this  term  was  very  small  and  was  neglected. 


W# 
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first  three  integrals  on  the  right  hand  side  of  equation  54 
was  done  analytically;  the  last  integral  was  evaluated  by 
numerical  methods.  This  was  required  due  to  the  complicated 
nature  of  the  dielectrics  in  the  interaction  region  (see 
figure  34) . The  electric  field  between  the  pole  tips  was 
assumed  to  be  uniform: 


E = — 

1 d 

where  U is  the  voltage  between  the  pole  tips. 

The  relative  energy  integral*  was  determined  to  be: 


(55) 


'Vp 


(56) 


where  the  integration  was  performed  over  the  volume  between 
the  pole  tips,  Vp.  For  the  specific  geometry  utilized  in 
the  experiment  and  with  U and  Grp  assumed  to  be  unity: 


j€rp*  El*  * dV  = 1 4-.  4-  volt* 


.cm 


(57) 


A/p 


The  stray  electric  field  in  the  cavity,  set  up  by  the 
pole  tips,  was  assumed  to  be  semicircular: 


E 


z 


U 

TT-X 


(58) 


where  x is  the  radius  from  the  center  of  the  gap.  The 


♦Relative  energy,  since  the  permittivity  of  free  space 
was  not  included. 
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relative  energy  integral  of  this  field  was  determined  to  be: 


= 2-u-  Si-  lw{2i^s^)|+2.(i?vRi-4)] 


where  the  integration  was  performed  over  the  cavity  volume, 
Vc.  With  the  assumption  of  Cy~c  equal  to  unity: 


€rt'E^dV=  10.13  volt2. cm 


The  electric  field  in  the  cavity  was  determined  by  ap- 
plying a coaxial  line  approximation: 


F„0vz)  - 


. u- sin{k-  (l-z)} 

2-  LN{§a]-SIW  {(<•(£-£)] 


relative  energy  integral  of  this  field  vas  determined 


to  be : 


^Ha-dV=  ^r-U' 


2.  Ch-dl-  C~R)-Sln{K-(h-dll 


where : 


^ _ 2*  ft*  (63) 

The  resonant  wavelength,  » for  the  cavity  used  in  the  ex- 
periment was  14.03  cm.  Using  this  value,  equation  62  was 
evaluated  to  be: 


95 


( : E • dV  = 1.(1  5 volt  .cm 

'c  3 


The  final  integral  of  equation  54  was  calculated  by 
utilizing  a relaxation  technique  to  numerically  solve  Lap- 
lace's equation*.  The  outputs  of  the  computer  program  were 
the  potential  and  electric  fields  in  the  interaction  volume. 
The  square  of  the  electric  field  multiplied  by  its  associated 
relative  permittivity  was  numerically  integrated  over  the 
total  bore  volume,  VB,  to  provide  the  solution  of  the  final 
integral  of  equation  54 : 


\ Ev-b'  E4  • c)V  - 25, volt2.< 

A/,» 


The  calculation  was  performed  also  for  a bore  volume  without 
the  dielectrics.  The  resulting  integral  was  evaluated  to 
be  8.482  volt  cm.  Therefore,  the  bore  region  with  dielec- 
trics contained  considerable  more  energy  than  the  one  with 
an  £r  = 1.  This  demonstrated  the  necessity  of  including  the 
dielectrics  in  computing  the  bore  region  contribution  to 
equation  54. 

The  solution  of  integral  65  did  not  include  the  con- 
tribution of  the  small  quartz  rod  used  as  a light  pipe  to 
the  photo-multiplier  tube.  This  contribution  was  considered 


gram. 


*See  Appendix  A for  a description  of  the  computer  pro- 
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to  be  negligible  since  this  rod  occupied  less  than  .13%  of 
the  cavity  volume. 


t 

Figure  35  depicts  the  log  of  the  square  of  the  elec- 
tric field  for  a bore  region  without  the  quartz  tubes  and 
the  carbon  tetrachloride.  Figure  36  depicts  the  same  with 
the  inclusion  of  the  dielectrics.  Figure  37  is  a cut  of  the 
mid-plane  of  figure  36.  As  depicted  in  figures  35  and  36, 
the  dielectrics  had  a very  pronounced  and  stepped  effect  on 
increasing  the  electric  field  in  the  bore  of  the  cavity. 

In  the  piasma  volume  of  the  discharge  tube,  the  highest  coup- 

\ 

ling  occurred  at  the  circumference  in  the  region  of  the  cav-  j 

ity  gap.  The  electric  field  did  not  fall  off  very  rapidly 
in  the  radial,  r,  direction  - indicating  a weak  dependence 
of  microwave  coupling  on  the  electron  density  radial  profile. 

The  field  did  fall  off  rapidly  in  the  axial,  z,  direction  - 
indicating  a good  spatial  resolution. 

The  final  tabulation  of  the  integrals  on  the  right 
hand  side  of  equation  54  resulted  in  a value  of  51.34  volt^ 
cm  for  the  normalizing  integral,  51. 

The  computation  of  the  integral  52  was  accomplished 
in  a manner  similar  to  the  computation  of  the  final  integral 
of  equation  54.  The  only  difference  was  that  the  relative 


energy  was  integrated  only  over  the  bore  region  inside  the 
inner  discharge  tube,  VC  The  integration  was  based  on  the 
following  relation: 


Ne'  £2'  dV  * Ne:\  FCrVE  (n2)*c*V 


(66) 


V 
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Figure  35.  Cavity  Electric  Field  without  Dielectrics 
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where  the  density  profile,  F(r),  was  assumed  to  be  either  a 
Bessel  (equation  53)  or  uniform  profile.  The  results  of  the 
computation  were  linear  solutions  c"  the  integral  for  the 
two  profiles: 

Neo' 

for  the  Bessel  profile. 

Nei\E2Cr,z)-d  V = 1.717‘N  eo  volt2 .cm  (68) 

/V' 

for  the  uniform  profile. 

Figure  38  depicts  the  axial  response  of  the  toroidal 
cavity  for  the  uniform  electron  density  profile;  the  norma- 
lized contribution  to  the  interaction  integral  is  plotted 
as  a function  of  position  along  the  axis  of  the  discharge 
tube . 

The  cavity  computations  were  experimentally  checked 
by  inserting  dielectric  rods  along  the  axis  of  the  cavity. 
Four  different  dielectric  rods  were  used,  two  of  quartz  and 
two  of  pyrex.  Table  I lists  the  results  of  the  experiments. 
In  each  case  the  relative  dielectric  constant  used  for  the 
calculated  resonant  frequency  shift  was  the  published  aver- 
age value  for  that  type  of  material  (Ref  32:E-48).  No  at- 
tempt was  made  to  independently  determine  the  dielectric 
constant  for  the  calibrating  rods.  The  results  indicated 
an  average  agreement  of  13.3%  between  the  measured  and 


) ^ ' E toVdV  = 0X1 1 * KJe0 

V'  volt2. cm 


(67) 
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calculated  resonant  frequency  shifts.  Also,  the  larger  the 
diameter  of  the  calibration  rod,  the  better  the  agreement 
was.  The  experimental  determination  of  cavity  performance 
was  merely  used  to  validate  the  model;  the  results  of  the 
experiments  were  not  used  to  modify  the  calculated  factors. 


Table  I.  Cavity  Performance 


Material 

Radius 

(cm) 

€r 

Af 

Measured 

Af 

Calculated 

% Difference 

Quartz 

.150 

3.93 

1.90 

1.56 

-17.9 

Pyrex 

.250 

4.74 

6.50 

5.64 

-13.2 

Quartz 

. 313 

3.93 

8.10 

7.01 

-13.5 

Pyrex 

.365 

4.74 

13.40 

12.34 

- 7.9 

Experimental  Results 

Figures  39  through  44  depict  the  measured  centerline 
electron  densities  as  functions  of  the  discharge  current  for 
pressures  between  2 and  7 torr.  For  the  high  field  form  of 
the  discharge,  the  Bessel  profile  coupling  factor  (equation 
67)  was  used  in  reducing  the  data.  For  the  low  field  form 
both  the  Bessel  and  the  uniform  profile  coupling  factors 
(equations  67  and  68)  were  used. 

Electron  densities  were  calculated  for  all  of  the  runs 
using  the  electron  drift  velocity,  tube  current,  electric 
field  and  pressure.  The  results  of  these  calculations  are 
depicted  as  dashed  lines  in  figures  39  through  44.  The 
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Figure  43.  Electron  Density  at  6 Torr 
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calculations  were  performed  using  the  following  equation: 


N 


eo 


cj.(  * VpfEW-ir-K-dr 

Jo  Neo 


(69) 


For  the  Bessel  profile: 


NeCr)  = K1 


eo 


Jo  TO 


(70) 


and  for  the  uniform  profile: 


NeW-  N 


eo 


(71) 


Using  the  Boltzmann  code  results  of  Chapter  VIII  (figures 
51  and  56) , the  drift  velocity  was  assumed  to  vary  as: 


V,  - A * (■§■) 


(72) 


where  the  coefficients  A and  B were  determined  by  regression 
analysis  over  the  applicable  E/N  ranges. 


E/N  is  assumed  to  be  a function  of  electric  field,  pressure, 
temperature  rise  and  radius  - in  accordance  with  the  follow- 
ing equation: 


£ 

N 


E 

P 


-12 

1.035X10  • 


290  + AT-  AT 


(73) 


where  the  number  density  was  calculated  using  the  ideal  gas 
law  and  a parabolic  approximation  for  temperature  as  a func- 
tion of  radius  (equation  22) . 

For  the  high  field  form  of  the  discharge  there  was 
good  agreement  between  the  experimental  and  calculated  elec- 
tron densities  at  low  pressures  and  low  discharge  currents. 

At  higher  currents  and  pressures  the  two  results  tended  to 
diverge.  This  was  probably  due  to  the  extrapolation  required 
to  determine  the  discharge  gas  temperature.  As  presented 
in  Chapter  IV,  the  gas  temperature  was  experimentally  mea- 
sured only  for  low  values  of  power  input;  high  power  temper- 
atures were  extrapolations.  Since  the  electron  density  was 
calculated  from  the  drift  velocity  which  was  E/N  dependent, 
the  accuracy  of  the  calculated  values  was  dependent  on  an 
accurate  knowledge  of  the  gas  temperature  to  determine  the 
neutral  species  number  density. 

For  the  low  field  form  of  the  discharge,  experimental 
and  calculated  values  are  depicted  for  both  Bessel  and  uni- 
form electron  density  profiles.  These  results  provide  upper 
and  lower  limits  for  the  electron  density.  Thompson  has 
demonstrated  that  for  discharges  in  which  the  negative  ion 
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densities  were  greater  than  the  electron  densities,  the  elec- 
tron density  profile  tended  to  be  uniform.  Since  the  low 
field  form  of  the  oxygen  discharge  has  been  shown  to  be  rich 
in  negative  ions,  the  actual  electron  densities  probably 
were  very  close  to  the  uniform  profile  results.  These  re- 
sults showed  very  good  correlation  between  experimental  and 
calculated  values. 

For  each  of  the  low  field  profiles,  two  experimentally 
determined  electron  density  curves  are  depicted  - one  for 
the  maximum  of  the  oscillation  and  one  for  the  minimum.  The 
minimum  of  the  oscillations  corresponded  to  the  higher  of 
the  electron  density  curves.  This  was  consistent  with  cur- 
rent continuity  considerations  since  the  drift-velocity, 
number-density  product  for  both  the  maximum  and  minimum 
points  must  be  equal.  Since  the  electric  field*  of  the  min- 
imum point  was  lower  than  that  of  the  maximum,  the  minimum 
point  electron  density  was  required  to  be  higher. 

The  electron  number  density  of  the  maximum  of  the  os- 
cillation in  the  low  field  form  was  underestimated.  This 
was  due  to  the  fact  that  the  electric  field  pulse  of  the  dis- 
charge oscillation  was  so  short  that  there  was  spatial  aver- 
aging due  to  the  cavity  gap  and  temporal  averaging  due  to  the 
measuring  electronics.  The  number  density  of  the  minimum  of 
the  oscillation  should  have  been  accurate  since  it  was  mea- 
sured over  a rather  flat  minimum  of  the  electric  field  os- 
cillation. 

*And  thus  the  drift  velocity. 
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VII . Mass  Spectra 


The  initial  objective  of  the  mass  spectral  experiments 
was  to  determine  the  relative  densities  of  the  neutral,  pos- 
itively ionized  and  negatively  ionized  species  as  functions 
of  pressure  and  current.  Of  particular  interest  was  the- 
negative  ion  spectra,  since  the  negative  ions  play  such  an 
important  role  in  the  low  field  form  of  the  oxygen  plasma. 

Due  to  equipment  problems,  only  the  neutral  spectra  were 
measured.  The  results  of  these  neutral  species’  measure- 
ments were  useful  in  the  validation  of  the  discharge  model- 
ing . 

In  the  first  section  of  this  chapter  the  experimental 
setup  is  described.  This  is  followed  by  a description  of 
the  calibration  of  the  mass  spectrometer.  Finally,  the  ex- 
perimental results  are  presented. 

Experimental  Setup 

With  the  exception  of  the  discharge  tube,  the  general 
experimental  setup  described  in  Chapter  II  was  utilized. 

Figure  45  depicts  the  discharge  tube  used  in  the  mass  spec- 
tra experiments.  The  length  and  cross  section  measurements 

■ 

of  the  tube  were  identical  to  those  of  the  tube  used  in  the 
previous  experiments.  However,  the  mass  spectra  tube  was 
constructed  without  electric-field  probe-ports  and  utilized 

* 

pyrex  glass  instead  of  quartz.  These  differences  were  dic- 
tated by  the  critical  construction  and  mounting  requirements 
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imposed  by  the  necessity  of  attaching  the  tube  to  the  mass 
spectrometer.  A sampling  port  was  mounted  in  the  side  of 
the  tube  at  a location  that  corresponded  to  half  the  dis- 
tance between  electric  field  ports  three  and  four  of  the 
previous  tube.  A stainless  steel  sampling  tube  was  inserted 
through  this  port;  this  tube  extended  to  the  edge  of  the 
discharge.  Through  this  tube  samples  of  the  discharge  plasma 
were  drawn  into  the  mass  spectrometer. 

Figure  46  illustrates  the  electrical  setup  for  the 
experiment.  The  sampling  tube  was  referenced  to  ground  po- 
tential, and  the  discharge  tube  anode  and  cathode  were  al- 
lowed to  "float."  In  order  to  utilize  grounded  current  and 
voltage  indicators,  a laboratory  designed  and  constructed 
isolator  was  used.  This  unit  employed  high  voltage  optical 
isolators  to  protect  the  indicators  from  damage. 

The  mass  spectrometer  utilized  was  an  Extranuclear 
EMBA  II.  This  system  used  a 22-centimeter  quadrupole  to 
provide  a M/AM  discrimination  better  than  1500.  Sensitivity 
of  the  system  was  enhanced  by  the  modulation  of  the  sample 
beam  by  a tuning-fork-chopper  prior  to  its  entrance  into  the 
ionizer  and  subsequently  using  a Princeton  Applied  Research 
HR-8  Lock-In  Amplifier  for  synchronous  detection  of  the  sig- 
nal from  the  electron  multiplier.  This  scheme  allowed  the 
effect  of  the  background  gases  in  the  spectrometer  to  be 
minimized.  The  output  of  the  mass  spectrometer  was  recorded 
on  a Hewlett  Packard  Model  7004  x-y  plotter. 
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System  Calibration 

In  order  to  be  able  to  measure  the  atomic  oxygen  con- 
centration of  the  plasma,  it  was  necessary  to  operate  the 
ionizer  at  a potential  which  would  not  cause  a significant 
dissociative  ionization  of  the  molecular  oxygen  but  would 
still  provide  a strong  ionization  of  atomic  oxygen.  To  de- 
termine this  ionizer  potential  the  mass  spectrometer  was 
tuned  to  atomic  oxygen,  and  its  output  was  recorded  as  the 
ionizer  potential  was  varied.  Runs  were  accomplished  with 
both  the  discharge  on  and  off.  The  results  of  these  runs 
are  depicted  in  figure  47.  An  ionizer  potential  of  21  volts 
was  chosen  for  the  remaining  runs. 

Setting  the  ionizer  at  21  volts  did  not  result  in  21 
ev  ionizing  electrons;  the  work  function  of  the  filament  had 
to  be  accounted  for.  In  order  to  determine  this  offset  the 
mass  spectrometer  was  tuned  to  molecular  oxygen;  and  with 
the  discharge  off,  the  system  output  was  recorded  as  a func- 
tion of  the  ionizer  potential.  Specifically  of  interest  was 
the  behavior  in  the  potential  interval  near  the  threshold 
of  ionization.  A linear  extrapolation  of  the  slope  of  the 
ionization  curve  intercepted  the  ionizer  potential  axis  at 
18.4  volts*.  However,  this  had  to  correspond  to  the  ioniza- 
tion potential  of  12.1  ev  for  molecular  oxygen.  Therefore, 
the  offset  potential  for  the  system  was  determined  to  be 
6 . 3 volts. 


rt 


*See  the  inset  of  figure  28. 
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Figure  48.  Ionizer  Performance  with  Molecular  Oxygen 
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Experimental  Results 


Figure  49  depicts  mass  spectral  scans  for  different 
discharge  tube  currents  at  a pressure  of  5 torr.  Of  inter- 
est was  the  sudden  increase  in  the  atomic  oxygen  output  as 
the  discharge  transitioned  from  the  low  to  high  field  forms 
between  40  and  50  ma.  Also,  there  was  a lack  of  measurable 
neutral  ozone  in  the  discharge.  The  signals  at  18(H20), 
28(N2),  and  44{C02)  represented  background  gases  of  atmos- 
pheric origin.  Although  there  were  undoubtedly  small  leaks 
in  the  discharge  tube  portion  of  the  experiment,  the  back- 
ground in  the  spectra  was  the  result  primarily  of  a very 
small  leak  in  the  back  chambers  of  the  mass  spectrometer. 
This  was  determined  by  observing  the  discharge  tube  gas  flow 
and  pressure  effects  on  the  background.  The  background 
peaks  were  a function  of  discharge  tube  pressure  and  not  of 
the  rate  of  gas  flow  into  the  discharge  tube. 

Figure  50  depicts  the  variations  of  molecular  and 
atomic  oxygen  with  discharge  current.  The  experimental  runs 
were  made  at  a pressure  of  5 torr.  It  was  noted  that  when 
tuned  to  molecular  oxygen,  the  mass  spectrometer  exhibited 
a long-term  drift  phenomenon;  the  signal  output  decreased 
as  a function  of  time.  This  drift  effect  was  not  noticed 
when  observing  other  species  and  was  therefore  assumed  to 
be  caused  by  an  oxidation  of  the  detector  electrodes.  The 
molecular  oxygen  curve  of  figure  50  was  corrected  for  this 
effect  by  assuming  a linear  drift  and  requiring  the  low  cur- 
rent start  and  end  points  to  coincide. 
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VIII.  Electron  Energy  Distribution 

In  order*  to  develop  a mathematical  model  of  the  oxy- 
gen plasma,  it  was  necessary  to  know  the  transport  proper- 
ties of  the  electrons  and  the  reaction  rates  between  the 
electrons  and  the  other  species.  The  reaction  rates  deter- 
mine the  generation  of  many  of  the  species  that  contribute 
to  the  properties  of  the  plasma,  for  example:  the  excited 
species,  the  ions,  and  the  electrons  themselves.  In  this 
study  the  Boltzmann  equation  was  used  for  the  numerical  cal- 
culation of  the  electron  energy  distribution  function,  and 
from  it  the  electron  transport  properties  and  reaction  rates 
were  calculated. 

The  chapter  starts  with  a brief  explanation  of  the 
theory  associated  with  the  use  of  the  Boltzmann  equation. 
Next,  selected  results  of  the  calculations  are  validated  by 
comparison  to  experiments  reported  in  the  literature.  Fin- 
ally, all  of  the  results  used  in  other  portions  of  this 
study  are  presented. 

Theory 

The  fundamental  kinetic  description  of  a plasma  is 
provided  by  expressions  for  the  velocity  distributions  of 
each  constituent  specie.  Because  of  the  relatively  high 
pressure,  low  fields  and  approximately  equal  masses,  the 
distribution  functions  for  the  ions,  atoms  and  molecules 
were  assumed  in  this  study  to  be  Maxwellian  with  a 
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characteristic  temperature  equal  to  that  of  the  gas*.  The 
electron  velocity  distribution  was  not  Maxwellian  and  was 
not  in  thermodynamic  equilibrium  with  the  other  species. 

The  electron  velocity  distribution  was  determined  by  a nuraer 
ical  solution  of  the  Boltzmann  transport  equation. 

Solutions  of  the  Boltzmann  equation  provided  descrip- 
tions of  how  the  electrons  were  distributed  in  phase  space** 
and  in  time.  The  inputs  to  the  Boltzmann  equation  were  the 
externally  applied  forces  and  the  characteristics  of  the 
background  gas  which  affected  the  collision-induced  changes 
to  the  distribution  function;  the  output  of  the  equation  re- 
vealed how  the  electrons  were  distributed  in  velocity  and 
position.  The  Boltzmann  equation  is  a very  general  mathema- 
tical relation;  however,  it  was  difficult  to  solve  without 
involking  simplifying  assumptions.  In  the  absence  of  a mag- 
netic field  the  Boltzmann  equation  for  the  electrons  can  be 
written  as  (Ref  42:34): 

S + ^ We'  fc(E. 

fe  = distribution  function  of  electrons  such  that 

nefdvxdvydvzdxdydz  denotes  the  average  number  of  elec- 
trons in  the  volume  between  x and  x + dx,  y and  y + dy 

*Assumed  to  be  300°K. 

**Phase  space  is  a six-dimensional  space  in  which  each 
particle  is  represented  by  both  position  and  velocity  coor- 
dinates. 
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z and  z + dz  and  with  velocity  components  between  vx 

and  vx  + dvx,  Vy  and  Vy  + dVy;  and  vz  and  vz  + dvz 

and  fdvxdvgdvzdxdydz  = 1. 

Xe  = x,y,z  = position  coordinates  in  six-dimensional 
phase  space. 

ve  = vx»vy»vz  = velocity  coordinates  in  six-dimensional 
phase  space. 

Vx  = Del  operator  with  respect  of  position. 

Vv  = Del  operator  with  respect  to  velocity. 

E = applied  electric  field. 

Since  the  plasma  studied  was  weakly  ionized,  only  the  colli- 
sions between  electrons  and  neutral  particles  were  consi- 
dered. The  collision  term  of  equation  74  can  then  be  writ- 


ten as  (Ref  8:122) 


Me) 

^ I Collisions 


fe* fN " fe-fN  )' 3 ' S M ' SIN  {*}-  (75) 


dx  • dy  • d3VN 

where  electron  velocity  distribution  functions 

before  and  after  a collision, 
u-  neutral  species  velocity  distribution  func- 
tions before  and  after  a collision. 

<0  = relative  speed  between  the  electron  and 
neutral  particle. 

Soo=  differential  scattering  cross  section. 


X=  scattering  angle. 
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V = azimuth  angle. 

d VN=  velocity  differential  volume  for  the  neutral 
particles. 

Equations  74  and  75,  when  combined,  form  the  general 
form  of  the  Boltzmann  equation.  No  closed  form  solution  to 
the  general  equation  was  known;  therefore,  an  approximate 
solution  was  obtained.  An  external  electric  field  causes 
the  electrons  to  acquire  a drift  velocity  opposite  to  the 
direction  of  the  field.  This  directed  velocity  component 
has  the  effect  of  causing  the  electron  velocity  distribution 
function  to  be  anisotropic.  For  low  values  of  electric 
field  strength  the  drift  velocity  is  small  compared  to  the 
random  velocity  of  the  electrons.  Therefore,  the  electric 
field  effects  were  treated  as  a perturbation,  and  the  elec- 
tron velocity  distribution  was  expressed  as  the  sum  of  the 
first  two  terms  of  a spherical  harmonic  expansion  (Ref  7: 
274)  : 


(76) 


where  fe^°^  is  the  isotropic  distribution  function  and  f e ^ 
is  the  directed  distribution  function. 

Equations  74,  75  and  76  were  combined  into  a computer 
program  titled  EED  for  solution  of  the  Boltzmann  equation*. 


*Program  EED  was  developed  by  William  Bailey  and  Wil- 
liam Long  for  use  by  the  personnel  of  the  Plasma  Physics 
Laboratory  of  the  Aerospace  Research  Laboratories. 
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The  program,  as  listed  in  Appendix  C,  was  used  in  this  study 
to  calculate  the  isotropic  and  directed  electron  velocity 
distribution  functions,  the  electron  transport  parameters, 
and  the  electron-other  species  energy  dependent  reaction 
rates.  The  required  inputs  to  the  program  were  the  reduced 
electric  field  (E/N)  and  the  various  electron  collision 
cross  sections  (Appendix  B) . Both  elastic  collision  cross 
sections,  in  which  the  internal  energies  and  identities  of 
the  target  particles  remained  unchanged,  and  inelastic  col- 
lision cross  sections,  in  which  the  state  and/or  identity 
of  the  target  particles  were  changed,  were  provided  to  the 
program.  For  the  inelastic  collisions  the  electron  energy 
after  the  collision  was  decreased  by  an  amount  equal  to  the 
threshold  energy  for  the  collision.  For  ionizing  collisions 
the  secondary  electrons  were  considered  to  be  created  with 
zero  energy. 

The  computer  program  EED  calculated  the  electron  energy 
distribution  components,  f e ^ and  fe^^,  and  utilized  them 
to  calculate  the  electron  transport  parameters  from  the  fol- 
lowing moments  (Ref  43;2): 


O 


00 

1 {& 

= V fe  (Ue^  • Ue-  due 


where  is  the  mean  electron  energy. 


Vec»  “ 10'  (UeV^e'clUe 


127 


Si 


r 


(77) 


(78) 


r.  




t 


i 


where  Ved  is  the  drift  velocity  of  the  electrons. 

o*>  , 

O / \2-  [ , , 2 ccft 

o 

where  De  is  the  diffusion  coefficient  for  the  electrons  and 
Qm  is  the  momentum  transfer  cross  section. 


I - .Ik  . JL 

U-CH0R  Ved  N 


where-  Ue-CHflR  is  the  characteristic  energy  of  the  electrons, 
The  program  used  the  calculated  electron  energy  distribution 
functions  td  calculate  the  electron-heavy  particle  reaction 
rates  for  each  of  the  cross  sections  input  to  the  program. 
These  calculations  were  performed  by  the  relation  (Ref  9:42) 

oo 

= 10-  (fe(U«)  • Q,2(Ue)-Ui-  JUe  (81) 


where  R12  is  the  reaction  rate  for  the  process: 


e+  Si 


e+  Sa 


Q12  is  the  cross  section  for  the  above  process  and  was 
assumed  to  be  isotropic. 

Comparison  with  Experiment 

The  Boltzmann  code,  EED,  was  run  for  input  reduced 
electric  fields,  E/N,  ranging  from  1 Td  to  100  Td.  To  vali- 
date the  low  energy  end  of  the  results,  the  calculated  drift 


128 


i 


( 


velocities  were  compared  with  the  results  of  drift  tube  mea- 
surements. Figure  51  depicts  the  results  of  this  compari- 
son. Experimental  data  from  Brown  and  Nelson  were  used  for 
the  comparison  (Ref  44:91;  Ref  45:4081).  Figure  51  indicates 
the  close  correlation  between  the  calculated  and  experimen- 
tal values. 

To  validate  the  high  energy  portion  of  the  results,  a 
comparison  was  made  between  the  calculated  molecular  ioniza- 
tion rates  and  the  experimental  rates  of  Corbin  (Ref  46:81). 
In  order  to  make  a direct  comparison  the  code  parameters  of 
E/N  and  ionization  rate  were  converted  to  the  parameters  re- 
ported by  Corbin,  fzo/E  °<f/P20  - Th  is  conversion  was 

made  using  the  following  relations: 

F2.0  _ N-R-Tzo  3.03  X 10  , 

T ~ ' — E — - — EAi <82) 


' I 


J 

N-R.T20 


|6  p. 

= 3.3  X\0  .JSi- 

\/ecl 


(83) 


where  P2q  ~ pressure  reduced  to  20°C  (torr) 

R = Universal  Gas  Constant 

cXi  = Townsends  first  ionization  coefficient  (cm“^) 

T20  = 29  3°K 

Ri  = ionization  rate  (cm^  sec"-*-) 

The  comparison  included  the  drift  velocities  calculated  by 
EED  at  the  high  energy  inputs.  Again,  as  i,n  the  case  of  the 
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low  energy  results,  the  comparison  indicated  a close  corre- 
lation between  the  EED  results  and  experiments,  as  depicted 
in  figure  52. 

The  electron  energy  distribution,  fe^°^,  at  an  E/N  of 

70  Td  was  compared  to  the  experimental  result  of  Rundle  for 

an  E/P  of  25.6  v/cm  torr*  (Ref  47:147).  Rundle's  experiment 

indicated  an  exponential  dropoff  in  the  distribution  of 
i f y 

-.0608E  ; an  analytic  fit  to  the  distribution  of  the  pre- 

sent study  resulted  in  a dropoff  of  - . 0745E-*-  • ^ . Both  ex- 
periment and  theory,  therefore,  indicated  the  distribution 
was  neither  Maxwellian  (E^)  nor  Druyvesteyn  (E^)  but  occu- 
pied a position  inbetween.  The  comparison  is  shown  in  fig- 
ure 53. 

Results 

Figures  54  and  55  depict  the  isotropic,  -ft°(ou)/Ue  • an<^ 
directed,  (ue)/Uc » energy  distribution  functions  for  various 
values  of  reduced  electric  field,  E/N**.  Since  the  Boltz- 
mann equation  was  solved  by  EED  for  the  steady  state  condi- 
tion, the  depicted  distribution  functions  apply  only  to  the 
high  field  form  of  the  discharge. 

Figure  56  is  a plot  of  the  electron  drift  velocity  as 
a function  of  the  reduced  electric  field.  Figure  57  presents 


*At  273°K  an  E/P  of  25.6  v/cm  torr  equals  an  E/N  of 
72.4  Td. 

**These  plots  are  a function  of  energy  and  not  of  velo- 
city; also,  a factor  of  the  square  root  of  the  energy,  , 
has  been  divided  out  of  the  distributions. 
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Figure  55.  Directed  Distribution 
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the  mean  and  characteristic  energies  as  functions  of  E/N. 

Figures  58  through  62  depict  examples  of  various  ex- 
citation, ionization,  and  attachment  rates*.  The  examples-' 
chosen  for  presentation  are  those  which  were  utilized  later 
in  the  calculations  of  the  plasma  chemistry**. 

Figure  63  is  an  estimate  of  the  dissociation  rate  of 
molecular  oxygen.  It  was  determined  by  summing  the  rates 
associated  with  the  pumping  of  the  6.0  and  8.4  electronic 
excitation  levels.  The  cross  sections  associated  with  these 
excitations  are  reasonably  well  known;  however,  the  physical 
transitions  they  correspond  to  are  still  in  doubt.  There  is 
general  agreement  that  the  8.4  ev  cross  section  leads  to  the 
B35_u_  electronic  level  of  molecular  oxygen;  however,  the 
Franck  Condon  principal  indicates  that  this  transition  leads 
to  a dissociative  electronic-vibrational  level  which  results 
in  a breakup  of  the  molecule.  Both  the  4.5  and  6.0  ev  cross 
sections  have  been  suggested  to  be  associated  with  the  exci- 
tation of  the  c'Zu.  electronic  level.  Of  these  two,  the 
6.0  ev  cross  section  was  postulated  to  lead  to  a dissocia- 
tive level  (Ref  48)  . Normally  the  c'SIu.  excitation  would 
be  represented  by  a single  cross  section.  However,  to  se- 
parate the  results  of  the  collisions  into  those  which  would 


♦Refer  to  Appendix  B for  a further  explanation  of  the 
cross  sections  that  lead  to  these  rates. 

**Refer  to  Chapter  IX  and  Appendix  E for  a description 
of  the  chemistry  calculations. 
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IONIZATION  RATE  (cm^/SEC) 


DISSOCIATIVE  ATTACHMENT  RATE  (cmj  / SEC) 


IX.  Plasma  Chemistry 


In  a molecular  gas  plasma  various  physical  and  chemi- 
cal processes  occur  and  generate  a large  number  of  other  spe- 
cies. Positive  and  negative  ions  are  created;  free  electrons 
are  generated;  molecular  bonds  are  broken  down,  creating  ato- 
mic species  and  other  molecules;  species  associate  to  form 
clusters;  and  species  are  raised  in  energy  to  various  exci- 
tation levels.  In  order  to  understand  the  phenomena  occur- 
ring in  the  oxygen  discharge,  it  is  important  to  be  able  to 
determine  the  concentrations  of  the  various  species  and  their 
excitation  levels  as  functions  of  the  applied  electric  field; 
the  gas  pressures,  temperature,  and  flow;  and  the  physical 
geometry.  This  chapter  presents  the  results  of  a chemistry 
code,  titled  Chemco,  which  was  used  to  calculate  the  species 
concentrations  and  selected  excitation  levels  of  the  oxygen 
plasma  under  electrical  discharge  conditions  corresponding 
to  the  high  field  experiments. 

Chemco  incorporated  reaction  rates  between  the  various 
species,  the  transport  parameters  of  the  various  species, 
and  the  effects  of  the  discharge  tube  geometry.  The  inputs 
to  Chemco  included  the  initial  oxygen  molecular  density  en- 
tering the  tube,  N;  the  reduced  electric  field,  E/N;  the  E/N 
dependent  reaction  rates;  and  the  gas  flow.  The  effects  of 
temperature  and  pressure  were  taken  into  account  by  specify- 
ing N.  The  non-E/N-dependent  reaction  rates  were  taken  from 
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the  open  literature.  The  E/N-dependent  rates  were  taken 
from  the  results  of  the  Boltzmann  Code,  EED , computations*. 
Since  EED  and  Chemco  were  not  coupled,  any  single  Chemco  run 
was  valid  for  only  the  chosen  value  of  E/N.  This  limited 
the  application  of  the  results  to  the  steady,  or  high  field 
form,  of  the  oxygen  discharge. 

Output  from  Chemco  were  various  charged,  neutral,  and 
excited  species  densities  and  the  effects  of  each  of  the  in- 
cluded processes  as  functions  of  E/N  and  N.  These  outputs 
were  used  to  provide  additional  insights  into  the  oxygen 
plasma  and  to  draw  conclusions  regarding  the  causes  of  the 
observed  phenomena. 

The  chapter  starts  out  with  a description  of  the  che- 
mistry code-Chemco.  This  is  followed  by  the  results  of  a 
large  number  of  calculations  using  the  code.  Some  conslu- 
sions  are  presented  with  the  results;  however,  the  major  con- 
clusions are  presented  in  Chapter  X along  with  the  conclu- 
sions from  other  chapters. 


Chemco  Description 


General  Approach.  The  basic  equations  of  Chemco  were 


the  differential  rate  equations  of  each  of  the  included  spe- 
cies. These  equations  were  simultaneously  integrated  as 
functions  of  time  rather  than  being  directly  solved  for 
their  steady  state  solutions.  This  method  was  chosen  be- 
cause the  coupled  nonlinear  algebraic  equations  that  resulted 


from  a steady  state  analysis  were  very  difficult  to  solve 
for  the  number  of  species  and  degree  of  coupling  between 
equations.  For  each  run  the  calculations  continued  for  a 
time  equivalent  to  the  residence  time  of  the  gas  in  the  dis- 
charge tube.  It  was  found,  however,  that  the  species  concen- 
trations for  the  conditions  of  interest  reached  equilibrium 
in  a very  short  period  of  time. 

Chemco  is  a one-dimensional  code;  the  only  dimension 
considered  is  time  or  the  equivalent  distance  along  the  axis 
of  the  tube.  The  radial  distributions  of  the  species  are 
not  specifically  calculated;  however,  assumed-radial-profile 
effects  are  taken  account  of  in  calculating  the  diffusion 
losses.  Throughout  the  Chemco  runs,  total  mass  and  charge 
are  calculated  and  checked  to  assure  conservation. 

Chemco,  as  used,  incorporated  a total  of  69  chemical 
and  physical  reactions  between  the  various  species.  These 
reactions  were  considered  to  be  the  dominant  ones  over  the 
pressure  and  electron  density  ranges  of  this  study.  Initial 
selection  was  made  from  a list  of  130  potential  reactions*. 
Nonutilized  reactions  were  rejected  because  they  either  in- 
volved a species  not  carried  in  the  code,  had  a low  reaction 
rate  in  an  absolute  sense,  had  a low  reaction  rate  relative 
to  another  channel  with  the  same  net  reaction,  or  were  deter- 
mined by  early  Chemco  runs  to  be  insignificant. 

In  addition  to  chemical  and  physical  reactions,  Chemco 

*See  Appendix  D. 
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includes  the  effects  of  ambipolar  radial  diffusion  of  all 
the  charged  species;  radial  diffusion  of  the  neutral  species 
with  the  exception  of  ozone,  O3;  gas  flow  effects  for  all 
species;  and  axial  drift  effects  for  all  charged  species. 

The  general  differential  rate  equation  that  is  itera- 
tively solved  in  Chemco  for  each  of  the  species  is : 

tv.  Gains  due  to  collisions  - Losses 

- s = due  to  collisions  - Diffusion  (84) 

dt  Losses  - Flow  Losses  - Drift 
Losses 

Species  Included  in  Chemco.  The  inclusion  of  all  pos- 
sible neutral,  charged  and  excited  species  present  in  an  oxy- 
gen plasma  is  not  practical.  The  number  of  species  that  can 
be  tracked  is  limited  by  the  size  of  the  computer,  the  com- 
puter time  available,  and  the  current  state  of  knowledge  of 
the  species  reaction  rates.  Therefore,  only  the  ten  domi- 
nant neutral  species  and  important  ions  present  in  the  plasma 
are  tracked  in  Chemco:  0,  O2 , (^(a'A)*  O3,  e,  0~,  0~2 , 0~3, 
0+ , and  0+2 • 

It  is  well  known  that  the  oxygen  molecule  readily 
forms  negative  ions  by  attachment  processes.  The  negative 
molecular  ions  included  in  Chemco  are  O2  and  0-j.  The  nega- 
tive atomic  ion,  0 , is  most  readily  formed  by  dissociative 
attachment;  it  is  included  in  Chemco.  The  only  positive 
molecular  ion  included  is  0^.  The  observation  of  ionized 
ozone,  O3,  has  not  been  reported  in  the  literature  and  was 
therefore  expected  to  be  unimportant  in  the  ion  balance; 
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thus,  it  is  not  included  in  Chemco.  The  positive  atomic  ion, 
0+,  is  generated  by  direct  ionization  of  the  atom  and  by 
dissociative  ionization  of  the  molecule;  it  is  included. 
Neutral  atomic  oxygen  is  predominantly  produced  by  electron 
impact  dissociation  and  dissociative  attachment  in  the  dis- 
charge; it  is  included  in  the  model.  Neutral  ozone,  O3,  is 
produced  by  several  chemical  reactions  and  has  been  observed 
to  occur  in  the  discharge;  it  is  included  in  the  code.  The 
only  electronically  excited  state  included  is  the  metastable 
O^a'A);  it  is  produced  primarily  by  electron  collision 
with  ground  state  O2.  Finally,  the  code  includes  the  pri- 
mary constituents:  molecular  oxygen  and  electrons. 

Species  Gain  and  Loss  Processes.  The  gain  and  loss 
processes  incorporated  in  Chemco  for  each  of  the  included 
species  are  listed  in  Tables  III  through  XXII.  The  letter 
M in  the  tables  refers  to  an  unspecified  specie  in  a three- 
body  collisional  process.  For  reference  and  utilization  in 
future  studies  a complete  listing  of  all  the  included  and 
nonincluded  rates  and  their  source  references  is  presented 


in  Appendix  D. 


Process  Number/Symbol 
8 
9 
13 
18 
19 
22 

23 

24 
29 
39 

51 

52 

53 

54 

55 
88 
93 

114 

119 

124 

136 

138 


139 


Process 


Table  IV.  Atomic  Oxygen  (0)  Losses 
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Process  Number /Symbol 
1 
3 
5 
7 

16 

26 

28 

30 

31 

32 

33 

34 

35 
104 
123 
127 
132 
134 

LDO 


Process 

0 + e — +■  0+  + 2e 

0 + 0 — > O2  + e 

0 + e 0 

0 + O2  + e — + 0 + O2 

0 + O2  — ^ O3  + e 

0 + O2  — ^ 0 + O2 

20  + O2  — ^ 2O2 

0 + 0+  + M — *•  0 

20  + M ->  02  + M 
0 + 2O2  ■+  O3  + O2 
20  + 02  ->  O3  + 0 
0 + 03  202 

0 + O3  -+  O2  + O2 

0 + 03  °2  + °2 

0 + O3  ->  202  + e 

0 + O2  + e 0 + O2 

20  + O2  -*•  02  + O2 
20  + O2  ->  03  + 0 
Diffusion  Losses 
Flow  Losses 


+ CN 


Table  V.  Molecular  Oxygen  (O2)  Gains 


Process  Number/Symbol 

Process 

3 

0 +0  — > O2  4-  0 

6 

0 4-  0^  — ^ O2  4-  0 

9 

0 + 02  — *■  O2  + 0 

11 

o~  + o+  + 0 -*■  o2  + 0 

23 

o2  + 0+  -*•  02  + 0 

26 

O2  + 0 02  + 0 

28 

02  + 20  ->  202 

29 

°2  + °2  °2  + 20 

31 

20  + M — >■  02  + M 

34 

0 + 0 2 — >■  202 

35 

0 + O3  o2  + o2 

36 

O2  4-  O2  202  ^ G 

38 

02  + M + e 02+M 

43 

°2  + °2  202 

45 

°2  + °3  °2  + °3 

54 

G 4-  O3  O2  4-0  + 0 

55 

02  + 0^  202  + 0 

71 

02  + 0 2 202 

82 

202  O2  + O2 

94 

02  + 2e  02  + e 

100 

°2  °2 

104 

0 + o^  o2  + 02 

105 

e + 02  — ► 0 2 + e 
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Process 


O2  + 0 — > O2  + 0 
O2  + 0^  - > O2  + o3 
~ ^ 2O2  e 

02  + O3  — > 2O2  + 0 

03  + 0 — > 2O2  + G 

20  + O2  - ► O2  + O2 
Diffusion  Losses  of  0 
Diffusion  Losses  of  O2 
Diffusion  Losses  of  O2  (a'A) 
Flow  Losses  of  (^(a'A  ) 

Drift  Losses  of  0^ 

Drift  Losses  of  O2 
Flow  Losses  of  0 
Drift  Losses  of  0+ 

Flow  Losses  of  O3 
Diffusion  Losses  of  0 


+ CM 


Table  XI . Molecular  Oxygen  ( 0 2 ) Losses 


Process  Number/Symbol 

Process 

4 

o"  + o2 

03 

+ e 

12 

0“  + O2  + M 

-> 

03  + M 

16 

°2  + 0 °3 

; + 

e 

24 

O3  + 0+ 

°3 

+ 0 

32 

0 + 202 

°3 

+ o2 

33 

20  + 0 2 

_ — . * 

°3 

+ 0 

117 

0 + o2 

°3 

+ e 

134 

20  + O2  -*» 

°3 

+ 0 

ldo5 

Diffusion  Losses  of  0i 

i m 

1 

Drift  Losses 

; of  03 

t 

Table  VIII.  Ozone  (03) 


Process 


O3  + e — ^ 0 + O2 

O3  + 0 O3  + 0 

O3  + e ->■  O2  + 0 

0^  + 0 202 

03  + 02  ->  O3  + 02 
O3  + e -»•  02  + 0 + e 

O3  + O2  ~ ® 

0 ^ + 0 — >•  O2  + O2 
O3  + O2  ® 

Flow  Losses 
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Process  Number/Symbol 

Process 

1 

e + 0 — > 2e  + 0+ 

2 

s + O2  - 2e  + O2 

3 

0~  + 0 -*•  e + 02 

4 

0 + O2  ~ 6 + 0 2 

16 

O2  + 0 — >■  e + 0 ^ 

36 

O2  + O2  - ^ e + 2O2 

51 

e + 0 — 2e  + 0 

88 

e + 02  2e  + 0 + 0+ 

117 

0“  + 0*  -*»  e + 03 

118 

6 + 202 

123 

°3  + 0 e + 202 

126 

e + O2  ->  2e  + O2 

140 

e + O2  ->  2e  + 0 + 0+ 

i 

o 

J8 

* 

Drift  Losses  of  0~ 

lMC>2 

Drift  Losses  of  O2 

lMC>3 

i 

L 

Drift  Losses  of  O3 

Table  X.  Electron  (e)  Losses 


Process  Number/Symbol 

Process 

5 

e + 0 0" 

6 

e + 03  0"  + 02 

7 

e + 0 + 02  -*»  0"  + 02  ' 

18 

e + 03  02  + 0 

19 

e + 0^  20 

21 

e + 0 + 03  ~ ^ O2  + 0 

22 

e + 0+  + M->  0+M 

38 

e + 03  + M ->•  02  + M 

41 

e + 2O2  O2  + O2 

53 

e + 02  ->■  0 + 0“ 

93 

2e  + 0+  -»•  e + 0 

94 

26  + O2  6 + O2 

* — * 

127 

e + 0 + 02  0 + 02 

* — 

130 

e + 0 + O2  — ^ O2  + 0 

137 

e + 0 2 + 0 2 *■  0 2 + 0 2 

* — 

139 

e + O2  0 +0 

141 

e + 02  + 02  °2  + °2 

142 

e + 02  + 02  — ► O2  + 02 

LDe 

Diffusion  Losses 

+ 

o 

3 

Drift  Losses  of  0+ 

LM02 

Drift  Losses  of  0-> 
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Table  XI.  Atomic  Negative  Ion  (0~)  Gains 


Process  Number/Symbol 


Process 


e + 0 


e + 0- 


e + 0 + 0- 


02  + 0 
e + 0o 


e + 0 + 0' 


e + O' 


0+0' 


0 + 0- 


0“  + 0- 


0+0 


0 + 0 


0+0 


* (N 


Process 


18 

21 

35 

41 

114 

124 

130 

137 

141 

142 


e + 03  02  + 0 

e + 0 + 02  02  + 0 

0 + 03  ->•  02  + o2 

6 + 202  02  + 0 2 

0 + o2  ->  02  + 0 

0 + o2  02  + 0 

G + 0 + 02  — 02  + 0 
e + 02  + 02  — 0 2 + 02 

g + 02  + 02  02  + o2 

* _ 

e + 202  02  + 0 


ro  * 


Table  XIV.  Molecular  Negative  Ion  (O2)  Losses 


Process  Number/ Symbol 

Process 

16 

°2  + 0 — 

03  + e 

23 

°2  + 0+ 

o2  + 0 

26 

°2  + 0 

0"  + o2 

29 

02  + O2  — 

02  + 20 

36 

°2  + °2 

202  + e 

43 

°2  + °2 

202 

45 

°2  + °3 
- * 

°3  + °2 

118 

o2  + o2 

202  + e 

LDOi 

Diffusion  Losses 

LM02 

1 

Drift  Losses 

Lvoi 

Flow  Losses 

Table  XV.  Ozone  Negative  Ion  (O3)  Gains 


Table  XVIII 


Atomic  Positive  Ion  (0+)  Losses 


Process  Number/ Symbol 

Process 

8 

0+  + 0" 

20 

11 

0 

+ 

1 

o 

+ 

+ 

o 

-> 

°2  + 0 ' 

22 

0+  + e + M 

-> 

0 + M 

23 

0+  + O2  -> 

°2 

+ 0 

24 

o+  + O3 

°3 

+ 0 

30 

0+  + 0 + M 

02  + M 

39 

0+  + O2 

°2 

+ 0 

93 

0+  + 2e 

0 + e 

136 

0 + 02 

02 

+ 0 

lD0+ 

Diffusion  Losses 

lMO+ 

LVO+ 

Drift  Losses 

Flow  Losses 

Table  XX.  Molecular  Positive  Ion  (02)  Losses 


Process  Number/Symbol  Process 


9 

02 

+ 

0"  -*■ 

0 + o2 

12 

0 J 

+ 

0“  + M 

- ► 0 3 + M 

19 

°2 

+ 

e -=*■ 

20 

29 

°2 

+ 

0 2 

02  + 20 

38 

°2 

+ 

e + M 

— >■  0 2 + M 

43 

°2 

+ 

°2 

202 

94 

02 

+ 

2e  -> 

02  + e 

LDC>3 

Diffusion  Losses 

lmo2 

Drift  Losses 

LV02 

Flow 

Losses 

Table  XXII.  Molecular  Metastable  (O^Ca'A))  Losses 


Process  Number/Symbol  Process 


71 

* 

°2 

+ 02 

202 

82 

20 

2 °2 

+ o2 

100 

°2 

o2 

105 

°2 

+ e — > 

02  + e 

107 

* 

°2 

5-  0 -> 

02  + o 

108 

02 

+ 0 ^ ^ 

°2  + °3 

117 

02 

+ 0“ 

03  + e 

118 

oj 

+ °2  "*■ 

202  + e 

119 

°2 

+ °3 

202  + 0 

124 

o2* 

+ 0~  -> 

°2  + 0 

126 

* 

°2 

+ e -*• 

«2  + 2e 

130 

* 

°2 

+ 0 + e 

-*■  o2  + 0 

134 

* 

°2 

+ 20  -*• 

°3  + 0 

136 

* 

o2 

+ 0+ 

0+  + 0 

138 

* 

°2 

+ e ->■ 

20  + e 

139 

* 

02 

+ e 

o + o" 

140 

* 

°2 

+ e 

0+  + 0 + 2e 

141 

_ * 
°2 

+ 02  + e 

0 

+ 

1 CN 
O 

t 

142 

202  + e — ^ 

°2  + °2 

Ldq*  Diffusion  Losses 

LV02 


Flow  Losses 


Diffusion  Losses  of  the  Charged  Species.  In  a plasma 
in  which  there  is  a multiplicity  of  charged  species  of  dif- 
fering mobilities  and  diffusion  coefficients,  there  is  an 
initial  tendency  for  these  species  to  diffuse  out  of  the 
plasma  at  differing  rates.  The  charge  separation  that  re- 
sults sets  up  an  electric  field  which  retards  the  diffusion 
of  the  faster  diffusing  species  and  enhances  the  diffusion 
of  the  slower  species.  The  resulting  coupling,  which  tends 
to  cause  all  species  to  diffuse  at  some  intermediate  rate, 
is  termed  ambipolar  diffusion.  Classical  ambipolar  diffu- 
sion theory  treats  only  two  species:  electrons  and  a single 
specie  of  positive  ions.  For  the  oxygen  plasma  a multipli- 
city of  species  must  be  considered,  including  electrons  and 
both  negative  and  positive  ions:  e,  0”,  0^,  O3,  0+,  and  0^. 
A three-species  solution  to  the  coupled  ambipolar  diffusion 
equations,  as  derived  by  Thompson,  is  used  for  the  chemistry 
code*  (Ref  24:818-819).  In  order  to  utilize  Thompson's  so- 
lution, the  five  ionic  species  carried  in  Chemco  were  re- 
duced to  two  by  combining  all  of  the  negative  ions  into  an 
averaged  negative  ion  and  all  of  the  positive  ions  into  an 
averaged  positive  ion**. 


*For  this  study  Thompson's  derivation  was  accepted  as 
presented  in  the  reference.  Attempts  at  an  independent  de- 
rivation were  unsuccessful. 

**As  will  be  shown  later  the  mobilities  of  the  differ- 
ent ions  are  sufficiently  close  together  to  allow  this  ap- 
proximation to  be  reasonable. 
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For  any  species,  S,  the  radial  diffusion  equation  is 


L 


i 

f 


(Ref  8:187) : 

f7  = D5*VNs  + Ns'/V^  (85) 


where  IT  is  the  flux  density, 

Ds  is  the  diffusion  coefficient; 

Ns  is  the  number  density; 

/A Js is  the  mobility; 

E is  the  space  charge  electric  field. 

Charge  neutrality  within  the  plasma  requires  that: 

N-h  = N_  4 N e.  <86) 

while  at  the  boundary: 

E - r * re  <s7» 

Combining  85,  86,  and  87  so  as  to  eliminate  the  electric 
field  term: 

~-Das-  VNs  (88) 


where  Das  are  the  ambipolar  diffusion  coefficients,  defined 
as : 


m 


( 


Da-  = D+' 


fl  + Y-i  I****) 


(■#:)■  {_;+ 


r>  ■ n 


Ae 


I -V  "V  A-  2*  c^'  X 

7 N 


i i w ■ « — W 

I + A±^4  2^ 
>le  >*« 


©£  and  Y are  defined: 


os  = 


NLr 

Ne 


and: 


Y = 


ti 


(90) 


(91) 


(92) 


(93) 


where:  £e  is  the  electron  characteristic  energy  in  ev. 

£/,  is  the  ion  temperature  in  ev*. 


For  a T^  of  300°K: 


and: 


£?=  3.88X10'  o-v 


Y=  25.8*  £c 


(94) 


The 


total  positive  ion  density  consists  of  the  sum  of  the  0 


ions  are  assumed  to  be  at  the  same  temperature. 
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and  02  densities: 


N+  - N0*  N, 


(95) 


Similarly,  the  total  negative  ion  density  is  described  as: 


N-  - N0-+  N0;+  Ny 


(96) 


Since  ion  mobility  is  the  dominant  factor  in  determin- 
ing the  diffusion  rate,  it  is  the  factor  that  is  averaged  by 


weighting  the  mobilities  of  the  constituents  of  N+  and  N : 


U - yUp+'N o+  + Moim  Nfoz 


(97) 


— y^°~  ^o~  /^-o z * NI07  N/o3  (98) 

N_ 


For  each  of  the  basic  species  the  mobility  is  (Ref  49:6) 


A 


liO-T 


2.73'  P 


■J* 


OS 


(99) 


Where yLXos  is  the  reduced  mobility,  T and  P are  the  tempera- 
ture and  pressure  of  the  background  gas  in  °K  and  Torr  re- 


spectively. In  terms  of  the  number  density  of  the  background 
gas,  equation  99  becomes: 


19 


A-i^rja-vM 


os 


(100) 


• * • • 


r 
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where  N is  the  background  gas  density  in  CM  . For  oxygen 
ions  in  their  own  background  gas  (Ref  50:477-487): 


JJ-  oo~ 

= 3.20 

CM2 

Volt-1 

SEC-1 

= 2.16 

CM2 

Volt-1 

SEC-1 

/Loot 

= 2.24 

CM2 

Volt-1 

SEC-1 

/^oo3 

= 2.55 

CM2 

Volt-1 

SEC-1 

The  mobility  of  the  0+  ions  is  not  listed.  For  this  study 
it  was  assumed  to  have  the  same  value  as  the  0 ions*.  Com- 
bining the  reduced  mobility  values  with  equations  97,  98, 
and  9 9 : 

= 2^r  ~(3,2' zi4~Nv)  <ioi) 

i^p-5-(3.2-N0-+2.l&Noi+2.55Kloi)  (102) 

The  electron  mobility , JJ,  e , was  determined  by  the  Boltzmann 
equation  code  (EED)  and  was  read  into  Chemco  as  a function 
of  E/N. 

The  averaged  diffusion  coefficient  is  related  to  the 


*This  was  a reasonable  assumption  since  reduced  mobi- 
lity is  inversely  proportional  to  the  mass  of  the  ion  and 
its  collision  cross  section  with  the  background  gas,  both  of 
which  should  be  approximately  the  same  for  the  o-  and  0+  ions. 
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averaged  mobility  by  the  Einstein  relation  (Ref  8:186): 


Ds  - (103) 

For  ion  temperatures  of  300°K: 

= 2.53  X l02*yU±  (104) 

The  diffusion  loss  rate  for  the  averaged  charged  spe- 
cies is  a function  of  the  ambipolar  diffusion  coefficients 
calculated  in  the  previous  paragraphs  and  of  the  radial  den- 
sity profiles  of  the  species.  For  the  averaged  ionic  spe- 
cies the  radial  density  profiles  were  assumed  to  be  zero  or- 
der Bessel  functions.  The  assumption  matches  closely  the 
experimental  results  reported  by  Thompson  (Ref  20:514). 

NtM  - Not  ‘ <105> 

where  No*,  is  the  centerline  density.  The  diffusion  loss,LDS, 
at  the  wall,  per  unit  length,  was  calculated  from: 

t 2 v 

U£  = (106) 

o e 
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The  loss  rates  for  the  ionic  species  were  assumed  to  be  re- 
lated to  that  of  the  averaged  ions  by: 


* 


The  averaged  loss  rate  from  the  volume  is: 


^ Lx^s  -rr-.D2  'Lds 


TT-R1 


Combining  equations  100,  109,  111,  112  and  113: 


<u>,  • 


<L.>_.  • “g  • »-■  <•», 


(113) 


(114) 


(115) 


(116) 


<U>o; 


(117) 


<L»>0-  =Jfrir!-^--<N>0- 


(118) 


The  electron  radial  density  profile  varies  from  generally 
Bessel  for  small  values  of  ©C  to  uniform  with  a rapid  drop- 
off near  the  wall  for  large  values  of  e<  (Ref  24:819-820). 
This  functional  relationship  in  the  vicinity  of  the  walls 
cannot  be  approximated  analytically.  Therefore,  the  electron 


176 


Mrfiibi 


diffusion  losses  are  determined  by  assuming  a zero  net  cur- 
rent flow  to  the  wall. 

<U)e  = <U>0  + <u>c-  <U>o_-  <u>of  <U>OJ  (119) 

Diffusion  Losses  of  the  Neutral  Species.  The  net  flux 
of  the  neutral  species  to  the  wall  is  (Ref  51:19): 

IT1  — R.w  s * Ns'  (120) 

•s  - 4-0- 


where  Rws  is  the  deactivation  or  recombination  coefficient 
and  is  the  mean  speed  of  the  species. 


(121) 


The  average  loss  per  unit  volume  is  related  to  the  net  flux 
by: 


__  2 •TP  Ro  . p 
ir-R  2 ** 

For  O2 (a ’ ^ ) at  300°K*: 

Rwor  =2X10^  sec"1 


(122) 


♦Refer  to  process  number  73  in  Appendix  D for  source 
reference. 
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<c>0,  = 4.4  6xio4  % 

Therefore:  = 0.4£>9*Mot.  (123) 

°2.  2 

For  0 at  300°K: 

l?wo  “ lO*  S&c  ' * 

<c>0  = &.3XI04  §fc 

Therefore:  <(  [_o\  - 3.32*  Mo  (124) 

Flow  and  Drift  Losses.  Since  the  gas  was  flowing 
through  the  tube  during  the  experiment,  there  was  a loss  of 
all  the  species  due  to  the  flow.  Also,  the  ionic  species 
drifted  in  the  electric  field  to  one  of  the  electrodes  where 
they  were  neutralized.  At  any  incremental  volume  in  the  po- 
sitive column  the  actual  gain  or  loss  of  a specie  due  to 
flow  or  drift  is  the  net  between  the  number  of  that  specie 
entering  the  volume  and  the  number  leaving.  This  value  is 
impossible  to  calculate  without  knowing  the  axial  variation 
of  the  species  densities  and  the  electric  field.  For  a 


*Refer  to  process  number  75  in  Appendix  D for  source 
reference. 
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perfectly  uniform  positive  column  the  losses  would  be  zero 
except  at  the  anode  and  cathode  boundaries  of  the  column. 

For  this  study  the  overall  flow  and  drift  losses  for  the 
discharge  are  calculated  and  are  assumed  to  be  distributed 
throughout  the  volume  of  the  positive  column.  This  assump- 
tion represents  an  upper  limit  on  these  losses  in  the  column 
Assuming  a gas  temperature  of  300°K,  the  experimental 
flow  rate  was  kept  constant  at  31  cm/sec.  This  equates  to 
a flow  loss  rate,  per  unit  volume,  of: 

LPS=^-NS  (125) 

For  a column  length  of  50  cm: 

Lfs  = o.a-Ns  <126> 

The  loss  of  each  of  the  species  is  returned  to  the  plasma 
as  molecular  oxygen  in  order  to  conserve  mass.  This  corre- 
sponds to  a replacement  of  the  species  flowing  out  of  the 
downstream  end  of  the  column  by  gas  flowing  in  at  the  up- 
stream end. 

The  drift  velocity  of  a specie  S,  Vds*  is  related  to 
the  electric  field  by: 


Vt>s  £■ 


(127) 
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The  mobility  of  a specie  S,^s  , 


is  related  to  its  re- 


duced mobility , yjios  , by: 


„ 2.6,9X10 

A n: 


Combining  equations  127  and  128: 


19  £ 

a II  - - * — 


(128) 


Vds=  2.C9Xio'yuos- 


(129) 


Using  the  reduced  mobilities  listed  in  Reference  50, 


'v'4d  B.  6 X-  IO  • ^ 


V„0.-  - s-.8xio,9-£ 


V»ot  = 6.02  Xio".  E 


Vjoj  - 6.85x10s.  | 


(130) 


(131) 


(132) 


(133) 


The  drift  loss  rate  per  unit  volume  of  the  ionic  species  is: 


U - ^ - Ws 


(134) 


^MS  ^50  * 


(135) 


*#g£»  ;■* . i ■ 


Actually,  in  Chemco  the  drift  and  flow  losses  are  com- 
bined for  the  ionic  species.  Due  to  experimental  geometry 
in  which  the  cathode  was  at  the  downstream  end  of  the  column, 
the  net  loss  to  the  positive  ions  is  the  sum  of  the  flow  and 
drift  losses;  the  net  loss  of  the  negative  ions  is  the  dif- 
ference between  the  flow  and  drift  losses. 

Results  of  Calculations 

Electron  Density.  Chemco  was  run  for  initial  molecu- 

1 (z  17 

lar  densities  between  3.2  x 10  and  5.6  x 10  oxygen  mole- 
cules per  cubic  centimeter  and  for  reduced  field  values, 

(E/N) , between  14  and  100  Townsends.  The  electron  densities 
determined  by  these  calculations  are  depicted  in  Figures  64 
and  65.  It  should  be  emphasized  that  any  particular  code 
calculation  was  performed  at  a fixed  value  of  E/N  and  is 
therefore  applicable  only  to  the  stable  form  of  the  discharge, 
the  high  field  form.  The  rapid  decrease  of  electron  density 
with  E/N,  as  illustrated  in  figure  64,  would  appear  to  be 
the  transition  of  the  discharge  from  a stable  to  a nonself- 
sustaining  form.  This  corresponds  to  the  experimentally  ob- 
served transition  between  the  high  field  and  low  field  forms 
of  the  discharge. 

Figures  64  and  65  represent  calculations  in  which  two- 
step  ionization  through  the  O2 (a^A  ) level  was  included.  In 
order  to  determine  the  magnitude  of  the  effect  of  the  two- 
step  ionization  on  the  electron  density,  a comparison  calcu- 
lation was  made  in  which  the  population  of  the  O2 (a^A  ) 
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level  was  set  at  zero.  This  calculation  was  made  for  an 
initial  molecular  density  of  1 x 10^  oxygen  molecules  per 
cubic  centimeter.  The  results  of  this  calculation  are  de- 
picted along  with  the  two-step  ionization  results  in  figure 
66.  Over  the  range  of  the  experimental  results  (40  to  60 
TD) , it  was  found  that  two-step  ionization  is  an  important 
process  in  determining  the  electron  density.  At  higher 
values  of  E/N,  the  two-step  process  becomes  less  important 

ie  1 

due  to  a decrease  in  the  (^(a  ^)  and  due  to  the  lower  ioni- 
zation threshold  of  the  metastable  becoming  less  important. 

Figures  67  and  68  compare  the  results  of  Chemco  with 
experiments  and  with  electron  current  continuity  calcula- 
tions. The  experimental  results  were  obtained  from  the  mi- 
crowave cavity  experiments  reported  in  Chapter  VI.  The  cur- 
rent continuity  calculations  were  determined  from  equation 
69.  The  chemistry  code  values  were  calculated  using  an  E/N 
taken  from  experiments.  In  determining  the  experimental  E/N, 
the  electric  field  was  measured  directly;  the  number  density 
was  calculated  by  determining  the  average  of  the  radial  tem- 
perature profile  (equation  22)  and  using  the  results  in 
equation  28; 


9-7 x io<0-  P 

( + 2.9o) 


(136) 
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Figure  66.  Effect  of  0o*(a'A  ) on  Electron  Density 
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where  Ale  is  the  temperature  difference  between  the  tube 
wall  and  the  center  of  the  discharge,  in  degrees  centigrade*. 
The  results  of  Chemco  were  interpolated  or  extrapolated  us- 
ing the  average  N and  E/N  to  determine  an  average  electron 
density.  This  average  density  was  then  multiplied  by  1.64 
to  obtain  the  centerline  electron  density;  1.64  is  the  peak 
to  average  ratio  for  a Bessel  profile. 

The  comparison  of  the  results  of  Chemco  to  other  cal- 
culations and  to  experiment  is  limited  by  two  factors.  A 
primary  limitation  is  the  comparison  of  a one  dimensional 
time  or  axial  variation  code  to  a two-dimensional  time  or 
axial  variation  and  radial  variation  experiment**.  The  code 
calculations  were  for  values  of  E/N  which  do  not  vary  in 
position  or  time.  The  stable- field  experimental  data  did 
not  change  as  a function  of  time,  but  it  is  expected  that 
the  quantities  are  functions  of  radial  position.  The  second 
limitation  is  the  linking  of  the  results  of  the  code  and  ex- 
periment through  the  reduced  field,  E/N.  The  experimental 
evaluation  of  this  critical  parameter  required  the  accurate 
determination  of  the  temperature  of  the  gas  in  the  discharge 
plasma.  As  discussed  in  Chapter  IV,  this  measurement  was 
difficult;  and  for  the  comparison  regime,  the  values  used 
were  extrapolations  of  measurements  made  at  lower  values  of 


♦Equation  24  withATc  = To-Tw. 

**The  plasma  parameters  were  assumed  to  be  symmetric 
around  the  tube  axis. 
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input  power.  Despite  these  approximations  the  results  of 
the  Chemco  calculations  compare  well  with  the  results  from 
electron  current  continuity  calculations  and  with  the  micro- 
wave  experimental  data. 

Ion  Densities.  Figures  69,  70,  and  71  depict  the  cal- 
culated charged  species  densities  as  functions  of  reduced 
electric  field  for  three  different  values  of  initial  molecu- 
lar oxygen  density.  For  the  experimental  parameter  space 
of  this  study,  the  dominant  positive  ion  was  calculated  to 
be  O^.  This  result  agrees  with  the  experimental  observa- 
tions of  Thov  ’ison  and  Keren  (Ref  12:52,  Ref  19:8-9,  Ref  24: 
820).  The  dominant  negative  ion  is  0~,  as  previously  mea- 
sured by  Thompson,  Lunt , and  Sabadil  (Ref  24:820,  Ref  25: 
1068,  Ref  26:299-308).  However  as  the  value  of  E/N  is  low- 
ered toward  the  point  of  transition  between  the  stable  and 
unstable  forms,  the  negative  molecular  ions,  O2  and  O3,  be- 
come significant,  especially  as  the  neutral  number  density 
is  increased.  In  fact,  at  a molecular  density  of  3.16  x 
1017  molecules  per  cubic  centimeter,  the  negative  molecular 
densities  are  approximately  the  same  as  the  atomic  ion  at 
the  transition  point. 

Neutral  Densities.  Figures  72,  73,  and  74  depict  the 
calculated  neutral  species  densities  as  functions  of  E/N  at 
three  different  initial  molecular  densities.  Of  note  is  the 
relatively  low  E/N  at  which  the  02(aiA)  density  reaches  a 
maximum.  At  E/Ns  above  this  value  the  density  of  O2  (a^A ) 
decreases,  reducing  the  importance  of  two-step  ionization 
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at  high  E/N  values.  The  O3  density  was  calculated  to  be 
several  orders  of  magnitude  lower  than  that  of  the  other 
neutral  species.  This  is  consistent  with  the  mass  spectro- 
meter experiments  during  which  O3  was  not  observed  within 
the  two  orders  of  magnitude  range  of  the  neutral  spectra 
plots* . 

The  change  of  atomic  oxygen  density  with  discharge 
tube  current  is  compared  with  experiment  in  figure  75.  An 
absolute  calibration  of  the  atomic  oxygen  density  was  not 
made  during  the  mass  spectrometry  experiments.  However,  a 
calibration  using  the  results  of  Chemco  indicated  a maximum 
dissociation  of  20%  for  the  5 torr  experimental  run. 

Principal  Electron  Processes.  Figures  76  through  81 
depict  the  principal  electron  gain  and  loss  processes  at 
three  different  initial  molecular  oxygen  densities.  The 
rates  are  normalized  to  either  the  total  gain  or  total  loss 
rate  to  show  relative  importance.  The  principal  gains  for 
the  electrons  are  from  detachment  processes  from  the  atomic 
and  molecular  negative  ions.  The  major  detaching  specie  is 
the  O^fa^A)  excited  molecule.  Atomic  oxygen  plays  a minor 
role  as  a detacher**.  The  principal  loss  mechanisms  for  the 
electrons  are  dissociative  attachment  from  the  O2 (a^A  ) and 


*Reference  Chapter  VII. 

**Chemco  runs  in  which  the  O2  (a-^A  ) population  is  set 
at  zero  results  in  the  same  high  to  low  transition  phenome- 
non. However,  in  those  cases,  the  detaching  specie  is  only 
atomic  oxygen;  and  the  transition  occurs  at  a much  higher 
value  of  E/N,  as  depicted  in  figure  84. 
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Figure  81.  Principal  Electron  Losses  at  N=3. 16X10  cm 


ground  states  of  the  molecule. 

As  the  point  of  transition  is  approached  by  the  lower- 
ing of  the  reduced  field,  the  relative  effect  of  direct  mole- 
cular ionization  from  the  O^Ca^A)  and  ground  states  becomes 
numerically  insignificant  to  the  other  processes.  The  at- 
tachment/detachment processes  numerically  dominate  the  elec- 
tron losses  and  gains.  However,  these  processes  merely 
transfer  electrons  back  and  forth  from  bound  to  unbound 
states;  and  any  true  losses  to  other  processes  are  balanced 
by  ionization.  Therefore,  although  numerically  insignifi- 
cant near  transition,  the  ionization  process  is  key  to  the 
maintenance  of  the  stable  discharge.  As  the  ionization  pro- 
cess rapidly  decreases  near  the  point  of  transition,  the 
electrons  stored  in  the  attachment/detachment  processes  be- 
come depleted  through  other  loss  processes;  and  the  discharge 
becomes  nonself-sustaining.  Due  to  the  presence  of  negative 
ions,  the  diffusion  process  is  insignificant  with  respect 
to  the  attachment  processes  in  the  range  of  initial  molecu- 
lar densities  considered  in  Chemco.  However,  the  diffusion 
loss  is  comparable  to  the  recombination  and  flow  losses. 
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X. 


Discussion  and  Conclusions 


Chapters  III  through  VII  of  this  study  reported  on  the 
results  of  the  experimental  observations  of  the  oxygen  dis- 
charge. Chapters  VIII  and  IX  presented  a computer  model  of 
the  discharge  and  compared  that  model  with  the  experiment's. 

In  this  chapter  the  results  of  these  previous  chapters  are 
used  to  explain  the  phenomena  associated  with  the  oxygen 
plasma.  Specifically,  the  discussion  focuses  on  the  high 
field  form  of  the  discharge  and  the  transition  from  the  high 
field  form  to  the  low  field  form.  Quantitative  data  are 
used  to  draw  conclusions  regarding  these  topics.  Some  qua- 
litative discussion  and  conclusions  are  also  presented  re- 
garding the  low  field  form  and  the  transition  from  the  low 
field  form  to  the  high  field  form. 

It  is  believed  that  the  results  developed  in  this  study 
are  much  more  comprehensive  than  those  previously  developed 
for  the  oxygen  plasma.  However,  there  still  remains  a sig- 
nificant amount  of  work  that  must  be  accomplished  to  fully 
understand  the  phenomena.  In  addition  to  oxygen,  it  is  be- 
lieved that  the  experimental  and  modeling  approaches  pre- 
sented in  this  study  are  also  useful  in  the  investigation 
of  other  electro-negative  plasmas. 


Discussion 

High  Field  Form.  Throughout  the  region  of  the  exper- 
imental investigation  the  high  field  form,  though  stable. 


"■r  • 


is  nevertheless  a plasma  in  which  the  negative  ions  play  im- 
portant roles.  As  the  point  of  transition  between  stable 
and  unstable  forms  is  approached,  the  attachment-detachment 
processes  dominate  the  electron  gains  and  losses,  leading 
to  a situation  in  which  attachment  is  approximately  balanced 
by  detachment  (Ref  30:11). 

Ne*  ka  - N_*  Kj  (137) 

where  Ka  is  the  attachment  processes,  Kd  is  the  detachment 
processes.  The  dominant*  attachment  and  detachment  processes 
for  values  of  E/N  between  60  Td  and  the  transition  points 
and  for  initial  molecular  densities  between  3.16  x 1016  and 
3.17  x 1017cm~3  are  listed  in  table  XXIII**.  Equation  137 
was  satisfied  very  closely  throughout  the  E/N  range  of  this 
investigation  for  the  processes  of  table  XXIII. 

Although  not  the  primary  source  of  electrons,  ioniza- 
tion is  important  for  replenishing  the  electrons  lost  by 
processes  other  than  attachment.  For  a stable  discharge, 
ionization  balances  these  other  losses. 

Kl-N0;Ne  - K/  N+*  NJe  4-  Ne+  Kv*Ne  <138) 


*Within  two  orders  of  magnitude. 

**Refer  to  figures  76  through  81  for  a graphical  des- 
cription of  the  dominant  processes.  Equation  137  is  true  only 
when  the  attachers  and  detachers  are  the  same  species:  in 
this  case  O2 (al^)  dominates  close  to  the  point  of  transition. 
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where  Ki  is  the  ionization  coefficient,  Kr  is  the  recombina 
tion  coefficient,  Kdf  is  the  diffusion  loss  coefficient, 
and  Kv  is  the  drift  loss  coefficient.  The  dominant*  ioniza 
tion  and  nonattachment  loss  processes  are  listed  in  table 
XXIV. 

Table  XXIII.  Dominant  Attachment  and  Detachment  Processes 


Attachment 

Process 

Number* 

Func. 

Dep. 

(30 

Typical 

Rates 
to  60  TD) 

e + °2  “*■ 

0 + 0" 

53 

E/N 

1.7 

3.2 

x 10"11  t 
x 10-11  to 

* 

G + 0 2 — ^ 

0 + 0" 

139 

E/N 

3.1 

4.2 

x l0"}}  to 
x 10-H 

e + 2O2  ~ ^ 

°2  + °2 

41 

None 

1.9 

x 10"30 

e + O2  + 0^ 

-=*■  O2  + 

o2 

141 

None 

1.9 

x 10-30 

Detachment 

e + 0 — *■ 

e + 02 

3 

None 

1.9 

x 10_1° 

0 + 02  -*• 

03  + e 

117 

None 

3 x 

10"1° 

°2  + 0 — ' *■ 

03  + e 

16 

None 

3.3 

x IQ"10 

°2  + 

202  + e 

118 

None 

2 x 

10-10 

‘Processes  listed  in 
N = 1 x 10l^cm-3. 

descending  order  of  importance  at 

*Withm  two  orders  of  magnitude 


Table  XXIV.  Dominant  Ionization  and  Loss  Processes 


Ionization 

Process 

Func. 

Typical 

Rates 

Number* 

Dep. 

(30  to  60  TD) 

* + 

e + O2  — C>2  + 2e 

126 

E/N 

1.4  x 10-33 

7.5  x 10“13 

to 

e + O2  — O2  + 2e 

2 

E/N 

1.6  x 10-15 

2.7  x 10-13 

to 

e + 0 -v  0+  + 2e 

1 

E/N 

7.8  x 10"18 

6.8  x 10"14 

to 

N_  + Electrode  N + e 

Drift 

E/N 

- 

Non-Attachment  Losses 

e + O2  20 

19 

None 

2.1  x 10~7 

e + 0^  + M -*►  02+M 

38 

None 

1 x 10~26 

e + N+  + Wall  N 

Diffusion 

Ne/N 

- 

e + N+  + Electrode  -*•  N 

Drift 

E/N 

- 

♦Processes  listed  in  descending  order  of  importance  at 

N = 1 x 10^7cm-3. 

Using  the  processes  of  table  XXIV,  equation  138  was  satis- 
fied within  t 16%  over  the  range  of  E/N  values  from  60  Td 
to  the  points  of  transition  and  for  initial  molecular  densi- 
ties between  3.16  x 1016  and  3.16  x 1017cm_3. 
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unstable  when  two  conditions  are  satisfied  (Ref  30:23-26, 
Ref  52:1030-1031).  The  first  condition  requires  that: 


Ur  s l/.  JLLN{kO_ 

? 3LN{%} 


(139) 


This  is  a necessary  condition  and  requires  that  the  attach- 
ment coefficient  be  an  increasing  function  of  E/N  and  of  a 
magnitude  comparable  to  the  ionization  coefficient.  This 
relationship  leads  to  a net  loss  of  electrons  during  a posi- 
tive fluctuation  in  E/N*.  In  order  to  maintain  current  con- 
tinuity in  such  a plasma,  the  discharge  must  adjust  itself 
to  a higher  E/N  which  further  decreases  the  electron  density. 
It  is  this  feedback  mechanism  that  leads  to  instability  and 
runaway.  Figure  82  depicts  the  ratio  of  the  two  quantities 
in  equation  139  as  a function  of  E/N  and  for  the  two  limit- 
ing conditions  where  all  of  the  ionization  is  from  the 

-ft  -I 

02(axA)  level  or  from  the  ground  level.  The  data  was  taken 
from  the  Boltzmann  Code  results  of  Chapter  VIII.  Even  with 

* "I 

a highly  populated  02(a  &)  level,  the  necessary  condition 
for  instability  is  only  satisfied  for  E/N  values  below 
67.4  TD. 

The  second  condition,  a sufficient  condition,  requires 
that  the  negative  ion  concentration  be  significant  compared 


*From  an  experimental  viewpoint  in  order  for  the  elec- 
tric field  to  change  as  a function  of  the  changing  nature 
of  the  plasma,  the  discharge  must  be  connected  to  a constant 
Current  source. 
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to  that  of  the  electrons.  This  is  required  so  that  negative 
ion  and  electron  processes  can  couple  during  a disturbance 
in  plasma  properties.  Specifically,  this  requires  for  in- 
stability that: 


M_  > Ne  <140> 

where  it  is  assumed  that  in  an  electronegative  multi-compo- 
nent plasma,  a summation  can  be  performed  over  all  the  nega- 
tive ion  species: 

N_  = N0-  4-  N0-  + N0-  (141) 

Figure  83  depicts  the  ratio  of  these  parameters  as  functions 
of  E/N  and  initial  molecular  density.  The  data  was  taken 
from  the  chemistry  code  results  of  Chapter  IX.  Using  the 
equality  of  the  electrons  and  negative  ions  as  the  suffi- 
cient condition  for  transition  from  stability  to  instability, 
the  E/N  values  at  the  transition  point  are  plotted  on  fig- 
ure 84  as  a function  of  initial  molecular  density.  Also 
plotted  on  this  figure  are  the  data  points  taken  from  the 
experiments.  There  is  considerable  scatter  in  the  experi- 
mental values,  due  to  difficulty  in  determining  the  actual 
gas  density.  Even  so,  the  data  indicates  that  the  calcu- 
lated transition  points  are  in  reasonable  agreement  with  ex- 
periments . 

Pressure  Effects.  Since  the  more  critical  condition 


211 


40 


60 

E/ N ( Td  ) 


80 


100 


Figure  83.  Electron-to-Negative  Ion  Ratio 
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for  the  transition  between  the  stable  high  field  forms  to 
the  unstable  low  form  of  the  discharge  is  the  equivalence 
of  the  electron  and  negative  ion  densities,  it  is  appropriate 
to  study  this  relationship  as  a function  of  initial  molecu- 
lar density  (pressure).  As  depicted  in  figure  84,  the  tran- 
sition value  of  E/N  is  a decreasing  function  of  the  initial 
molecular  density.  This  implies  that  the  pressure  effect 
is  caused  by  processes  which  increase  the  electron  density 
or  decrease  the  negative  ion  densities  as  the  initial  mole- 
cular density  is  increased.  Normally,  this  results  from 
three-body  collisional  processes;  however,  there  are  no  do- 
minant three-body  processes  in  Chemco  which  could  cause  the 
effect.  There  are,  however,  two  noncollisional  processes 
which  can  account  for  the  pressure  effects  on  the  transition 
values  - the  diffusion  and  drift  losses  of  the  electrons. 

In  order  to  determine  the  effects  of  electron  diffu- 
sion and  drift  losses,  consider  an  expansion  of  equation 
138: 


Ki-N0-Ne=  Kr-Ne-N0,  + +■  Kv*  Me 


(142) 


where  NeN02  is  substituted  for  N+  and  Kr  takes  on  a new 
✓ 

value,  -Kr.  Equation  142  can  be  simplified  to: 


M JO  _ Kv  _ k?F 

Kr  Kr"  Wo2  Kp/*  Mo^ 


(143) 
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The  only  pressure  sensitive  terms  in  equation  143  are  the 
diffusion  and  drift  terms.  As  the  pressure  increases,  these 
terms  decrease;  and  the  transition  point  values  of  E/N  de- 
crease. In  order  to  determine  the  magnitude  of  the  effect 
of  these  terms  on  the  transition  point,  Chemco  was  modified 
to  neglect  the  electron  diffusion*  and  to  compensate  for  the 


pressure  effects  of  the  electron  drift**.  The  results  of 
these  modifications  are  illustrated  in  figure  84.  As  shown, 


C 


the  electron  diffusion  only  has  a minor  effect  on  the  tran- 
sition point  pressure  dependence;  the  major  effect  is  due 
to  the  drift  loss.  In  fact,  the  particular  compensation 
term  that  was  used  overcorrected  the  functional  dependence 
and  i esulted  in  a change  in  the  slope  of  the  curve.  No  at- 
tempt was  made  to  determine  a better  compensation  term  since 
the  objective  was  only  to  determine  the  pressure  dependent 

processes  in  Chemco. 

* 1 

Effect  of  0?(a  A)«  For  many  years  one  of  the  unknown 
aspects  of  the  oxygen  discharge  has  been  the  role  O^fa^A) 
plays  in  determining  the  characteristics  of  the  discharge. 
Although  no  direct  measurement  was  made  of  the  density  of 
O2 (a*A)  during  the  experimental  phase  of  this  investigation***, 

♦Actually  all  of  the  charged  particle  diffusion  terms 
were  neglected.  The  negative  ion  diffusion  losses  are  so 
low  as  to  have  negligible  effect  on  the  electron  to  negative 
ion  ratio. 

♦‘Compensation  term:  V:  ' = - • Vs 

* . 3.16  x lOl? 

♦♦♦The  (^(a-^A)  densities  that  were  predicted  by  the  cal- 
culations were  approximately  equal  to  those  reported  by 
other  experimenters  for  similar  discharge  conditions  (Ref  53)  . 


I 

I 
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a comparison  between  the  experimentally  observed  character- 
istics of  the  discharge  and  the  results  of  the  model  of  the 

ic  -I  * 1 

discharge  which  includes  O2  (a-'-A)  indicates  that  O2  (a  &) 
plays  a prominent  role  in  the  discharge.  It  affects  the 
characteristics  of  the  stable  high  field  form  in  two  differ- 
ent ways:  first,  two-step  ionization  is  required*  to  achieve 
the  electron  densities  observed  during  the  experiments;  and 

^ -1 

second,  the  detaching  effect  of  O2  (a-xS)  is  required  to  match 
the  high-to-low  E/N  transition  values  seen  in  the  experi- 
ments. 

Low  Field  Form.  The  unstable  low  field  form  of  the 
discharge  was  studied  experimentally  in  this  investigation 
but  was  not  modeled.  A considerable  expansion  of  the  che- 
mistry code  would  have  been  required  to  allow  a time  varying 
value  of  E/N  and  a fixed  electrical  current.  The  results 
of  the  high  field  model  can  be  used,  however,  to  draw  infer- 
ences as  to  the  processes  involved  in  the  low  field  form. 

At  the  point  of  transition  between  forms,  it  was  shown  that 
the  attachment  and  detachment  processes  are  very  strong,  and 
indeed  dominated  the  gain  and  loss  of  the  electrons.  The 
logarithmic  derivative  of  the  attachment  coefficient  with 
respect  to  the  reduced  electric  field  is  greater  than  that 
of  the  ionization  coefficient.  It  is  reasonable  to  assume, 
therefore,  that  as  E/N  is  further  decreased  and  the  discharge 
becomes  unstable,  the  instability  can  be  characterized  as 

— 

*In  addition  to  direct  ionization. 
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an  ionization-attachment-detachment  instability  (Ref  54). 

In  such  an  instability,  as  a local  variation  in  E/N 
takes  place,  the  strong  attachment  rate*  causes  a local  de- 
crease in  electron  density.  In  order  to  maintain  a constant 
discharge  current,  the  value  of  E/N  must  decrease.  This  pro- 
cess continues  until  either  the  ionization  rate  dominates 
or  until  there  is  a sufficient  buildup  in  detaching  species 
to  recover  the  electrons  through  the  detachment  process. 

The  study  of  the  stable  discharge  showed  that  the  latter 
process  occurs  at  a lower  value  of  E/N  than  the  former. 

Once  the  detachers  are  produced,  they  have  a rather  long 
life**;  and  during  the  time  they  are  removing  free  electrons 
from  the  negative  ions,  the  reduced  field  returns  to  a low 
value.  Once  the  electrons  are  detached  and/or  the  detachers 
have  decayed  or  have  chemically  combined,  the  process  starts 
over  again. 

Low-to-High  Field  Transition.  The  transition  from  the 
unstable  low  field  form  to  the  stable  high  field  form  is 
postulated  to  occur  when  the  detacher  specie  density  has  in- 
creased to  the  point  where  the  electron  density  is  contin- 
uously greater  than  the  negative  ion  density.  The  fact  that 
the  low-to-high  and  high-to-low  transitions  exhibit  hyster- 
isis  phenomena  indicates  that  different  physical  processes 
could  be  involved  in  triggering  the  transitions.  For  the 


♦Compared  to  ionization. 
♦♦Especially  the  O2 (a*A) • 
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low-to-high  transition,  a preliminary  study  (Ref  55)  of  the 
changes  to  the  electron  energy  distribution  function  that 
occurs  with  dissociation  indicates  a shifting  of  energy  into 
the  high  energy  "tail"  of  the  distribution.  This  was  stu- 
died by  inserting  cross  sections  for  atomic  oxygen  into  the 
Boltzmann  Equation  Code  and  running  the  code  for  various 
percentages  of  dissociation.  A study  of  the  forward  rates 
that  were  developed  from  the  code  runs  indicated  that  disso- 
ciation could  change  the  electron  energy  distribution  func- 
tion sufficiently  to  be  initially  self-enhancing.  This  pro- 
cess, therefore,  could  provide  a positive  feedback  for  the 
generation  of  atomic  oxygen.  Atomic  oxygen  is  a good  de- 
tacher but  is  not  considered  to  be  the  prime  detacher  near 

"k  -i 

the  point  of  transition;  (^(a^A)  is.  At  the  transition  val- 
ues of  E/N  the  effect  of  the  changing  electron  energy  dis- 
tribution function  on  the  generation  of  O2  (a-^A)  is  small. 

Conclusions 

The  following  are  the  conclusions  from  the  experimen- 
tal and  computational  portions  of  the  study.  For  the  most 
part,  this  summary  is  qualitative  only.  Since  the  quantita- 
tive results  are  generally  dependent  on  many  variables,  the 
reader  should  refer  to  the  appropriate  portion  of  the  study 
report  for  the  numerical  results. 

A summary  of  the  major  conclusions  is  presented  first. 
This  is  followed  by  a more  detailed  listing  of  the  results 
derived  from  the  experimental  and  modeling  phases  of  the 
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study.  All  significant  observations  are  listed.  Some  of 
these  have  been  previously  reported  by  other  investigators 
and  are  listed  here  as  corroboration.  Other  results  are 
original  with  this  study  and  are  identified  by  asterisks. 

Summary  of  Major  Conclusions.  The  primary  emphasis 
of  this  study  was  the  development  and  experimental  verifica- 
tion of  a model  of  the  stable,  high-field  form  of  the  oxygen 
discharge.  This  model  was  then  used  to  determine  the  char- 
acteristics of  the  stable  discharge  and  the  criteria  for 
transition  to  the  unstable,  low-field  form  of  the  discharge. 

The  chemistry  of  the  stable  high-field  form  of  the 
oxygen  discharge  is  dominated  by  attachment  and  detachment 
processes.  Electrons  are  continuously  attached  to  and  de- 
tached from  the  negative  ions.  The  principle  negative  ion 
is  0";  however,  the  molecular  negative  ions  C>2  and  O3  become 
significant  as  the  pressure  is  increased.  Both  the  metasta- 
ble molecule,  O2 (a^A)  , and  the  atom,  0,  act  as  electron  de- 
tachers. The  former,  however,  was  determined  to  be  the  pre- 
dominant detacher.  Although  the  negative  ion  densities  are 
high  in  the  stable  discharge,  their  combined  density  is 
smaller  than  that  of  the  electrons.  With  the  exception  of 
very  high  values  of  E/N,  the  positive  molecular  ion,  0^, 
dominates  the  atomic  ion,  0+. 

Although  attachment/detachment  processes  dominate  the 
electron  chemistry,  they  tend  to  only  transfer  the  electrons 
from  bound  to  free  states;  any  permanent  electron  losses 
must  be  balanced  by  ionization.  In  the  stable  discharge. 


the  ionization  processes  are  sufficiently  strong  to  dominate 
recombination  and  other  nonattachment  losses.  Two-step  ion- 
ization from  the  ground  level,  O2 (X^£)  , through  the  metasta- 
ble level,  O2  (a^-A)  / of  the  molecule  is  an  important  source 
of  electrons  in  the  stable  discharge. 

The  transition  from  the  stable  to  the  unstable  form 


of  the  discharge  occurs  when  the  reduced  electric  field 
value,  E/N,  is  lowered  to  the  point  where  the  logarithmic 
derivative  of  attachment  with  respect  to  E/N  is  greater  than 
that  of  ionization  and  the  combined  negative  ion  densities 
equals  that  of  the  electrons.  Pressure  sensitive  processes 
which  directly  or  indirectly  affect  the  electron  density 
also  affect  the  point  of  transition  between  the  stable  and 
unstable  forms.  The  electron  diffusion  and  drift  losses 
cause  a decrease  in  the  value  of  the  transition  E/N  as  the 
pressure  is  increased.  Also,  the  addition  to  the  discharge 
of  another  detaching  specie,  CO,  decreases  the  value  of  E/N 
at  which  the  transition  occurs. 

Detailed  Conclusions  from  the  Experiments.  The  oxygen 
discharge  was  observed  to  exist  in  two  forms,  the  high  field 
form  and  the  low  field  form.  The  transition  between  the 


forms  is  a function  of  pressure  and  current.  The  high  field 
form  generally  exists  at  high  currents  and  low  pressures. 

The  low  field  form  was  observed  at  low  currents  and  high 
pressures.  There  is  also  evidence  of  a high  field  form  that 
exists  at  very  low  currents  and  is  not  a function  of 


pressure* . 

The  transition  between  forms  exhibits  a hysteresis 
effect.  The  discharge  tends  to  persist  in  its  form  as  the 
current  is  swept  from  high-to-Iow  or  low-to-high  values. 

Near  the  point  of  transition  between  forms,  both  forms 
tend  to  simultaneously  exist  in  the  positive  column  with  the 
high  field  form  on  the  anode  side  of  the  interface  between 
forms  and  the  low  field  form  on  the  cathode  side.  The  inter- 
face between  forms  usually  occurs  at  a discontinuity  in  the 
discharge  tube  geometry,  such  as  an  electric  field  probe. 

The  impedances  of  the  high  and  low  field  forms  do  not 
change  as  a function  of  gas  flow  velocity*.  The  transition 
point  between  forms  does  change  however;  as  the  flow  velo- 
city is  increased,  the  transitions  occur  at  lower  currents; 
and  hysterisis  becomes  more  pronounced*. 

The  additive  of  an  inert  gas,  argon,  does  not  change 
the  impedances  of  the  high  and  low  field  forms*.  The  point 
of  transition  between  forms  is  only  slightly  affected  by  the 
addition  of  argon;  the  transition  current  increases  as  argon 
is  added*. 

The  addition  of  an  electron  detaching  gas,  carbon 
monoxide,  (CO) , to  the  discharge  has  marked  effects  on  the 
discharge  characteristics*.  The  addition  of  CO  increases 
the  low  field  impedance  to  the  point  where  the  difference 
between  forms  essentially  ceases  to  exist;  the  impedance  of 
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the  high  field  form  is  not  af fee  Led*.  As  CO  is  added  the 
transition  point  between  forms  decreases  to  lower  currents, 
and  the  hysteresis  effects  become  less  pronounced*. 

The  high  field  form  of  the  discharge  was  observed  to 
be  stable,  with  no  perceivable  oscillations  or  noiselike  be- 
havior. 

The  low  field  form  of  the  discharge  was  observed  to 
be  nonstable,  exhibiting  either  periodic  or  aperiodic,  noise- 
like oscillatory  behavior.  The  periodic  oscillations  exist 
only  for  certain  frequencies  or  modes*.  The  discharge  impe- 

i 

dance  is  a function  of  the  oscillation  frequency*.  The  op- 
tical output  from  the  periodic  low  field  form  is  highly  mo- 
dulated.  The  photomultiplier  output  indicates  very  short 
pulses  of  light  with  relatively  long  dark  periods  between 
pulses.  Electron  density  measurements  indicate  that  the 
electron  density  is  at  a maximum  during  the  dark  period  and 
at  a minimum  during  the  light  pulses*.  The  oscillation  fre- 
quency ranges  from  20KHz  to  90KHz.  The  oscillations  exhibit 
a small  amount  of  normal  dispersion  with  phase  velocities 
between  2.6  x 10^  and  3.5  x 10^  cm/sec  and  group  velocities 
between  2.5  x 10^  and  3.1  x 10^  cm/sec. 

The  temperature  of  the  discharge  is  represented  by  an 
analytic  function  of  pressure  and  power  input  to  the  dis- 
charge*. 


*Original  observation. 
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Mass  spectra  data  indicates  that  atomic  oxygen  concen- 
tration increases  suddenly  and  significantly  as  the  discharge 
transitions  from  the  low-to-high  field  forms  of  the  dis- 
charge*. In  both  forms  there  is  a lack  of  measurable  ozone. 

Detailed  Conclusions  from  the  Modeling.  The  isotropic 
portion  of  the  electron  energy  distribution  exhibits  a roll- 
off characteristic  intermediate  between  that  of  Maxwellian 
and  Druyvesteyn  distributions. 

For  the  range  of  E/N's  that  correspond  to  the  experi- 
ments, the  dominant  positive  ion  was  calculated  to  be  0^. 

At  high  values  of  E/N  the  0+  density  becomes  dominant*. 

The  dominant  negative  ion  was  calculated  to  be  0 . As 
the  E/N  is  lowered  toward  the  point  of  transition  from  the 
high-to-low  forms  of  the  discharge,  the  number  densities  of 
the  O2  and  0^  ions  become  significant  - especially  at  the 
higher  values  of  initial  molecular  density  (pressure)*. 

The  calculated  atomic  oxygen  density  indicates  a high 
level  of  dissociation,  especially  at  high  values  of  E/N*. 

The  ozone  density  was  calculated  to  be  several  orders 
of  magnitude  lower  than  the  other  neutral  species  carried 
in  Chemco* . 

•k  "l 

The  metastable  &)  density  was  calculated  to  reach 

a maximum  at  relatively  low  values  of  E/N  and  then  decrease 
as  the  E/N  is  increased*. 

Over  the  range  of  40-60  TD,  it  was  determined  that 
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two-step  ionization  through  the  (^(a  A)  state  is  an  impor- 
tant contributor  in  matching  calculated  to  experimental 
electron  densities*.  At  higher  values  of  E/N  the  two-step 
process  becomes  less  important  due  to  a decrease  in  the 
O2  (aVi)  density  and  due  to  the  lower  ionization  threshold 
of  the  metastable  becoming  less  important*. 

The  principal  electron  gain  processes  are  electron  de- 
tachment from  the  negative  ions*.  Both  O2 (a^A)  and  0 act 
as  the  detachers  with  the  predominate  process  involving  the 
metastable  molecule  detaching  the  electron  from  the  negative 
atomic  ion*. 

The  principal  electron  losses  are  by  attachment  pro- 
cesses, both  from  the  ground  state  and  the  metastable  mole- 
cules . 

Since  the  attachment/detachment  processes  only  trans- 
fer electrons  back  and  forth  from  bound  to  free  states,  any 
permanent  losses  are  balanced  by  ionization.  As  the  ioniza- 
tion decreases  with  decreasing  E/N,  these  other  loss  pro- 
cesses eventually  lead  to  a nonself-sustaining  discharge 
(transition  from  high-to-low  field  forms) . 

One  of  the  conditions  for  instability  - the  dominance 
of  the  derivative  of  the  attachment  coefficient  with  respect 
to  E/N  over  that  of  the  ionization  coefficient  - is  satis- 
fied for  values  of  E/N  below  67.4  TD*. 

The  other  condition  for  instability  - the  equivalence 
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of  the  negative  ion  and  electron  number  densities  - predicts 
reasonably  well  the  experimentally  determined  transitions 
from  the  high-to-low  field  forms*.  Calculations  were  made 
using  both  02(axA)  and  0 as  the  dominant  detacher.  The 

calculations  provide  the  better  match  to  experiments*. 

The  pressure  effect  of  the  transition  between  forms 
was  determined  to  be  due  to  the  electron  diffusion  and  drift 
loss  processes*.  As  the  pressure  is  increased,  these  pro- 
cesses have  the  tendency  of  increasing  the  electron  density 
and  thus  decreasing  the  value  of  E/N  at  which  the  transition 
takes  place*. 


♦Original  observation. 
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XI . Limitations  of  the  Study  ; 

and 

Recommendations  for  Future  Work 
Limitations  of  the  Study 

Throughout  this  report  the  limitations  associated  with 
the  experiments  and  with  the  model  have  been  specified.  It 
is  appropriate  at  this  point  to  summarize  those  limitations 
in  order  to  provide  an  introduction  into  the  recommendations 
for  further  work. 

Experiments.  The  major  limitation  of  the  experimental 
portion  of  the  study  was  the  difficulty  in  determining  the 
gas  density  exactly.  Temperatures  could  be  measured  only  at 
low  values  of  input  power.  The  gas  densities  at  high  input 
power  were  determined  by  an  extrapolation  process. 

Models . The  major  limitations  with  the  models  were: 

1.  The  electron  energy  distribution  code  was  run  only 
for  the  ground  state  molecular  oxygen  cross  sections.  It 
did  not  include  the  effects  of  other  atomic  and  molecular 
species  on  the  electron  energy  distribution  function.  The 
two  species  which  were  present  in  the  discharge  in  large 
numbers  were  the  metastable  oxygen  molecule,  (a^A)  > and 
the  oxygen  atom.  The  inclusion  of  these  in  the  Boltzmann 
Code  could  affect  the  transport  and  forward  generation  rates 
at  the  high  E/N  values  where  their  densities  were  relatively 
significant. 


For  the  E/N  values  at  which  the  model  was  compared  to 
the  experiments,  the  rates  calculated  by  EED  would  probably 
not  change  greatly  in  magnitude.  Since  the  Chemco  code  only 
modeled  steady  conditions,  the  results  would  not  change  sig- 
nificantly with  minor  changes  in  the  rates.  A time  varying 
code,  however,  could  be  affected  by  a minor  but  rapid  change 
in  the  rates. 

2.  The  rates  utilized  in  the  chemistry  model  were  for 
one  value  of  reduced  field.  Once  the  E/N  was  set  on  a par- 
ticular run,  it  could  not  vary  and  thus  model  the  oscilla- 
tory mode  of  the  discharge  or  take  into  account  external  cir- 
cuit effects.  Therefore,  only  the  steady  state  discharge 
was  fully  modeled. 

3.  The  chemistry  model  did  not  include  calculations 
of  the  radial  distribution  of  the  energies  and  densities  of 
the  different  species.  Radial  density  profiles  were  assumed 
in  order  to  estimate  diffusion  losses.  Calculations  of  the 
radial  distributions  and  energies  would  lead  to  a more  exact 
determination  of  diffusion  losses. 

Recommendations  for  Future  Work 

Experiments . Additional  experimentation  is  required 
to  fully  compare  the  oxygen  positive  column  with  mathemati- 
cal models  and  validate  the  models.  Specifically,  experi- 
ments should  be  repeated  using  a more  exact  method  of  deter- 
mining the  gas  density.  This  is  important  because  the  model 
and  experiment  are  linked  through  the  E/N  and  N parameters. 
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Also,  it  would  be  extremely  beneficial  to  be  able  to  deter- 
mine the  absolute  densities  of  the  negative  ions,  atomic 
oxygen,  and  the  O2  (a^)  metastable  molecules,  since  these 
species  are  important  in  the  discharge  in  the  region  of 
transition. 

Models.  The  chemistry  model  should  be  upgraded  to  in 
elude  the  calculation  of  the  radial  distributions  of  the 
species  and  their  energies.  The  model  should  also  allow 
for  the  species  concentrations  to  feedback  through  a solu- 
tion of  the  Boltzmann  equation  into  a determination  of  the 
electron  energy  distribution  function  and  thus  the  forward 
generation  rates.  This  should  allow  the  chemistry  code  to 
model  the  unstable  form  of  the  discharge  in  addition  to  the 
stable  form.  Each  of  these  recommendations  would  have  a 
very  large  impact  on  the  size  of  the  computer  program  and 
on  the  amount  of  computer  time  required  to  make  a run. 

If  one  is  interested  in  high-pressure  electron-beam 
sustained  discharges,  the  chemistry  code  would  have  to  be 
modified  to  model  them.  Two  changes  would  be  required: 
different  chemical  reactions  would  have  to  be  included,  and 
provisions  would  have  to  be  made  to  incorporate  a constant 
source  of  high  energy  electrons.  The  chemical  reactions 
chosen  for  Chemco  were  the  ones  that  were  dominant  for  the 
pressure  range  of  the  study.  For  higher  pressures  the  do- 
minant rates  would  have  to  be  re-established  using  the  list 
ing  of  Appendix  D. 
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Appendix  A 


Microwave  Cavity  Proqram 


The  program  listing  presented  in  Table  A - I was  the 
one  utilized  to  compute  the  electric  field  in  the  bore  and 
interaction  regions  of  the  microwave  cavity  described  in 
Chapter  VI.  The  solution  involved  the  utilization  of  a re- 
laxation technique  to  solve  Laplace's  equation  for  the  poten 
tial  field.  The  electric  field  was  then  determined  by  numer 
ically  taking  the  spatial  derivative  of  the  potential  field. 
The  program  was  developed  in  collaboration  with  Captain  M.R. 
Stamm  of  the  Plasma  Physics  Laboratory*  of  the  Aerospace  Re- 
search Laboratories. 

The  geometry  of  figure  A-l  was  utilized  in  setting  up 
the  program: 


///////////////A 

€r4(QUARTZ)  = 3.78 
€r3(CCU  ) = 2.17 

' 0-2  (QUARTZ)  = 3.78 


Gi  (VACUUM)  = 1 


OF  GAP 


OF  TUBE 


Figure  A-l.  Bore  Region 


*Now  the  High  Power  Branch  of  the  A.F.  Aero  Propulsion 
Laboratory. 


Cylindrical  symmetry  was  assumed  around  the  axis  of  the  tube, 
and  mirror  symmetry  was  assumed  on  either  side  of  the  gap. 
Therefore,  the  problem  was  reduced  from  three  to  two  dimen- 
sions in  order  to  simplify  its  solution.  The  mesh  was  con- 
structed of  71  divisions  in  the  radial  direction  (1.4  cm) 
and  101  divisions  in  the  z direction  (4  cm) . This  resulted 
in  a cell  size  of  approximately  .02  cm  x .04  cm.  The  boun- 
dary condition  on  the  pole  was  assumed  to  be  .5  volt,  corre- 
sponding to  a U of  1 volt.  The  boundary  along  the  z = 0 
axis  was  assumed  to  be  zero.  The  boundary  across  the  gap 
was  assumed  to  have  a linear  variation  between  0 and  .5  volt. 
The  boundary  at  z = 4 cm  was  assumed  to  be  zero.  The  z - 0 
boundary  was  adjusted  such  that  the  first  mesh  point  on  the 
negative-r  side  of  the  z axis  corresponded  in  value  to  the 
last  mesh  point  on  the  positive-r  side  of  the  axis,  for  all 
values  of  z.  Only  one  mesh  point  was  carried  on  the  nega- 
tive-r side  of  the  z axis.  Values  at  each  of  the  mesh 
points  were  calculated  by  successive  solutions  of  Laplace's 
equation.  The  calculation  was  terminated  when  the  maximum 
change  in  the  successive  values  at  any  mesh  point  was  less 
than  the  specified  minimum  of  5 x 10  volt. 
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Table  A 


Microwave  Cavity  Program  (continued) 


00  1 1=1,13 


A(I,l)=n. 


, ) =v 


A ( I * JO ) =V 


I E ( I . GE . I w ) N = jw-  1 


00  1 3=2,  A. 


A ( I , J) = .G«v 
00  2 J = Jw  , J‘jM 


A < 1 * , J) =V 
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ICOUNT  =0 
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"HE5"= 


IMG  rdh  POTENTIAL,  A,  AT  EACH  MESH  POINT 
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Appendix  B 
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Electron  Impact  Cross  Sections 

In  this  appendix  the  electron  impact  cross  sections 
utilized  in  the  Boltzmann  equation  program,  EED* , are  pre- 
sented. All  of  the  cross  sections,  with  the  exception  of  v 

rotational  excitation,  were  determined  from  experimental 
data.  The  rotational  excitation  cross  section  was  derived 

from  a theoretical  development.  Dr.  Arthur  Phelps'  consul-  , 

tation  was  extremely  helpful  in  choosing  the  set  of  electron 
impact  processes  to  include  in  the  Bolzmann  code  and  in  de- 
termining the  values  of  the  cross  sections  (Ref  45) . The 
information  Dr.  Phelps  provided  was  based  on  his  experiences 
in  running  Boltmann  codes  and  matching  experimental  transport 
parameters . 

Momentum  Transfer  Cross  Section  (Figure  B-l) 

The  momentum  transfer  cross  section  was  taken  from  Re- 
ference 56:72.  The  listing  was  for  the  data  without  the  4.5 
ev  resonance.  At  energies  below  1 electron  volt,  the  cross 
section  was  reduced  slightly  in  order  to  agree  with  the  mi- 
crowave data  of  Reference  57  and  the  drift  velocity  data  of 
Reference  45. 


*Refer  to  Appendix  C. 
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Rotational  Excitation 


— 16  9 

A rotational  cross  section  of  .15  x 10  cmz  was  used 
with  a threshold  of  .001  ev  (Ref  58:1674). 

Low  Energy  Vibrational  Excitation  Cross  Sections  (Figure  B-2) 
Low  energy  vibrational  excitation  is  accomplished  pri- 
marily by  resonance  scattering  from  the  vibrationally-ex- 
cited,  temporary,  negative-ion  states.  Direct  excitation  of 
the  vibrational  levels  of  the  ground-state,  neutral  molecule 
is  roughly  on  the  order  of  1%  of  the  excitation  through  the 
resonance  channels  (Ref  59:1619,  1621).  The  data  used  for 
the  low  energy  vibrational  excitation  cross  sections  were 
taken  from  Reference  59.  Since  the  cross  sections  were  so 
narrow,  their  magnitudes  were  adjusted  to  make  their  energy- 
integrated  values  match.  Once  this  was  accomplished,  the 
cross  sections  were  doubled  to  allow  the  low  energy  transport 
parameters  calculated  by  the  program  EED  to  match  experimen- 
tal values.  This  adjustment  was  recommended  by  Dr.  Phelps 
(Ref  48)  and  was  also  confirmed  by  adjusting  a cross-section 
multiplication  factor,  running  the  program  EED  and  then  com- 
paring the  drift  velocity  against  drift  tube  data. 

The  thresholds  for  the  vibrational  excitations  were 
taken  from  Reference  60:  V = 0 1,  1190  ev;  V = 0 2, 

.380  ev;  V = 0 -*■  3,  .570  ev;  and  "V  = 0 -*■  4,  .750  ev. 

High  Energy  Vibrational  Excitation  Cross  Sections  (Figure  B-3) 
The  low  energy  vibrational  excitation  cross  sections 
were  via  the  ground  state,  temporary  negative-ion  resonance. 
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Vibrational  excitation  of  the  ground  state  neutral  molecule 
has  been  observed  to  also  occur  via  electronically  excited, 
temporary,  negative-ion  resonances.  The  data  for  these  high 
energy  vibrational  cross  sections  were  provided  by  Dr.  Phelps 
(Ref  48) . He  derived  the  values  by  normalizing  the  differ- 
ential cross  sections  of  Reference  61:970  to  the  integrated 
values  provided  by  Reference  62. 

The  energy  losses  for  these  cross  sections  were  the 
same  as  used  for  the  corresponding  low  energy  vibrational 
excitations . 

Electronic  Metastable  Excitation  Cross  Sections  (Figure  B-4) 
The  cross  sections  for  excitation  of  the  ground  state 
molecule,  02(X^]Tg)»  to  the  two  low  metastable  levels, 

O2  (a^Ag)  and  O-^b^Xg)*  were  taken  from  Reference  62:1491. 

The  threshold  used  for  the  02(aJ&g)  excitation  was 
.980  ev;  for  the  O2 (bxXg)  excitation,  1.672  ev. 

Electronic  Excitation  Cross  Sections  (Figure  B-5) 

The  cross  sections  for  the  4.5,  6.0,  and  8.4  ev  thresh- 
old excitations  were  provided  by  Dr.  Phelps  (Rev  48). 

Dissociative  Attachment  (Figure  B-6) 

The  cross  section  for  dissociative  attachment*  from 
the  ground  level  state  was  taken  from  Reference  63:1486. 

In  the  absence  of  measured  data,  the  cross  section  for 

*e  + 02(x3[g)  0 + 0" 
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Figure  B-4.  Electronic  Metastable  Excitation  Cross  Sections 
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Figure  B-5.  Electronic  Excitation  Cross  Sections 
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Figure  B-6.  Dissociative  Attachment  Cross  Sections 
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ELECTRON  ENERGY  (ev) 


* 1 

dissociative  attachment  from  the  (^(a  Ag)  state  was  assumed 
to  be  the  same  shape  and  magnitude  as  the  process  from  the 
ground  state.  However,  the  energy  values  were  shifted  down 
by  .98  ev,  the  difference  in  energy  between  the  two  states. 
This  assumption  was  based  on  the  similiarity  of  the  poten- 
tial energy  curves  for  the  two  states.  In  order  to  avoid 
having  the  dissociative  attachment  cross  section  from  the 
a^Ag  state  influence  the  electron  energy  distribution  calcu- 
lation, the  cross  section  values  were  decreased  by  two  or- 
ders of  magnitude  when  they  were  entered  into  EED;  the  reac- 
tion rates  that  were  calculated  by  the  program  for  the  pro- 
cess were  then  increased  by  two  orders  of  magnitude  to  com- 
pensate. The  influence  of  the  excited  state  dissociative 
attachment  cross  section  on  the  distribution  function  was 
eliminated  due  to  a lack  of  knowledge  of  the  relative  popu- 
lation of  the  excited  level.  If  this  population  were  known, 
the  cross  section  would  have  been  included;  and  a more  exact 
calculation  of  the  distribution  function  would  have  been 
possible. 

The  threshold  for  the  ground  state  process  was  4.400 
ev;  for  the  excited  state  process,  it  was  3.420  ev. 

Molecular  Ionization  (Figure  B-7) 

• • 

The  molecular  ionization  cross  section  for  the  ground 
state  process  was  taken  from  Reference  64:1471-1472. 

The  cross  section  for  the  excited  state  process  was 
derived  and  handled  in  a manner  similar  to  that  discussed 
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for  the  dissociative  attachment  process. 

The  threshold  for  the  ground  state  ionization  was 
12.200  ev;  for  the  excited  state  process,  it  was  11.220  ev. 

Dissociative  Ionization  (Figure  B-8) 

The  dissociative  ionization  cross  section  was  taken 
from  Reference  65. 

The  threshold  energy  for  the  process  was  18.5  ev. 

Atomic  Ionization  (Figure  B-8) 

The  atomic  ionization  cross  section  was  taken  from  Re- 
ference 66.  The  cross  section  values  utilized  in  the  EED 
program  were  decreased  by  two  orders  of  magnitude  to  avoid 
influencing  the  distributions.  The  ionization  rates  that 
came  out  of  the  code  runs  were  increased  by  two  orders  of 
magnitude  to  compensate.  The  reason  for  the  de-emphasis  of 
the  cross  section  is  the  same  as  for  the  dissociative  attach- 
ment  from  the  02  (a  Ag)  state. 

The  threshold  energy  for  the  atomic  ionization  process 
was  13.618  ev. 
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Boltzmann  Equation  Program,  EED 

The  program  listing  presented  in  Table  C - I was  the 
one  utilized  to  calculate  the  electron  energy  distributions, 
transport  properties,  and  reaction  rates  described  in  Chap- 
ter VIII.  The  program  was  developed  by  William  Bailey  and 
William  Long  for  use  by  the  personnel  of  the  Plasma  Physics 
Laboratory*  of  the  Aerospace  Research  Laboratories.  The 
only  portion  of  the  program  provided  by  the  author  was  the 
consistent  set  of  cross  sections  described  in  Appendix  B. 
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♦Now  the  High  Power  Branch  of  the  A. F.  Aero  Propulsion 
Laboratory. 
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Table  C 


Boltzmann  Equation  Program,  EED 


ROGHrM  FED  (INPUT. OUTPUT, PUNCH»DATA,TAPE5=DATA, PLOT) 


DIMENSION  XR ( 30 ) : 


MMUn/CRCSS/  A (5) , IMX (5) , IMG  


COMKON/fUN/  ETA  (200,30)  ,QR  (200)  ,DF  (200)  ,EOVRN,XMR, ID (30) 


EP(30)  ,0(200)  , I TP. IOV.UU,TG.TV,F,XNU(30) ,XNV(30) ,IVO,IVI,AN 


COMMON/SUN/  bTA(20Q) , DOE (200) ,DEL(200)  ,RT<50) * I T V ( 30 ) 


C0MM0N/5CALER/LX, LY,XP(5) .NAME (5)  " 


DATA  IPPT, IPN/1 ,0/ 


DATA  F »TV,TG/0.  ,0.  ,293,/ 


m 


»«»»«»«» 


***  NUMERICAL  SOLUTION  OF  BOLTZMANN  EQUATION  **<>< 


EED  NUMERICALLY  SOLVES  THE  COLLISIONAL  B0LT7MANN  EQUATION  IN  THE 


M APPROX  I HAT  TON  FOR  THE  ELECTRON  ENERGY  DISTRIBUTION  FUNCTION. 


M THIS  IS  CALCULATED  THE  MEAN  ENERGY,  DRIFT  VELOCITY,  AND 


Y LOSS  RATES. 


READ  IN  BASIC  PARAMETERS 


READ*,  NEN.IPRT.IPN 
NUMBER  OF  E/N,S  TO  BE  CALCULATED 


>RlNT  OUT  FO  AND  FI  IPRT=1 


I PRT=0  NOT  REQUESTED 


liiiKMiMwsmwflaniraaH 


IPN=0  NO  PUNCH 


CCgW.L-]**. 


m Trn.(u»MuilTim^«i:mi:am 

u.Mi.iHMiwuMin.iwmmmm— 
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Table  C - I Boltzmann  Equation  Program,  EED  (dontinued) 


TE  OUT  PUMPING  KATES  AND  ENERGY  LOSSES  FOR  ADDED  GASES  t 


EIN=VD«E0VWN*1 .E-16 


PRINT  490,  OLOSS/EfN  1 


490  FORMAT <1H0»  * FRACTIONAL  RECOIL  ENERGY  LOSS* t F 1 2. 4) 


EE=0. 


IBT=1 


DO  04*  irv=l,lNG 


ITP=1*X (ITM)-1+I0T 


DO  84b  J=IPT,ITP 


XR(J)=0. 


IF(IDU)  .Gt. MOM-1)  GO  TO  046 


EH(J)/UU-FLOaT ( 10 ( J ) ) 


ot j)  *i » j)« <xnu( j)*dof  « i n { j)  *d-xnv( j)*i-0F  ( i ) )*n .-p/2.) 


I'M  J)  ♦ 


■Mf  • 

ESI 


J)=XR(J)*ETA(I,J)*(XNU(J)*OOF(I)-XNV(J)*(P*DOF(I-IO(J)-1)*(1.-  i 


)*00MI-I0<  J)  ) ) ) 


45  CONTI NU 


46  XRU)=XRU>*S.«?9E-9/<P<ITM)*nU/XNERM  i 


05 . NAME ( I TM)  : 


, * NET  PUMPING  RATES  CM3/SEC  FOR  *,A10> 


XR(I) *I  = JffT«ITP> 


5T.ITP 


X«(J7*0£P(U)*XP<ITM)  i 


Ml 

mSM 

PVUMaElllE 


i,a*fjiiLfuujr.iii  i«bw  h 1 1 via  m 

■tawiiifHwiiBMniig— ■■■ 


h !'■  « tn  ir«  »i i iui^t  wmiamflua 

i<  * lar-a  m urt 


LWiKmH 


;=MMiumiiHaui 


Table  C - I Boltzmann  Equation  Program,  EED  (continued) 


wars®' 


Table  C - I Boltzmann  Equation  Programm,  EED  (continued) 


1 I0V=I-1 


IF(IOV.LE.l)  ST0P2 


C0UM.t0.3n)  60  TO  ] 6 


tTM.LT ..01)  60  TO  16 


IF  <SGl\*SGNO.t_T.O.)  FEM=FE“*.S 


5GN0=SG\ 


OEFL=ALOG10(OEFM*1.) 


IF (OtTL.GT.FtM)  S=0EFL/FEM 


00  12  1 = 1.10)/ 


1 2 OOF  ( I ) = 0()F  ( I > **  < 1 . - 1 . /S ) *0F  < I > **  < 1 . /S) 


SL0Pt=00F ( IUV) /OOF ( IOV-1 ) 


00  1<*  1 = I0V»  190 


OOF  ( 1+ 1 ) =00F ( I ) *SL0PE 


IOV=200 


GO  TU  2 


6 HETUKN 


END 


SUBROUTINE  VI 6 


COMMO.M/FUN/  ETA (200,30) ,OP(200) .OF  <?00>  .EOVRN , X*R  . 10  ( 30  ) 


$ UEP(30)  .0(200)  , TTP,I0V,DU.TC-,TV,F,XNU<30)  ,XNV<  30)  , I VO,  I VI.  am 


E 1 =DtP ( I VO ) 


E2=AN*E1/(1.-2.*0E) 


VMIN=b0. 


IF  (AN. ME. n.)  VMI.M=E1*TG/  (2.*E2MV)  *.5 


M=0. 


00  10  1=1.30 


V=FLOAT(I-l) 


V = V* (El - (V-l . H»E?) 


IF (V.GT.VMN)  GO  TO  5 


XV=EXH-V«E1/TV*V«  (V-l.)  *E2/TG) 


Appendix  D 


Chemistry  Rates 

Table  D - I is  a listing  of  the  chemical  and  physical 
reaction  rates  collected  during  the  course  of  this  study. 

The  last  column  in  the  table  specifies  the  reactions  utilized 
in  Chemco  (See  Chapter  IX) . Nonutilized  reactions  were  re- 
jected because  they  did  not  involve  species  carried  in  the 
code,  were  slow  in  an  absolute  sense,  were  slow  relative  to 
another  channel  with  the  same  net  reaction,  or  were  deter- 
mined by  early  Chemco  runs  to  be  insignificant.  The  letter 
M in  the  reaction  column  refers  to  an  unspecified  specie  in 
a three-body  collisional  process. 

The  rate  dimensions  for  single  body  processes  are 
sec-''-;  for  two-body  processes,  cm^sec-^;  and  for  three-body 
processes,  cm^sec-^. 
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Table  D - I.  Reaction  Rates 
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Table  D - I.  Reaction  Rates  (Continued) 


Rates  (Continued 


o 

O 

o 

O 

o 

O 

O 

o 

in  c 

G 

c 

o 

O 

u 

u 

u 

U 

U 

u 

— O 

O 

o 

CO 

KD 

> — *H 

*H 

•rl 

N 

N 

N 

N 

N 

N 

N 

N 

rH 

rH 

'T  >i  -P 

■P 

-P 

■P 

-P 

■P 

-p 

-p 

■P 

-P 

•P 

r 

1 

1 T5  O 

V 

o 

rH 

rH 

1— 1 

H 

rH 

rH 

1 1 

r — 1 

o 

O 

in 

ro 

O O <1> 

0) 

0) 

o 

o 

o 

o 

o 

o 

o 

o 

rH 

rH 

i 

l 

H ffl  w 

in 

tn 

m 

CO 

m 

cn 

tn 

tn 

tn 

tn 

o 

o 

X 

X 

rH 

rH 

X <D  W 

to 

tn 

6 

e 

g 

g 

g 

g 

g 

g 

rH  tn 

tn 

tn 

o 

O 

o 

o 

o 

o 

0 

o 

CN 

in 

X 

X 

h tro 

O 

O 

Vi 

VI 

VI 

VI 

Vi 

n 

VI 

p 

• 

• 

• p 

p 

p 

tn 

X 

tn 

tn 

tn 

In 

tv. 

tv. 

CN 

rH 

CN 

CN 

rH  CJ 

a 

O 

- — » 

< — *. 

^ . 

, * 

• 

• 

• 

• 

• 

• 

• 

• 

tn 

V) 

cn 

U) 

tn 

W 

tn 

V) 

0) 

Q> 

QJ 

o 

<i) 

0) 

tu 

0) 

X 

X 

X 

x 

X 

x 

x 

X 

X 

X 

.g 

s 

> 

S 

O’ 

O’ 

tr> 

O’ 

o 

o 

o 

0 

•H 

•H 

•H 

*H 

PI 

PI 

.-1 

pi 

X 

X 

K 

Sw^ 

" 

CN 

* 

s. 

o 

(1) 

<D 

aj 

a) 

a) 

<D 

<u 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

CN 

—1 

CN 

cn 

CN 

CO 

O 

0) 

tn 

Vi 

a) 

X! 

cm  O 


O’  + O’  O’  + O’ 

W < ^ 


(Charge  Ex. 


AD-A066  196  AIR  FORCE  INST  OF  TECH  MR I OHT -PATTERSON  AFB  OHIO  SCH— ETC  F/6  7/2 
DISCHARGE  PROCESSES  IN  THE  0XT6EN  PLASMA. (U) 

OEC  ?•  J M OETTMER 

mi/os/mm-i 


UNCLASSIFIED 


NL 


Table  D - I.  Reaction  Rates  (Continued) 


E <D 
<d  c <U  'O 
w -h  x:  o 
d o o 


1 

X 

X 

X 

X 

CO 

on 

rH 

CN 

eg 

in 

in 

in 

rH 

rH 

in 

VO 

r- 

in 

rH 

CN 

CN 

00 

00 

00 

rH 

rH 

00 

•*3* 

in 

in 

VO 

VO 

• • 

•• 

•• 

•• 

• • 

• • 

• • 

p 

•• 

•• 

•• 

• • 

•• 

•• 

VO 

on 

ON 

ON 

rH 

rH 

ON 

w 

ON 

ON 

ON 

ON 

ON 

ON 

VO 

VO 

VO 

CN 

CN 

VO 

w 

VO 

VO 

VO 

VO 

VO 

VO 

<u 

X) 

c 

0 

o 

u 

•H 

rH 

o 

rH 

■p 

rH 

rH 

rH 

N 

rH 

o 

| 

00 

i 

CO 

O 

i 

-p 

i 

o 

O 

ON 

ON 

ON 

a> 

O 

• — ! 

O 

rH 

rH 

o 

rH 

O 

rH 

CN 

rH 

rH 

w 

rH 

1 

rH 

| 

i 

rH 

0 

rH 

1 

| 

i 

i 

i 

o 

O 

o 

m 

o 

O 

o 

O 

o 

w 

X 

rH 

X 

rH 

rH 

X 

X 

rH 

rH 

rH 

rH 

rH 

in 

E 

O 

m 

X 

in 

X 

X 

in 

o 

ON 

X 

X 

X 

X 

X 

P 

< 

• 

• 

• 

p 

• 

U 

z 

r- 

CN 

CN 

in 

00 

CN 

In 

CN 

co 

CO 

rH 

rH 

rH 

CN 

O 

+ 

Q 

+ Cr> 

rH 

p 

Kl 

' 

co 

rH 

CN 

O 

ja 

O 

o 

*• — 

o 

+ 

CN 

+ 

+ 

O 

ai 

+ 

CN 

P 

O 

+ 

rH 

+ o> 

■* 

CN 

e 

< 

a) 

Us/ 

1 

O 

+ 

0) 

* 

E 

<u 

O 

rH 

+ 

id 

Q) 

+ 

X) 

+ 

o 

+ 

+ 

+ 

tr> 

CN 

CO 

+ CM 

+ 

4 

CM 

< 

CN 

o 

O 

o 

o 

+ ro 

CM 

o 

f 

O 

rH 

o 

4 

o 

o 

o 

w 

(d 

T 

+ 

+ 

+ 

LD 

' — ’ 

A 

P< 

4 

CM 

4 

4 

CO 

+ 

+ 

2 

+ CM 

T 

A 

A 

A 

co 

T 

O 

T 

r 

T 

o 

o 

o 

o 

T 

T 

O 

CN 

o 

IN 

+ 

m 

CN 

CO 

+ 

O 

s 

+ 

+ CM 

O 

o 

O 

o 

f 

4 

A 

O 

o 

A 

+ 

'O 

T 

T 

T 

+ 

+ 

1 

+ 

tn 

+ 

+ 

+ 

+ tM 

0) 

+ 

+ 

w 

w 

CN 

E 

o 

O 

o 

c/3 

rH^* 

rH 

D 

O 

a 

a) 

a) 

rH 

a 

<d 

a 

C/3 

C/3 

X) 

+ 

+ 

CN 

+ 

+ 

rH 

rH 

rH 

rH 

-P 

+ 

+ 

+ 

CM 

— 

s-* 

CM 

0 

+ 

4 

+ 

4 

+ 

+ CM 

<u 

O 

O 

o 

O 

O 

O 

o 

z 

a) 

o 

o 

O 

o 

ai 

4) 

CM 

O 

rH 

CM 

CO 

•** 

in 

VO 

r- 

00 

ON 

o 

rH 

CM 

CO 

r* 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

ON 

ON 

ON 

ON 

ON 

Table  D - I.  Reaction  Rates  (Continued) 
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Table  D - I.  Reaction  Rates  (Continued) 
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Notes : 


*References  67,  68,  69,  and  75  are  compilations  of 
rates  and  as  such  do  not  represent  primary  sources. 

1 - Rate  entered  into  Chemco  as  a function  of  E/N. 

2 - Computed  in  EED  but  not  entered  in  Chemco. 

3 - Rotational  excitation. 

4 - Momentum  transfer. 

5 - Electronic  excitation. 

6 - Low  energy  vibrational  excitation. 

7 - High  energy  vibrational  excitation. 

8 - Cross  section  not  used  in  EED  or  Chemco. 

9 - No  rate  found. 

10  - See  Appendix  B for  cross  section  modification. 

11  - Rate  not  available.  Reaction  assumed,  so  as  to 
allow  reactions  involving  O2  (a^Ag)  in  Chemco. 
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Appendix  E 


Chemistry  Code,  Chemco 

The  program  listing  presented  in  Table  E - I was  the 
one  utilized  to  calculate  the  species  concentrations  and  re- 
lative reaction  contributions  described  in  Chapter  IX.  The 
basic  program  was  developed  by  Michael  Stamm  for  use  by  the 
personnel  of  the  Plasma  Physics  Laboratory*  of  the  Aerospace 
Research  Laboratories.  The  portions  of  the  program  provided 
by  the  author  were:  the  initial  conditions;  chemical  rates 
and  reactions;  chemical  difference  equations;  diffusion, 
flow  and  drift  loss  equations;  and  the  input/output  instruc- 
tions. 


*Now  the  High  Power  Branch  of  the  A.F.  Aero  Propulsion 
Laboratory. 
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Table  E - I Chemistry  Program,  CHEMCO 


rnfMCo 


rf'i'-‘Q*'/<n»  /o,p 


oppm  .hippot 


OM/SCA 


nru-MMP*:  kthmo 


rO,-i'-Q,t/t)ftTKr«f’/C|1  ,t»m.t»?.PP4  .PI  ,P 


*10.918  . : 1 8.  pi  r .919,921  ,9??. 9 


r<  3 ] .•->•»?.  p n o . pi  o , 988  , W3<.  . ror.  . F<  3*  . P'4  ] , 93*  . f,iq  , W4  3 . 04 


ps?.i  y>.M7 1 .p.op.rbb.p 


PI  12.9107, p lnfSKj  1 7,91 1B.9123.W1 ?4,P1?6,PI?7.91?8.P1?9 


s p]  3f-.Pi'ii  .n  'p.i-i  33.pi34.pio8.9i  ^,0] 37.FnB.unq.pu 


6 9Y.P/ 


5, 87. '•o. 89.811  ,81 


1 $19.8^1  .82?.S?3,8?4,S24,S?«.S?9.830,S31 .S32.S33.S34 


54S , »S l .S52.SS4.S71 .S82.S55 


3 S-rB.Sl  r,n  ,s;  04.81  OS,  8 107,  SI  08, 8 


813'*.  81  38.8) 3a, 81  39.814n.SX.  SI  41  .86.818,841  , SS3  , S'93  , S9/ 


5 r»l  14. 81  19.87,8/  


'>■'  <rHiP.,)Q?p.npN).'"io?N].po3N'.nF'.n 


l'V>: P- ' , '/OP  , VP7P  . VON . van*.  , VQ3M.VP  . V03 . f XC I T r . VF  , V0?D , VF  . i L F a 


!»  CH/  1 H 1 

. IP?, 

1 H 3 , 1H4 

H8,  1H7 

II-.'ITTAL  CO’IOT  T 1 0N8 


OFLT  = l .Fl'l 


CYC  *,  ' = 


rortCl 


Ft  "CTt’O 


PPFP 


soion 


2 


8(1 0 0 


'•  .vF.XOTCH 


vjTn»  » 8Ki  » ft.F-m 


08  r ).F-P*5M 


1 .F-8*8I 


0?F  = |.K.7*5‘i 


03  r 1. 


O'V  = 1 .(■  - l n»8'l 


,S*(r)*C'>rP)-0?0-02\-3?F-I  .5*  C 03*03*>l> 


94 


o^ST  ®KL vsn  mcuCAtt** 

THIS  PA®  K BEST  Q ^ 

^0*  OCtf*  WBWaW®  W 

Table  E - I Chemistry  Program,  CHEMCO  (continued) 


KB^i«kiui^iliH^ikiiilj8!nSIB3!HTBSHin?fl 


01  X = PS? 


PIM  = PK« 


F ( fouFra 


= of  ♦o?p-r.M-o?t',-03N 


<F»* 


1 r ( '->-  0D  (PPPP.5 


•■AT  ( ' t 1 ] ) 


PI'i  >T  S?S.  VP  .FOvFF>M.  Vi.XDTCH 


1 <6.3> 


“Si51I4.FECMAP,yFM,R]?^.R) 39 


V!  =*F10.3*  OTOHRs* 


ICO  CONTINUE 


IR..VG  = < 


r,M.'.  •’JrTf'F'i  (of IT  .FPACT«FF,ACT0n»CHFMDT,CHEP,T,CHfe'MCYC,rHFPKAX  ,KK) 


T I F ♦ CHF'i  r 


CYU.F  '.  = C YCl  F F ♦ 1 . 


IF  ( C VC!.F  S . T .OYC'^AX)  ISTOP 


IF  ( 1 i 'F.ct.T  V±*t ) I STOP=  1 


XU) 


X<?> 


X (.1) 


»i'.|  = C3 


jBO*  oOPi  FUKH1SHKD  TO  MC  - 


Table  E - I Chemistry  Program,  CHEMCO  (continued) 


—i  lu  u i>3n-gm-si 


BOC  -g.^'UV-SI  B»Sia<-?.oS\Q*S??*S?3*S7fr*?.»S?9* 


sot.rqi  ♦s?«^1.<s7«?ih*<;76*?.»sg^*s?p«-: 


itmumiM 


JfT^lTM^A-kkHAmAL-gir 


S0^0=s?*s<:-*s<  ♦«;]  1 ♦ S73*S?f>»S?3«S?9»<;Vi*?.*';3<.*S3q*S3B*S3fl»g.»<;fr3* 


1 ♦?<oSB?*S9a»S)  n()*Sl  0^*9 1 05*S  1 07*91  OB»?.osi  1 B* 


gg.o'ill  0 ♦?.<*•  I H^^ftOaP^OOgN^OOgr^VOPO^VOZP^VOgN*  .5»ytl*  .5«VQP* 


9Q7t=Sg«S'.*S1S*r.?1  ♦?.»S32*S33*S39*S**l*S5?*S53*SBB*S 


■=Ci‘*Sl  ?*91f,»S?4»S37*S33*Sn7*S)  34+003N* VOIN 


9*V03 


ECXXU1I 


SO.ll.sSB  + Sl  3*S\B  + S3<.*SA5*SS4*SS5*S 


1 ♦ 


.1  «9?*S  3»S/+*916*S36*SS1*S*B->S1  ] 7*9nS*Sl?3*Sl?6«S140*VOtv!*VQ?N« 


1 V03r 


IHIS  PACE  IS  BEST  QUALITY  PRACTICABU 
KBOJI  OOPY  PURUISHED  TO  DDQ  — — 


Table  E - I Chemistry  , rogram,  CHEMCO  (continued) 


SEl  = SS  + 96*S7*Sl«*S19*s;?l*3??*S3f!»S41*S53*S93*S94* 


!S)A?*nF*vnp*vf)?p  


^kLZUkkmUl 


SOUI.  =S3  + 3^*S<US9*S1  1+S1?*S13*S15*SS1.  *S  \ 1 1 7 *S1  ?A  *L)ON*  VON 


bO?l'f'i  = Sl  9*S?l+S3S  + S41*SlU  + Sl?^*S13n*Sl4l 


S0^:i=SlH*l.P3*<;?F,*SP9*S36*S93*S45*Sl  1R*00PN*VQ?N 


S03Nn=Sl  3*S15*S45 


S03M_=S?4*S3S*9I ?3*C03N*V03N 


SOK  =9«*91 1 *S??*S23*S?4*330*S39*S93*S1 36*00P*VQP 


1 ♦ 1. 


SP?P('-  = c?  + S30  + s3O*si  ?S*S136 


Sn2PL  = S“*Sli:*S1  9«St!9t538*S43*S94*D0?P*V0?P 


1 ♦ 1. 


S02P&=S104*S1  \ t 


1 ♦ 1. 


SO?f.L=rWl  ♦?.o^PP*S1  00*S10S*S107*S\0  3*sn7*sna*S119*S1  ?4*S1P6  + 


1 SI  J-) ♦ s 1 *S1  1 34  + si  39*S140*S141  ♦r)0?0*V0?D 


? ♦ 1. 


•7S4  = P.oSft/S(.!F, 


=‘.9/S. 


/soo 


P./S'-G 


osi9=e.osi9/sor. 


/S'lo 


0S?3  = ^ 


ms  MCE  IS  BEST  woin  mcnoBM 
yuyg  oupi  s ufuusHjau  10  w>.Q 


THIS  PACE  IS  BEST  QUALITY  PRACTICABLE 
FMOJI  COPY  maUSHED  TO  DDC  ^ 


Table  E - I Chemistry  Program,  CHEMCO  (continued) 


Fsn  7 = 


ES1P7  = 


ES  \?f>  = 


Ff.UC  = 


EV  f;W 


Evr>3f-  = 


ESS  = S 


E 


ES130 


CIVS'- 


O'  SI  3 


SS1/SFG 


SflP/SEO 


= S11 7/SFG 


SI IP/SEG 


= S123/SEG 


= S 1 ?^/?FG 


= S)'*0/SFG 


VO'./SEG 


vn?'l/Er 


= vn.-n.i/SEG 


SS/SEl. 


So/SFI. 


S7/SFL 


= S’S/StL 


S 1 c /Sr  I 


= S?]/St'L 


= 5??/SEL 


S3F/SEL 


= S41/SEL 


= SS3/SEL 


SS7/SEL 


= S9-/SFL 


= S13G/SEL 


= S)  3-V/SFL 


= S1M/SEL 


Of /SE  L 


GC°/SFL 


= vc.etvs^l 


SS/SG'JG 


= SG/SOI-T, 


= S7 /SOf G 


= S?r>/SOr-;G 


- SS3/SOVC 


v = S 1 3-i/SO-iG 


/50Nt. 


= SP/SOmL 


C 
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itHTS  PACE  IS  BEST  QUJ&1TT  JRACIIC^  - 


Table  E - I Chemistry  Program,  CHEMCO  (continued) 


o?Pn?r> 


0?pv?p 


02. »]  re. 


0201  1 2 


0 2 OS  7 I 


0205P 


0?L>1  Cf> 


oprios 


020107 


02.0)  OK 


opoi  is 


P2D1 ] 9 


02G129 


020 1 2K 


02(1  no 


0201 39 


020136 


0201. 39 


02!M  3^ 


020 1 90 


no?p/502P 


= W0.?P/9O2P 


C 1 09/90200, 


2/SO? 0 


= S71/SO20I 


2.*Sn2/S0?0 


= S 1 Q0/SQ2DI. 


= Si  03/9020 


1 97/5020 


= S19P/SO20L 


= e 1 ) //S020L 


=S1 1S/SO20L 


= S119/SP20L 


: S129/S02PL 


SI 26/S020L 


= S 1 30/S020L 


= S139/S020L 


5139/5020L 


: 5 1 .3:!/5O20l. 


= SJ39/S02DL 


: <-190/0020 


020)91 

= S191/SC20I. 

020020 

= 002D/S02DL 

G20V2O 

r VO2D/S020L 

c 

PPI'IT  Ml)  1 * 0, SOP,, SOL 


ro'i  o-n  nno.^o  -oiFiF.i n.?«  gain  =*ipifio.p*  loss  =<mdifio.2) 


■<I  1 M.12 


/,6X,*SM  S9  S13  SIS  S19 


29  S29  S39  S51  S52 


16.oSM,;>V).uS13,0SlR.0S19,0S22.0S23,0S29,0S29,0S39,0S51* 


S5  I 


/,SX,*SS9  SS5  SMS  S93 


S 1 29  S 1 3fc  SI  3fl  S139 


) 


559,0s 


9»*  tOlJuA  ,000 


(.9 


WWiH.T;n,R>W.ai 


wa 


1 9, OS 1 1 9, OS  1 29 , OS  1 36, OS 13H ,05 1 39* 


itraitETZimw 


/,6X,<* 


aa 


S39 


0S33.OS39, 


, <*S  1 09 


ralSFA®ISBK^^m^CIlCABU 

yao*  oopy  yuRiiisHED  to  dgc  ^ 


Table  E - I Chemistry  Program,  CHEMCO  (continued) 


pm 


013 


H8BI1B 


PPL-.T  RlO.F.SFG.SF 


FiMwftT  (1W0>«F  =*1P1F1« 


PPIM  p 


FOP'-  AT  ( 


S3 


S] AO*) 


* 1 P 1 F 1 0 


* 1 P1F 1 0 


1 126.FSR0 


F 1 2 FP'<vat  (/,ax.*VON 

the 

N 

PI.  I .T  <5H..!:V>'.FV0?M.F  V0.3N 

ppi  ■ r o i 7 

FPP  >AT  (/,AX.» 


Si  SO*) 


VCM-  AT  *S1  39 


PSTM  9 1 ft.  Kill  39. FS 141  .FOF.EVOP 


PPL  T «tS,0;  .SOK'G.SON 


Fr>- i'A.T  MHO.*OM  =*1P1F.  10.2* 


VO?P* 


Ft'”  v A T (/,*,X.*SS 


1 <)*> 


P-Mf-T  'Jlft.O'SS.O^SA.OMS/.ONSPe.ONSSS.ONSlSP 


ppi<  r si7 


FCPf- AT  ( /,a,t.*S3 


s 


? r>0  ;*> 


PP  I T Slf-.OOS.S.O'IS'+.ONSO.ONSP. 


,C  SI 1 7 .A'SI PA.ONnOM  


pur  t «ih 


rpp «AT  t/.ftX.»VOP») 


pp  1 ' T 91ft.  OMwOM  


EI011R2ETFW17OTHJI301 


PQIM  fiiq.oz^.soprCuSOgN 


FQP.-:AT  (1HU»Q?H  =»1P1F10.?»  GAIN  =»1P1E10.2» 


PPI'-T  P.20 


FO'-M.Ar(/,T,X.«Sin  S?1 


S 1 4 1 « ) 


PPI*'T  9H.,Q2MS18,02NS21  ,9?NS35,02NS41  .Q?N1  I 4 ,0?N1  24  ,Q?f  ! 


°PIVT  82 1 


FOP  SAT  ( /.6X,»S1S F23 S26  S29  S36 


S SUfl  no?M 


PRINT  91  A.U2r  S16»02N.i523,02NS26,02NS29,02N636.02NS43.0?NS45.0?N 


10?Wr2N»02NV/2N 


PPI*  1 822*03r..S03NG,5>93NL 


TV1 


AT  (1H0,*03" 

F.  10.2*  GAIN 

= 

* 1 P 1 E 

format  (/,6X,»S24 S35 S123  D03N  V03N*) 


PP  T »•  T «?3 


FOP’-.-AT(/,f.>U*S13  SIS  S45*) 


P PI  <1  916.03'"S13,03NS15*03NS45 


PPI  -'T  R24 


PP1  r.T  916,03!!S?4,03NS35,03N123,03M33N,03NV3K' 


PW1:.1  8?S.OP,SOPG,SOPt 


o)P1E10,?»  LOSS  =*1P1F10.2> 


PKT'T  826 


MTH'AT l/,4X,*Sl  S88  Sl40*> 


PS1 ,0PS88.nPS 


PPI  'T  827 


**S8  Sll  S22 


S 1 36  OOP  VOP*) 


V i 6 » OhSH  »0PS1 l ,OPS22.C)PS23,OPS24,0PS30,OPS39,0PS93,OPS136.OP 


OOP • OPVOP 


Pi'T  n?n »02P*  S0?PG*S02PL 


FOP.'  Al  (1HO,*0?P  =»1P1E10.2*  GAIN  =o}P]F.)0.2*  LOSS  =*1P1E10.2) 


POINT  829 


IT l/,GX,*S2  $30  S 39  S126  S136*) 


916,02PS2,02PS30*02PS39,02P126*02P136 


63  C 


T(/.fX,*S9  SI?  S 1 9 S29  S38  S43 


VC?P*1 


02P02P,0 


ramaaarairaimraiap^^ 


wiarjiHiffiiiMii 

Bjiamm1 

KiSmww; 


■Mil  m«fTnra?inir^Kng^^ 


Si  26  S 


, 021)  1 00, 0201  OS,  02D 107, 020 108, 0201 17,020118 


IHIS  PAGE  IS  BEST  QUALITY  PRACTICABLE 

raOBOOfYFURUISHEDlODOC  ' 


Table  E - I Chemistry  Program,  CHEMCO  (continued) 


OT  = AMINM  (DDFLCHtK)  ,DELCH(K)  ,DT) 


CONflMUt 


'JPDf  If 


no  220  K= I »KK 


OCMtK) 


OC^nOT (M's  CMOOT (K) 


OChtK)  = CHDOT (K ) *DT 


CH(K)  = rH(K>  ♦ 


IF tCH(K)  .LF  . . 1 ) CH ( K 


CO.'.T  I Nl'E 


Tl^R  = TIKE  ♦ OT 


CHE  MT  = T I Mt 


CHLVDT 


CHt'VCVC  = CHEMCYC  ♦ 1 


-liffTTHa 


IF ( ( 1 . ♦ 1 »f-10)*TIKE.LT .OELT)  GO  TO  200 


Pr TUP" 


ENi 


SUBROUTINE  OEPIV  (X.  XPOT) 


comv:0 n/n/'SRrt/  e.h 


C OFRI7  TAKES  THE  VARIABLES  AND  RETURNS  THEIR  DERIVATIVES. 


C 


0 1 •■FAST Of-  X ( 1 0 ) , XUOT  (10) 


A/HATkrK/Pl ,P93,P2,R94,P3.P53.R4.R6,R5,R7,Rf».R9.Rl 1 ,R12 


R1?.R1S.P16»R!G»R19»W21 , R22.R23.R24. R25.R 1 14.R28.R29.R30 


^31,10?.  VI  19,t->33.H55.R34,K35.P:36.K41  .R36.R39.RA3.R45. 


51  ,>5E.P5i.,P71  .R62.R66.R100.R10A.R105. 


.31  r,  f.  9I0R.FM  1 7.01 10.Rl?3.R124.RT26.lU27.f:l2fi.ol?q. 


•■l  30.  M 31  ,R1  32.D133,W134,Rl35,H13f,,fil37,R13B,O139,BlA0.PX.!<141  , 


waraEftsanai 


S'^'/Sl ,S?.S3.S4,S5,S7,S6.Sq,Sll»S12.S13.S15.S16. 


22.SR3.324.S26.S20.S29.530.S3I.S32.S 


S32.S5 


‘.>R*,S1  Oflol  0‘..S105.Sloy,S10P.Sl  1 7,  SUB.  SI  23,  SI  24.  SI  26. S 


36. SI 36,S139.Sl40,SX .5141 *Sf>*Sl  B.S4 1 , S53 , S93  »S9A ,S1 12, 


WIHMiaiMWI 


mnsz 


«.f*r  /P.EECHAR.XEE.DOP, D02P.00U.0O2N.DO3N.DE  .D0.D02D.SN. 


V02P, VON, V02N.VD3N, VO, V03, EXCITE ,VE.VO?D,VF .ALFA 


X (21 


X (3) 


X (A) 


THIS  PAGE  IS  BEST  QUALITY  PRACTICABU 
Kao*  COPY  yUEOtlSHKD  TO  DOC  «■"' 


Table  E - 


Chemistry  Program,  CHEMCO  (continued) 


pyof:v')»  (n; 


P<9»0'  »CPP 


»7r-'«-n  p*n 


>«r>fv<>()?p0OM 


P 1 S»Of 


g?)»r*0»n? 


P??<iFeqp<KVA_ 


q?'»oF?v*£P 


PP40P3- »0P 


o?A#n?v#o 


P?o*0?\*0?P 


tno«-oP»o*oM 


*0*0*0>* 


W3?*c-*u?*o? 


p 7 lononon? 


P3F*0?'ioO'-» 


pTQe.OP*OP 


1 op*r.?*  <o?*o?n 


P4  3*02f  *0?P 


RF3*E*0P 


W5**0M#F 


q<E,s*<1?*03 


P71  *'.)2*0?n 


‘.«2*02n*n?o 


P{i3*F*fc'»0P 


-V  k 

THIS  PACE  IS  BEST  QUALITY  PRACTICABLE 
FHOH  COPY  TURMISHED  TO  - 


Table  E - I Chemistry  Program,  CHEMCO  (continued) 


R94*E*E*o?p 


=p] oo*n?n 


= wi  04*03*0 


= R 1 ns*u?p*E 


= 07*020*0 


R104*o?fi*03 


?*n?*r 


Rn4*0M*02 


rn  i 7*00*020 


opo 


= Rl] 9*020*03 


= R 1 23*03*'*0 


• R)?4*020*0N 


= Rl  ?P*.)?P*E 


= R 1 30*O?n*E*O 


= Kl34*i)20*0*0 


= R136*0?I>*0P 


= R13R*O20*E 


= Rn9*0?0*E 


= c!140*O?0*E 


.=-  014  1 *021'*  <C2+0?0) *E 


RX*Oiv 


»Y*02N 


PZ*03\ 


liMHH 


0P+02R 


uaeok 


.2*00*2.24* 

)/<XP*SN* 

Inal 

00) 

19*  (3.2*00*2. 16*02N»2. 99*030)  /<XN*SN*1  . 


XPJ  = P.5H1: -2*/P-l 


A'lA  = FF  COAR*2t>  • 79 


AL^A  = XN/(E*1. E-100) 


AX  1 =-  (l,  *06%  a*2.*ALEA*GAMA) 


VE**  ( l . ♦ A|_r  «)  /XEK 


4X1  = ALM*x:i<-' /Xt.** 


*AX1*  n . ♦ t X3  j/T  11.*ALFa*GAMA)*(1.*AX2*AX3).) 


lii/n-wm  xviyftcT^iTj  JSEyi  ftn  Ti  fUdll 


XPP*i-.x)  /ll  ,*AX2*AX3) 


>1  no) 


SO***”** 


Uht  2h«OAW»0 o/  (Sfj*XI\M*l  .E-l  no) 


i.  P « ( ) A I ! * 0 2 N / < 5 N * X (»’  M ♦ 


n*0  A»m*03N/ 


r ♦'*02 


WraiSWKwa*  ,*-^  - 


Vita 


John  W.  Dettmer  was  born  on  27  September  1935  in  Cin- 
cinnati, Ohio.  He  was  graduated  from  Purcell  High  School  in 
1951  and  attended  the  University  of  Cincinnati,  receiving  an 
Electrical  Engineering  Degree  and  a commission  in  the  U.S. 

Air  Force  in  1956.  He  entered  active  duty  on  17  December 
1956  and  attended  pilot  training.  After  being  awarded  his 
wings,  he  flew  B-47Es  for  the  Strategic  Air  Command  from 
Davis-Monthan  AFB,  Arizona  until  December  1963.  He  then  at- 
tended the  air  Force  Institute  of  Technology  where  he  re- 
ceived a Master  of  Science  in  Electronics  Engineering  in 
1965.  He  then  served  as  a Communications  Engineer  with  the 
Ground  Electronics  Engineering  and  Installation  Agency  at 
Wheeler  AFB,  Hawaii  until  1966,  at  which  time  he  attended 
the  Aerospace  Test  Pilot  School  at  Edwards  AFB.  After  re- 
ceiving the  Rest  Pilot  specialty  in  1967,  he  was  assigned 
to  the  Test  Operations  Division  of  the  Manned  Orbiting  Labor- 
atory Program  located  at  Los  Angeles  A.F.S.  In  1968  he  was 
assigned  as  a Plans  Officer  with  the  Seventh  Air  Force  Head- 
quarters located  at  Tan  Son  Nhut,  RVN.  From  1969  until  1971 
he  served  as  a test  pilot  and  operations  officer  for  the 
Cargo  Test  Branch  of  the  4950  Test  Wing  at  Wright-Patterson 
AFB,  Ohio.  In  1972  he  again  entered  the  Air  Force  Institute 
of  Technology,  becoming  a Doctoral  Candidate  in  1973.  He 
was  then  assigned  to  Aerospace  Research  Laboratories  (ARL) 
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to  pursue  his  doctoral  research.  Upon  disestablishment  of 
the  organization  in  1975  until  August  1978,  he  was  again  as- 
signed to  the  4950  Test  Wing  as  the  Assistant  to  the  Comman- 
der for  the  Airborne  Laser  Laboratory  (ALL)  and  as  Commander 
of  the  Wing's  Detachment  at  Kirtland  AFB.  He  was  then  as- 
signed to  the  Weapons  Laboratory  as  Director  of  Test  and 
Operations  for  the  ALL.  During  the  period  after  the  dises- 
tablishment of  ARL,  he  pursued  his  doctoral  research  on  a 
part-time  basis. 

John  W.  Dettmer  is  a member  of  Eta  Kappa  Nu  and  the 
American  Institute  of  Electronic  and  Electrical  Engineers. 

He  is  married  to  the  former  Erin  Lee  Ross  of  Tucson,  Arizona; 
and  they  have  four  children. 


292 


Security  CLASSIFICATION  of  THIS  PAGE  (When  Ueta  Entered) 


l 


REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 

' REPORT  NUMBER  2.  GOVT  ACCESSION  NO. 

AFIT/DS/78-3  

3.  RECIPIENT'S  CATALOG  NUMBER 

- ■ 

A.  TITLE  (and  SublltU) 

Discharge  Processes  In  The  Oxygen  Plasma 

5.  TYPE  OF  REPORT  ft  PERIOD  COVERED 

f-;  idi 

6.  PERFORMING  ORG.  REPORT  NUMBER 

7.  AUTMORfiJ 

John  W.  Dettmer,  Colonel,  USAF 

Air  Force  Weapons  Laboratory 

Kirtland  AFB,  New  Mexico  87117 

8.  CONTRACT  OR  GRANT  NUMBER(») 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Air  Force  Institute  of  Technology 

Wright  Patterson  AFB,  Ohio  45433 

10.  PROGRAM  ELEMENT.  PROJECT.  TASK 

AREA  ft  WORK  UNIT  NUMBERS  ; 

N/A 

11  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Same  as  9. 

12.  REPORT  DATE 

November  1978  ' 

13.  NUMBER  OF  PAGES 

14.  MONITORING  AGENCY  NAME  ft  AODRESSf/f  different  from  Controlling  Oftice) 

Same  as  9. 

15.  SECURITY  CLASS,  (of  thla  report) 

Unclassified 

15a.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 

IS.  DISTRIBUTION  ST ATEMENT  (ol  thla  Report) 

Approved  for  public  release,  distribution  unlimited. 

17.  DISTRIBUTION  STATEMENT  (ol  the  abatract  entered  In  Block  20,  It  different  from  Report) 

10.  SUPPLEMENTARY  NOTES  APPftr 

Oxygen  1C'  RCLCASE  AFR  190  17. 

FCB  P 1977  jostoV.  nipps/^avr',  USAF* 

Director  of  Information 

19.  KEY  WORDS  ( Continue  on  reverae  aide  II  neceeaary  and  Identify  by  block  number) 

Oxygen  discharge,  glow  discharge,  unstable  discharge,  oxygen  plasma,  solutions 
to  Boltzmann  equation,  plasma  chemistry  model,  O2  metastable,  negative  Ion, 
two  step  ionization. 

20.  ABSTRACT  (Continue  on  reverae  aide  If  neceaaary  and  Identify  by  block  number) 

The  electric  glow  discharge  in  oxygen  was  experimentally  studied  and  mathemati- 
cally modeled  for  pressures  of  one  to  ten  torr  and  current  densities  to  35  mil- 
liamperes  per  square  centimeter.  Experimentally  it  was  observed  that  the  oxygen 
discharge  existed  in  either  a high-electric-field  or  a low-electric-field  form. 

The  high  field  form  was  present  at  high  currents  and  low  pressures,  the  low 
field  form  at  low  currents  and  high  pressures.  The  high  field  form  exhibited 
stable  characteristics  while  the  low  field  form  was  oscillatory.  Solutions  of 

DD  ,Er»  1473  EDITION  OF  I NOV  6S  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Wimi  0«f»  F.nfr.d) 


H 


SECURITY  CLASSIFICATION  OF  THIS  PAG6(H,han  Data  Entered) 


20.  ABSTRACT  (continued) 

the  Boltzmann  transport  equation  indicated  that  the  electron  energy  distribu- 
function  exhibited  a characteristic  intermediate  between  Maxwellian  and 
Auyvesteyn.  Chemistry  model  results  showed  the  dominant  positive  ions  to  be 
™ and  0".  The  transition  from  the  stable  to  unstable  forms  was  determined  to 
be  related  to  two  criteria:  the  dominance  of  the  derivative  of  the  attachment 
coefficient  with  respect  to  E/N  over  that  of  the  ionization  coefficient  and  the 
equivalence  of  the  negative  ion  and  electron  number  densities.  The  chemistry 
model  results  indicated  that  two-step  ionization  from  the  metastable  02(a'A) 
and  collisions  between  the  metastable  and  negative  ions  were  important  electron 
sources. 


SECURITY  CLASSIFICATION  OF  THIS  PAGEflWlFn  Dmlm  Entmrmd) 


